Phonation threshold pressure and onset frequency in a two-layer physical model of the vocal folds

The influence of vocal fold geometry and stiffness on phonation onset was experimentally investigated using a body-cover physical model of the vocal folds. Results showed that a lower phonation threshold pressure and phonation onset frequency can be achieved by reducing body-layer or cover-layer stiffness, reducing medial surface thickness, or increasing cover-layer depth. Increasing body-layer stiffness also restricted vocal fold motion to the cover layer and reduced prephonatory glottal opening. Excitation of anterior-posterior modes was also observed, particularly for large values of the body-cover stiffness ratio. The results of this study were also discussed in relation to previous theoretical and experimental studies.     

Physical mechanisms of phonation onset: a linear stability analysis of an aeroelastic continuum model of phonation

In an investigation of phonation onset, a linear stability analysis was performed on a two-dimensional, aeroelastic, continuum model of phonation. The model consisted of a vocal fold-shaped constriction situated in a rigid pipe coupled to a potential flow which separated at the superior edge of the vocal fold. The vocal fold constriction was modeled as a plane-strain linear elastic layer. The dominant eigenvalues and eigenmodes of the fluid-structure-interaction system were investigated as a function of glottal airflow. To investigate specific aerodynamic mechanisms of phonation onset, individual components of the glottal airflow (e.g., flow-induced stiffness, inertia, and damping) were systematically added to the driving force. The investigations suggested that flow-induced stiffness was the primary mechanism of phonation onset, involving the synchronization of two structural eigenmodes. Only under conditions of negligible structural damping and a restricted set of vocal fold geometries did flow-induced damping become the primary mechanism of phonation onset. However, for moderate to high structural damping and a more generalized set of vocal fold geometries, flow-induced stiffness remained the primary mechanism of phonation onset.     

Simulation of vocal fold impact pressures with a self-oscillating finite-element model

Vocal fold impact pressures were studied using a self-oscillating finite-element model capable of simulating vocal fold vibration and airflow. The calculated airflow pressure is applied on the vocal fold as the driving force. The airflow region is then adjusted according to the calculated vocal fold displacement. The interaction between airflow and the vocal folds produces a self-oscillating solution. Lung pressures between 0.2 and 2.5 kPa were used to drive this self-oscillating model. The spatial distribution of the impact pressure was studied. Studies revealed that the tissue collision during phonation produces a very large impact pressure which correlates with the lung pressure and glottal width. Larger lung pressure and a narrower glottal width increase the impact pressure. The impact pressure was found to be roughly the square root of lung pressure. In the inferior-superior direction, the maximum impact pressure is related to the narrowest glottis. In the anterior-posteriorfirection, the greatest impact pressure appears at the midpoint of the vocal fold. The match between our numerical simulations and clinical observations suggests that this self-oscillating finite-element model might be valuable for predicting mechanical trauma of the vocal folds.     

Spatiotemporal classification of vocal fold dynamics by a multimass model comprising time-dependent parameters

A model-based approach is proposed to objectively measure and classify vocal fold vibrations by left-right asymmetries along the anterior-posterior direction, especially in the case of nonstationary phonation. For this purpose, vocal fold dynamics are recorded in real time with a digital high-speed camera during phonation of sustained vowels as well as pitch raises. The dynamics of a multimass model with time-dependent parameters are matched to vocal fold vibrations extracted at dorsal, medial, and ventral positions by an automatic optimization procedure. The block-based optimization accounts for nonstationary vibrations and compares the vocal fold and model dynamics by wavelet coefficients. The optimization is verified with synthetically generated data sets and is applied to 40 clinical high-speed recordings comprising normal and pathological voice subjects. The resulting model parameters allow an intuitive visual assessment of vocal fold instabilities within an asymmetry diagram and are applicable to an objective quantification of asymmetries.     

Visualization and Quantification of the Medial Surface Dynamics of an Excised Human Vocal Fold During Phonation

