| | An Investigation of Cricoarytenoid Joint Mechanics Using Simulated Muscle ForcesAccepted 16 April 2002. Abstract Summary: Rotational and translational stiffnesses were calculated for arytenoid motion about the cricoarytenoid joint. These calculations were obtained from measurements on five excised human larynxes. For each larynx, known forces were applied to the arytenoid cartilage, and three markers were tracked as a function of applied forces. Assuming rigid body motion, arytenoid translations and rotations were computed for each applied force. Translational stiffnesses were obtained by plotting force versus displacement, and rotational stiffnesses were calculated by plotting torque versus angular rotation. A major finding was that the translational stiffness along the anterior-posterior direction was three times as great as the translational stiffnesses in the other two directions. This nonisotropic nature of the stiffnesses may be an important consideration for phonosurgeons who wish to avoid subluxation of the cricoarytenoid joint in patients. The computed rotational and translational stiffnesses currently are being implemented in 2D and 3D models. These stiffness parameters play a vital role in prephonatory glottal shaping, which in turn exerts a major influence on all aspects of vocal fold vibration, including fundamental frequency, voice quality, voice register, and phonation threshold pressure.
INTRODUCTION  The prephonatory shape of the glottis exerts a major influence on all aspects of vocal fold vibration, including fundamental frequency, voice quality, voice register, and phonation threshold pressure. Several key factors that influence glottal shape include the mechanical properties of the laryngeal tissues, muscular contractile forces, and the stiffness of several joints. The present study is an investigation of cricoarytenoid joint (CAJ) mechanics. Specifically, the mechanics of the CAJ are studied as a function of simulated muscle forces. Many prior investigations on excised larynxes have documented the motion of the arytenoids cartilage about the CAJ.1., 2., 3., 4., 5., 6., 7., 8., 9. Fortunately, most of these studies report similar findings. They agree that the principal motion of the arytenoid cartilage is a rocking about a roughly longitudinal axis of the cricoid cartilage, as well as a sliding motion along this axis. With a few notable exceptions, most of these studies have been qualitative in nature. Furthermore, all of these previous studies have been exclusively kinematic, i.e., CAJ movements were studied without regard for the forces generating such movements. Tomographic observations4 on living subjects have also confirmed that the arytenoid rocks on the cricoid facet during vocal fold movement. Recently, Selbie et al.10 conducted a quantitative study on the geometry of the CAJ facet surfaces. Based on these geometries, they also reported that the optimal axis for movement on these surfaces corresponded to a rocking of the arytenoid on the cricoid. The present investigation is an important extension of former CAJ studies because it quantifies and analyzes CAJ movement as a function of known muscle forces. With this information, effective rotational and translational stiffnesses are computed for the CAJ, which are critical for modeling the prephonatory shape of the vocal folds. For example, ongoing studies have begun to incorporate this information into 2D and 3D biomechanical models of CAJ motion.11 Thus, effective rotational and translational stiffnesses are computed for both 2D and 3D models of the CAJ. In the analysis, rigid body motion is assumed. Although rotation and translation of the arytenoid is allowed, any deformation of the arytenoid cartilage is assumed to be negligible (within its physiological range of motion). This was also the assumption of the Selbie et al10 study, which used CAJ facet geometry to infer motion of the arytenoid cartilage. METHODS Five male larynges were obtained from the autopsy unit at the University of Iowa hospitals and clinics. At the time of death, the males ranged in age from 44 to 82 years, and in weight from 110 to 237 lbs. Specific data for each subject are shown in Table 1. No history of voice disorders was known for any subject and no abnormalities were apparent on any larynx. | | |  | | Age (years) | Weight (Ibs) | |
|---|
| Larynx 1 | 44 | 110 | |
| Larynx 2 | 56 | 198 | |
| Larynx 3 | 82 | 165 | |
| Larynx 4 | 67 | 196 | |
| Larynx 5 | 56 | 237 | |
| Mean | 61.10 | 181.2 | |
| SD | 14.3 | 47.3 | | | | |
Immediately after harvest, each larynx was refrigerated and stored in saline solution. Experiments were performed within 5–7 days of harvest. Each larynx was dissected and mounted. Initially, the trachea was shortened to approximately 2 inches. All extrinsic laryngeal muscles were removed, and the superior portion of the thyroid cartilage was removed to the level of the ventricular folds. During the final phase of dissection, the ventricular folds were removed to allow an unobstructed view of the true vocal folds. Throughout the dissection process, care was taken to leave as much laryngeal tissue intact as possible. However, some extraneous tissue surrounding the muscular process had to be removed to place sutures, and simulate the action of laryngeal muscles. Mounting consisted of fixing both the cricoid and thyroid cartilages. The cricoid cartilage was fixed by placing a rigid plastic tube through the opening of the trachea and into the subglottal space. Screws were placed through the cricoid cartilage to secure the tube. The hard plastic tube was clamped in order to hold the larynx in a vertical position. The thyroid cartilage was sutured to a fixed crossbar, thus fixing its position as well. Through the experiment, vocal fold tissues were periodically moistened with saline solution to minimize stiffness changes caused by dehydration. In order to document the motion of the arytenoid as a function of simulated muscle forces, three microsutures were placed on the arytenoid as markers. Specifically, a microsuture was placed on each of the following locations: (1) the vocal process, (2) the muscular process, and (3) the apex, as shown in Figure 1. Because the arytenoid was assumed to move as a rigid body, a minimum of three points was needed to define its motion. The positions of the three microsutures were measured before and after application of each of the simulated muscle forces using a 3D digitizer, the MicroScribe-3DX, with a reported resolution of 0.13 mm, from Immersion Corporation (801 Fox Lane, San Jose, California 95131 USA; Telephone: 408-467-1900). As indicated in Figure 1, the orientation of axes assumed throughout this study was as follows: the z-axis denoted the inferior-superior direction, the x-axis denoted the medial-lateral direction, and the y-axis denoted the anterior-posterior direction. In keeping with standard conventions used in the physical sciences, counter-clockwise rotations about the coordinate axes were positive, and clockwise rotations were negative. Six laryngeal muscle forces were simulated, including two components of the posterior cricoarytenoid (PCA) muscle, with separate horizontal and oblique compartments; the lateral cricoarytenoid (LCA) muscle; two components of the thyroarytenoid (TA) muscle, with separate vocalis and thyrovocalis compartments; and the interarytenoid (IA) muscle. The forces were simulated by attaching one end of a suture to the arytenoid (roughly corresponding to the muscle insertion location), coursing the suture in the direction of the muscle fibers, and attaching the remaining end of the suture to a weight and pulley system. The average directions of these applied forces (averaged over all five larynges) are shown in Table 2. Figure 1 illustrates how a suture was attached to the arytenoid cartilage to simulate the horizontal compartment of the PCA. The “zero” position (or the rest position) was defined to be the position at which no external forces (weights) were applied. | | | | Muscle | nx | ny | nz | |
|---|
| PCA horizontal | −0.696 | −0.665 | 0.270 | |
| PCA oblique | −0.396 | −0.657 | −0.642 | |
| LCA | 0.026 | 0.965 | −0.261 | |
| TA vocalis | −0.086 | 0.995 | 0.057 | |
| TA thyrovocalis | 0.116 | 0.985 | −0.128 | |
| IA | −0.921 | −0.383 | 0.068 | | | | |
The applied forces roughly approximated the magnitude and direction of physiological musculature forces on the arytenoid cartilage, as judged by the maximum active stress in laryngeal muscles (∼100 kPa) multiplied by the cross-sectional areas.12., 13. In a separate study, precise directions of muscle fiber bundles were quantified.14 Because it was the goal of this study to calculate effective rotational and translational stiffnesses of the CAJ, a variety of muscle forces were applied. The applied forces and resultant displacements were decomposed into rotational and translational components, enabling the calculation of rotational and translational stiffnesses. Through application of a large variety of muscular forces, substantial data were collected from which to infer the stiffnesses. In particular, for each of the 6 muscles, 4–5 force/weight conditions were applied, ranging from zero to a maximum weight, as shown in Table 3. | | | | Muscle | Weight (g) | |
|---|
| PCA horizontal | 50 | |
| PCA Oblique | 50 | |
| LCA | 50 | |
| TA vocalis | 100 | |
| TA thyrovocalis | 80 | |
| IA | 20 | | | | |
THEORY In our analysis, the description of arytenoid motion was performed with rigid body analysis. As noted previously, most of the studies of CAJ motion have indicated that the arytenoid rocks or slides on the cricoid. Rocking is a specific type of rigid body movement which, qualitatively, might be described as a rotation about a moving axis. Mathematically, rocking is usually implemented as an equation of constraint. In the following empirical analysis, although rocking or sliding may be the principal type of motion observed, no attempt is made to a priori constrain the arytenoid cartilage to such motion. Assuming the arytenoid to move as a rigid body, the current position coordinates x′ of any point i on the arytenoid may be expressed as a rotation R of their initial positions xi (about the origin) followed by a translation vector D:
If, instead, one desires the rotation to occur about the center of mass of the arytenoid, the equation may be written as:
where xcm represents the center of mass of the arytenoid. For a 2D biomechanical model, the rotation matrix R represents a rotation θ z about the z axis and is written as:
In this 2D model, x′i, xi, xcm, and D are two-dimensional vectors with x and y components. In a more general three-dimensional model, these vectors are three dimensional and the total rotation R is considered to be three independent rotations about the three coordinate axes, i.e., R =RxRyRz, where
In this experiment, for each force applied to the arytenoid, the initial positions xi and the final positions x′i were measured at the three marked locations on the arytenoid (the vocal process, the muscular process, and the apex). The center of mass of the arytenoid xcm was also estimated, although any arbitrarily chosen reference point would have been equally valid for the present study. Thus, for each applied force, only the rotation matrix R and the translation vector D were unknown. By defining x̄′ito be the average value of the position coordinates of the three arytenoid markers after an applied force, and x̄i to be the average value of the three arytenoid markers before the applied force, the two-dimensional translation vector D could be expressed, for any angle θz as:
By computing the translation vector D in this way, 15 Equation 2 could be minimized with respect to just one variable, θ z. Specifically, using the FMIN routine in MATLAB 5.3 Release 11, the mean square difference between the left-hand side and right-hand side of Equation 2 was minimized as a function of θ z. For the three-dimensional problem, Equation 5 was also used to calculate the translation vector D. However, in this case, Equation 2 required minimization with respect to three variables: θ x, θ y, and θ z. This was done using the FMINS routine in MATLAB (similar to the FMIN routine, except that the FMINS routine allows optimization over more than one variable). Each force applied to the arytenoid brought it to a new state of static equilibrium, which was quantified by measuring the new position coordinates for the three arytenoid markers. For each state, the sum of all external forces on the arytenoid was set to zero, and the sum of all external torques or moments on the arytenoid was set to zero. For the two-dimensional problem, the CAJ attachment to the arytenoid was modeled as three separate components, including: (1) a linear spring in the x-direction, (2) a linear spring in the y-direction, and (3) a linear torsional spring about the z-axis. Thus, the equations of equilibrium were written as:
Similarly, for the three-dimensional problem, the CAJ attachment to the arytenoid was modeled as six separate components, with a linear spring and a linear torsional spring associated with each coordinate axis. In particular, the equilibrium equations were written as:
The vector rattached/cm was measured from the arytenoid center of mass to the suture attachment of the applied force to the arytenoid. Fapplied was the applied force vector and Fapplied,x, Fapplied,y, Fapplied,z, the x, y, and z components of the force vector; kx, ky, and kz were the translational stiffnesses along the x, y, and z axes; and κθx, κθy, and κθx were the torsional stiffnesses about the x, y, and z axes. The magnitude of the total force was taken as the weight attached to the mass/pulley system. For each applied force, the direction of the force was determined by digitizing two locations: (1) the suture attachment of the applied force to the arytenoid, and a location on the suture about 1 cm away from the attachment. A unit vector capturing the direction of each force was determined by subtracting the first point from the second, and then normalizing. One might question whether it was necessary to conduct separate optimizations for both 2D and 3D models. Couldn't the 2D model be considered merely a projection of the 3D? In theory, perhaps yes. However, in practice, the optimization problem was different in 2D than in 3D. In particular, the z data were ignored in the 2D optimization. Thus, because the 2D and 3D models optimized over different data sets, the 2D model was not merely a projection of the 3D model. RESULTS AND DISCUSSION After decomposing the forces and displacements into angular and translational components, plots of force (or torque) versus displacement (or angular rotation) were made for each of the equations of Equation 6. In the 2D analysis, the translational (kx, ky) and rotational (κθx) stiffnesses were then computed by determining the slope of the best-fit line through the data with a least-square error procedure, as shown in FIGURE 2, FIGURE 3, FIGURE 4 for larynx 5. FIGURE 2, FIGURE 3, FIGURE 4 show the data from all the larynxes on one plot. Although the data may appear noisy when the data from all the subjects are plotted together, in general, this is not the case when data from a single subject are used, as illustrated in FIGURE 2, FIGURE 3, FIGURE 4. Thus, FIGURE 2, FIGURE 3, FIGURE 4 suggest a relatively large intersubject variability, which is often a common feature of biomechanical data. To avoid pooling incomparable data, translational and rotational stiffness were computed, separately, for each of the five larynxes. To explicitly examine inter-subject variability, mean values and standard deviations were also computed, as summarized in Table 4. The slopes drawn in FIGURE 2, FIGURE 3, FIGURE 4 merely illustrate the mean slopes computed in Table 5. | | | | | kx (N/m) | ky (N/m) | kθz (Nm/rad) | |
|---|
| Larynx 1 | 93 | 250 | 0.0117 | |
| Larynx 2 | 86 | 225 | 0.0056 | |
| Larynx 3 | 93 | 311 | 0.0103 | |
| Larynx 4 | 82 | 331 | 0.0026 | |
| Larynx 5 | 57 | 146 | 0.0046 | |
| Mean | 82.2 | 252.6 | 0.0070 | |
| SD | 14.9 | 73.6 | 0.0039 | | | | |
| | | | | kx (N/m) | ky (N/m) | kz (N/m) | kθx (Nm/rad) | kθy (Nm/rad) | kθz (Nm/rad) | |
|---|
| Larynx 1 | 75 | 267 | 94 | 0.0044 | 0.0039 | 0.0112 | |
| Larynx 2 | 90 | 217 | 132 | 0.0078 | 0.0031 | 0.0054 | |
| Larynx 3 | 92 | 302 | 50 | 0.0042 | 0.0034 | 0.0102 | |
| Larynx 4 | 41 | 156 | 46 | 0.0038 | 0.0055 | 0.0024 | |
| Larynx 5 | 64 | 132 | 34 | 0.0018 | 0.0022 | 0.0049 | |
| Mean | 72.4 | 214.8 | 71.2 | 0.0044 | 0.0036 | 0.0068 | |
| SD | 21.0 | 71.8 | 40.9 | 0.0022 | 0.0012 | 0.0037 | | | | |
In these data, perhaps the most noticeable feature is that ky, the stiffness along the y-axis (the anterior-posterior direction), is roughly three times kx, the stiffness along the x axis (the medial-lateral direction). This would seem to imply that, while some arytenoid translation may occur in connection with arytenoid adduction, relatively little arytenoid translation occurs in connection with lengthening and shortening of the folds. Alternately, it could simply mean that larger anterior-posterior force components are used to generate arytenoid translation along the anterior-posterior length of the folds, about three times as large as in the medial-lateral direction. In either case, our data suggest that medial-lateral gliding is less constrained than anterior-posterior gliding. In some cases, exponential or nonlinear fits to the data might be more appropriate than a linear fit. However, in most cases, the linear fit did capture most of the variance of the data. Thus, the biomechanical constants provided a reasonable initial description of arytenoid motion about the CAJ. Similarly, rotational and translational stiffnesses from the 3D analysis were calculated from the slopes of the force-displacements/torque-angular displacement plots in FIGURE 5, FIGURE 6, FIGURE 7, FIGURE 8, FIGURE 9, FIGURE 10. Again, composite data are shown in FIGURE 5, FIGURE 6, FIGURE 7, FIGURE 8, FIGURE 9, FIGURE 10, and an example of single subject data is shown FIGURE 5, FIGURE 6, FIGURE 7, FIGURE 8, FIGURE 9, FIGURE 10. The translational and rotational stiffnesses of each of the five larynxes, as well as mean values and standard deviations, are summarized in Table 5. In general, the three-dimensional model tells a similar story to the two-dimensional model. Inparticular, the values of kx, ky, and κθz are similar for both models, although not identical for the reason mentioned earlier (e.g., z data not included in the 2D optimization). Again, the anterior-posterior translational stiffness, ky, is about three times greater than the stiffnesses in the other two directions. In terms of rotational stiffnesses, the stiffness about the z axis, κθz is approximately 55%–90% greater than the stiffnesses about the other two axes. SUMMARY Rotational and transitional stiffnesses of the CAJ were calculated for both 2D and 3D models of the CAJ. To calculate these stiffnesses, arytenoid motion was studied as a function of applied forces. Although many previous kinematic studies of arytenoid motion have been conducted, this is the first investigation to study arytenoid motion as a function of applied forces. A major finding was that the translational stiffness along the anterior-posterior direction was roughly three times larger than the translational stiffnesses in the other two directions, suggesting that the primary constraint of the ligaments aroundthe CAJ is to keep the arytenoid firmly attached to the cricoid during vocal fold elongation. Medial-lateral and vertical gliding was less constrained. This nonisotropic nature of the stiffnesses may be an important consideration for phonosurgeons who wish to avoid subluxation of the CAJ in patients. These rotational and translational stiffnesses are critical for modeling the prephonatory posturing of the vocal folds, which exert a major influence on many aspects of vocal fold vibration. Preliminary investigations yielded reasonable results when these CAJ stiffness values were inserted into our models of vocal fold posturing.11 In particular, laryngeal muscles were able to both adduct and abduct the vocal folds. We note that the effective stiffnesses calculated in this investigation were computed based on a linear fit to the data. However, to implement the data with finite element models, it may be necessary to use higher-order (strain-dependent) stiffnesses. This would allow for less constraint at small displacements and greater constraint at large displacements, a condition often encountered in biomechanics to avoid injury. Acknowledgements This work was supported by grants R01 DC03072 and P60 DC00976 and from the National Institute on Deafness and Other Communication Disorders. The preliminary investigations of Dr. Mike Bayerl during his postdoctoral research fellowship at the University of Iowa were also helpful in the present study. References 1..
1.
Sonesson B
.
Die funktionelle anatomie des cricoarytenoidgelenkes
.
Zeitschrift für Anatomie und Entwicklungsgeschichte
. 1959;121:292–303
.
MEDLINE 2..
2.
Frable MA
.
Computation of motion at the cricoarytenoid joint
.
Arch Otolaryngol
. 1961;73:551–556
.
MEDLINE 3..
3.
von Leden H
, Moore P
.
The mechanics of the cricoarytenoid joint
.
Arch Otolaryngol
. 1961;73:541–550
.
MEDLINE 4..
4.
Ardran GM
, Kemp FH
.
The mechanisms of the larynx, part I: the movements of the arytenoid and cricoid cartilages
.
Br J Radiol
. 1966;39:641–654
.
MEDLINE |
CrossRef
5..
5.
Broad DJ
.
Kinematic considerations for evaluating laryngeal cartilage motions
.
Folia Phoniatr
. 1968;20:269–284
.
6..
6.
Maue WM
, Dickson DR
.
Cartilages and ligaments of the adult human larynx
.
Arch Otolaryngol
. 1971;94:432–439
.
MEDLINE 7..
7.
Sellars IE
, Vaughn CL
.
A biomechanical investigation of cricoarytenoid joint kinematics
.
International Series on Biomechanics
.
Biomechanics VIII-A
. Vol. 4A.
Champaign, IL: Human Kinetics Publishers; 1981;
.
8..
8.
Sellars IE
, Sellars S
.
Cricoarytenoid joint structure and function
.
J Laryngol Otol
. 1983;97:1027–1034
.
MEDLINE 9..
9.
Neuman TR
, Hengesteg A
, Lepage RP
, Kaufman KR
, Woodson GE
.
Three-dimensional motion of the arytenoid adduction procedure in cadaver larynges
.
Ann Otol Rhinol Laryngol
. 1994;103:265–270
.
MEDLINE 10..
10.
Selbie SW
, Zhang L
, Levine WS
, Ludlow CL
.
Using joint geometry to determine the motion of the cricoarytenoid joint
.
J Acoust Soc Am
. 1998;103:1115–1127
.
MEDLINE |
CrossRef
11..
11.
Hunter E. Three Dimensional Biomechanical Model of Vocal Fold Posturing. [PhD thesis]. Iowa City, IA: The University of Iowa; 2001.
12..
12.
Cox KA
, Alipour F
, Titze IR
.
Geometric structure of the human and canine cricothyroid and thyroarytenoid muscles for biomechanical applications
.
Ann Otol Rhinol Laryngol
. 1999;108:1151–1158
.
MEDLINE 13..
13.
Tayama N
, Chan RW
, Kaga K
, Titze IR
.
Functional definitions of vocal fold geometry for laryngeal biomechanical modeling
.
Ann Otol Rhinol Laryngol
. 2002;111:83–92
.
MEDLINE 14..
14.
Mineck CW
, Tayama N
, Chan RW
, Titze IR
.
Three-dimensional anatomic characterization of the canine laryngeal abductor and adductor musculature
.
Ann Otol Rhinol Laryngol
. 2001;109:505–513
.
MEDLINE 15..
15.
Spoor CW
, Veldpaus FE
.
Rigid body motion calculated from spatial coordinates of markers
.
J Biomechanics
. 1980;13:391–393
.
∗ Department of Speech Pathology and Audiology, National Center for Voice and Speech, The University of Iowa, Iowa City, Iowa, USA † Department of Otolaryngology, The University of Iowa, Iowa City, Iowa, USA  Address correspondence and reprint requests to David A. Berry, PhD, UCLA School of Medicine, Division of Head and Neck Surgery, 31-24 Rehabilitation Center, 1000 Veteran Avenue, Los Angeles, CA 90095-1794, USA
PII: S0892-1997(03)00026-2 doi:10.1016/S0892-1997(03)00026-2 © 2003 The Voice Foundation. Published by Elsevier Inc. All rights reserved. | |
|
FULL TEXT00026-2/journalimage?img=PIIS0892199703000262.gr1.lrg.gif&fig=FIG1&kwhquery=&issn=0892-1997&ishighres=false&allhighres=false&free=yes','journalimage');)
00026-2/journalimage?img=PIIS0892199703000262.gr2.lrg.gif&fig=FIG2&kwhquery=&issn=0892-1997&ishighres=false&allhighres=false&free=yes','journalimage');)
00026-2/journalimage?img=PIIS0892199703000262.gr3.lrg.gif&fig=FIG3&kwhquery=&issn=0892-1997&ishighres=false&allhighres=false&free=yes','journalimage');)
00026-2/journalimage?img=PIIS0892199703000262.gr4.lrg.gif&fig=FIG4&kwhquery=&issn=0892-1997&ishighres=false&allhighres=false&free=yes','journalimage');)
00026-2/journalimage?img=PIIS0892199703000262.gr5.lrg.gif&fig=FIG5&kwhquery=&issn=0892-1997&ishighres=false&allhighres=false&free=yes','journalimage');)
00026-2/journalimage?img=PIIS0892199703000262.gr6.lrg.gif&fig=FIG6&kwhquery=&issn=0892-1997&ishighres=false&allhighres=false&free=yes','journalimage');)
00026-2/journalimage?img=PIIS0892199703000262.gr7.lrg.gif&fig=FIG7&kwhquery=&issn=0892-1997&ishighres=false&allhighres=false&free=yes','journalimage');)
00026-2/journalimage?img=PIIS0892199703000262.gr8.lrg.gif&fig=FIG8&kwhquery=&issn=0892-1997&ishighres=false&allhighres=false&free=yes','journalimage');)
00026-2/journalimage?img=PIIS0892199703000262.gr9.lrg.gif&fig=FIG9&kwhquery=&issn=0892-1997&ishighres=false&allhighres=false&free=yes','journalimage');)
00026-2/journalimage?img=PIIS0892199703000262.gr10.lrg.gif&fig=FIG10&kwhquery=&issn=0892-1997&ishighres=false&allhighres=false&free=yes','journalimage');)