SUMMARY: The purpose of this investigation was to investigate physical mechanisms of vocal fold vibration during normal phonation through quantification of the medial surface dynamics of the fold. An excised hemilarynx setup was used. The dynamics of 30 microsutures mounted on the medial surface of a human vocal fold were analyzed across 18 phonatory conditions. The vibrations were recorded with a digital high-speed camera at a frequency of 4,000 Hz. The positions of the sutures were extracted and converted to three-dimensional coordinates using a linear approximation technique. The data were reduced to principal eigenfuctions, which captured over 90% of the variance of the data, and suggested mechanisms of sustained vocal fold oscillation. The vibrations were imaged as the following phonatory conditions were manipulated: glottal airflow, an adductory force applied to the muscular process, and an elongation force applied to the thyroid cartilage. Over the range of variables studied, only the variation in glottal airflow yielded significant changes in subglottal pressure and fundamental frequency. All recordings showed high correlation for the distribution of the dynamics across the medial surface of the vocal fold. The distribution of the different displacement directions and velocities showed the highest variations around the superior region of the medial surface. Although the computed vibration patterns of the two largest empirical eigenfunctions were consistent with previous experimental observations, the relative prominence of the two eigenfunctions changed as a function of glottal airflow, impacting theories of vocal efficiency and vocal economy.     

The effects of rehydration on phonation in excised canine larynges

Experiments using excised canine larynges were conducted to study the restoration of vocal efficiency in dehydrated larynges. Excised larynges were dehydrated with warm, dry air to the point that airflow through the approximated vocal folds would not entrain the folds to produce phonation. The dehydrated vocal folds were then bathed in a saline solution. The rehydrated larynges were then remounted on the bench apparatus that enabled phonation with a constant humidified airflow, and measurements were made of phonation threshold pressure, glottal airflow, and amplitude. Hydration resulted in significantly increased efficiency and decrease in phonation threshold pressure. The findings confirm clinical impressions that hydration is critical in the physiology of normal phonation.     

Vibration in a self-oscillating vocal fold model with left-right asymmetry in body-layer stiffness

This study compares the phonatory behavior of an asymmetric vocal fold model to that of each individual vocal fold model in a hemi-configuration. Although phonation frequencies of the two folds in hemi-configurations had a ratio close to 1:3, a subharmonic synchronization between the two folds was not observed in the asymmetric model. Instead, the vibratory behavior was dominated by the dynamics of one fold only, and the other fold was enslaved to vibrate at the same frequency. Increasing subglottal pressure induced a shift in relative dominance between the two folds, leading to abrupt changes in both vibratory pattern and frequency.     

A lumped mucosal wave model of the vocal folds revisited: recent extensions and oscillation hysteresis

This paper examines an updated version of a lumped mucosal wave model of the vocal fold oscillation during phonation. Threshold values of the subglottal pressure and the mean (DC) glottal airflow for the oscillation onset are determined. Depending on the nonlinear characteristics of the model, an oscillation hysteresis phenomenon may occur, with different values for the oscillation onset and offset threshold. The threshold values depend on the oscillation frequency, but the occurrence of the hysteresis is independent of it. The results are tested against pressure data collected from a mechanical replica of the vocal folds, and oral airflow data collected from speakers producing intervocalic /h/. In the human speech data, observed differences between voice onset and offset may be attributed to variations in voice pitch, with a very small or inexistent hysteresis phenomenon.     

Model-based classification of nonstationary vocal fold vibrations

Classification of vocal fold vibrations is an essential task of the objective assessment of voice disorders. For historical reasons, the conventional clinical examination of vocal fold vibrations is done during stationary, sustained phonation. However, the conclusions drawn from a stationary phonation are restricted to the observed steady-state vocal fold vibrations and cannot be generalized to voice mechanisms during running speech. This study addresses the approach of classifying real-time recordings of vocal fold oscillations during a nonstationary phonation paradigm in the form of a pitch raise. The classification is based on asymmetry measures derived from a time-dependent biomechanical two-mass model of the vocal folds which is adapted to observed vocal fold motion curves with an optimization procedure. After verification of the algorithm performance the method was applied to clinical problems. Recordings of ten subjects with normal voice and ten dysphonic subjects have been evaluated during stationary as well as nonstationary phonation. In the case of nonstationary phonation the model-based classification into "normal" and "dysphonic" succeeds in all cases, while it fails in the case of sustained phonation. The nonstationary vocal fold vibrations contain additional information about vocal fold irregularities, which are needed for an objective interpretation and classification of voice disorders.     

Assessment of local vocal fold deformation characteristics in an in vitro static tensile test

Voice quality is strongly dependent on vocal fold dynamics, which in turn are dependent on lung pressure and vocal fold biomechanics. Numerical and physical models are often used to investigate the interactions of these different subsystems. However, the utility of numerical and physical models is limited unless appropriately validated with data from physiological models. Hence a method that enables analysis of local vocal fold deformations along the entire surface is presented. In static tensile tests, forces are applied to distinctive working points being located in cover and muscle, respectively, so that specific layer properties can be investigated. The forces are directed vertically upward and are applied along or above the vocal fold edge. The resulting deformations are analyzed using multiple perspectives and three-dimensional reconstruction. Deformation characteristics of four human vocal folds were investigated. Preliminary results showed two phases of deformation: a range with a small slope for small deformations fading into a significant nonlinear deformation trend with a high slope. An increase of tissue stiffness from posterior to anterior was detected. This trend is more significant for muscle and in the mid-anterior half of the vocal fold.     

Resonance properties of the vocal folds: in vivo laryngoscopic investigation of the externally excited laryngeal vibrations

The study presents the first attempt to investigate resonance properties of the living vocal folds by means of laryngoscopy. Laryngeal vibrations were excited via a shaker placed on the neck of a male subject and observed by means of videostroboscopy and videokymography (VKG). When the vocal folds were tuned to the phonation frequency of 110 Hz and sinusoidal vibration with sweeping frequency (in the range 50-400 Hz) was delivered to the larynx, three clearly pronounced resonance peaks at frequencies around 110, 170, and 240 Hz were identified in the vocal fold tissues. Different modes of vibration of the vocal folds, observed as distinct lateral-medial oscillations with one, two, and three half-wavelengths along the glottal length, respectively, were associated with these resonance frequencies. At the external excitation frequencies below 100 Hz, vibrations of the ventricular folds, aryepiglottic folds and arytenoid cartilages were dominant in the larynx.     

Computation of physiological human vocal fold parameters by mathematical optimization of a biomechanical model

With the use of an endoscopic, high-speed camera, vocal fold dynamics may be observed clinically during phonation. However, observation and subjective judgment alone may be insufficient for clinical diagnosis and documentation of improved vocal function, especially when the laryngeal disease lacks any clear morphological presentation. In this study, biomechanical parameters of the vocal folds are computed by adjusting the corresponding parameters of a three-dimensional model until the dynamics of both systems are similar. First, a mathematical optimization method is presented. Next, model parameters (such as pressure, tension and masses) are adjusted to reproduce vocal fold dynamics, and the deduced parameters are physiologically interpreted. Various combinations of global and local optimization techniques are attempted. Evaluation of the optimization procedure is performed using 50 synthetically generated data sets. The results show sufficient reliability, including 0.07 normalized error, 96% correlation, and 91% accuracy. The technique is also demonstrated on data from human hemilarynx experiments, in which a low normalized error (0.16) and high correlation (84%) values were achieved. In the future, this technique may be applied to clinical high-speed images, yielding objective measures with which to document improved vocal function of patients with voice disorders.     

Restraining mechanisms in regulating glottal closure during phonation

Recent experimental studies showed that isotropic vocal fold models were often blown wide apart and thus not able to maintain adductory position, resulting in voice production with noticeable breathy quality. This study showed that the capability of the vocal fold to resist deformation against airflow and maintain adductory position can be improved by stiffening the body-layer stiffness or increasing the anterior-posterior tension of the vocal folds, which presumably can be achieved through the contraction of the thyroarytenoid (TA) and cricothyroid (CT) muscles, respectively. Experiments in both physical models and excised larynges showed that, when these restraining mechanisms were activated, the vocal folds were better able to maintain effective adduction, resulting in voice production with much clearer quality and reduced breathiness. In humans, one or more restraining mechanisms may be activated at different levels to accommodate the varying degree of restraining required under different voice conditions.     

Phonation Threshold Power in Ex Vivo Laryngeal Models

This study hypothesized that phonation threshold power is measureable and sensitive to changes in the biomechanical properties of the vocal folds. Phonation threshold power was measured in three sample populations of 10 excised canine larynges treated with variable posterior glottal gap, variable bilateral vocal fold elongation, and variable vocal fold lesioning. Posterior glottal gap varied from 0 to 4 mm in 0.5 mm intervals. Bilateral vocal fold elongation varied from 0% to 20% in 5% intervals. Vocal fold lesion treatments included unilateral and bilateral vocal fold lesion groups. Each treatment was investigated independently in a sample population of 10 excised canine larynges. Linear regression analysis indicated that phonation threshold power was sensitive to posterior glottal gap (R² = 0.298, P < 0.001) and weakly to vocal fold elongation (R² = 0.052, P = 0.003). A one-way repeated measures analysis of variance indicated that phonation threshold power was sensitive to the presence of lesions (P < 0.001). Theoretical and experimental evidence presented here suggests that phonation threshold power could be used as a broad screening parameter sensitive to certain changes in the biomechanical properties of the larynx. It has not yet been measured in humans, but because it has the potential to represent the airflow-tissue energy transfer more completely than the phonation threshold pressure or flow alone, it may be a more useful parameter than these and could be used to indicate that laryngeal health is likely abnormal.     

The effect of head position on glottic closure in patients with unilateral vocal fold paralysis

The present study was designed to assess the effect of head position on glottic closure as reflected in airflow rates (open quotient and maximum flow declination rate), in patients with unilateral vocal fold paralysis. Ten patients, 2 males and 8 females ranging in age from 40 to 75, with a mean age of 57.3, served as subjects. Airflow measures were taken during sustained phonation of two vowels (/i/ and /a/) in 3 head positions (center, right, left). Vowels /i/ and /a/ were produced at subject's comfortable pitch and loudness, with random ordering of both vowel order and head orientation. Subjects were trained to focus eye gaze on right and left markers (70-degree angle) and a central marker at eye level directly in front of the subject. Theoretically, if turning the head during phonation alters the laryngeal anatomic relationship by bringing the vocal folds in closer proximity to one another, then airflow rate should lessen. Our results indicate that head position does not improve glottic closure in these patients, which is in contrast to previously published research.(1) Our results question the utility and underlying theoretical construct for the use of head turning as a therapeutic technique for improvement of voice in patients with unilateral vocal fold paralysis.     

Physiologic and acoustic differences between male and female voices

Comparison is drawn between male and female larynges on the basis of overall size, vocal fold membranous length, elastic properties of tissue, and prephonatory glottal shape. Two scale factors are proposed that are useful for explaining differences in fundamental frequency, sound power, mean airflow, and glottal efficiency. Fundamental frequency is scaled primarily according to the membranous length of the vocal folds (scale factor of 1.6), whereas mean airflow, sound power, glottal efficiency, and amplitude of vibration include another scale factor (1.2) that relates to overall larynx size. Some explanations are given for observed sex differences in glottographic waveforms. In particular, the simulated (computer-modeled) vocal fold contact area is used to infer male-female differences in the shape of the glottis. The female glottis appears to converge more linearly (from bottom to top) than the male glottis, primarily because of medial surface bulging of the male vocal folds.     

A theoretical study of the effects of various laryngeal configurations on the acoustics of phonation.

Simulation of glottal volume flow and vocal fold tissue movement was accomplished by numerical solution of a time-dependent boundary value problem, in which nonuniform, orthotropic, linear, incompressible vocal fold tissue media were surrounded by irregularly shaped boundaries, which were either fixed or subject to aerodynamic stresses. Spatial nonuniformity of the tissues was of the layered type, including a mucosal layer, a ligamental layer, and muscular layers. Orthotropy was required to stabilized the vocal folds longitudinally and to accomodate large variations in muscular stress. Incompressibility and vertical motions at the golttis played an important role in producing and sustaining phonation. A nominal configuration for male fundamental speaking pitches was selected, and the regulation of fundamental frequency, intensity, average volume flow, and vocal efficiency was investigated in terms of variations around this nominal configuration. Parameters which were varied consisted of geometrical factors such as length, thickness, and depth, factors for shaping the glottis, as well as tissue elasticities, tissue viscosities, and subglottal pressure. Since nonlinear stress-strain properties were not included, subglottal pressure did not produce a pronounced effect upon fundamental frequency under these somewhat edealized conditions F0 rasing correlated strongly with increased tension in the ligament, and somewhat with increasing tension in the vocalis. F0 lowering correlated with increase in vocal fold length when the tensions were held constant, but not with increase in vocal fold thickness. Vocal intensity and efficiency are shown to have local maxima as the configurational parameters are varied one at a time. It appears that oral acoustic power output and vocal efficiency can be maximized by proper adjustments of longitudinal tension of nonmuscular (mucosal and ligamental) tissue layers in relation to muscular layers. Quantitative verification of the "body-cover" theory is therefore suggested, and several further implications with regard to control of the human larynx are considered.