Determination of superior surface strains and stresses, and vocal fold contact pressure in a synthetic larynx model using digital image correlation

Stresses and strains within the vocal fold tissue may play a critical role in voice fatigue, in tissue damage and resulting voice disorders, and in tissue healing. In this study, experiments were performed to determine mechanical fields on the superior surface of a self-oscillating physical model of the human vocal folds using a three-dimensional digital image correlation method. Digital images obtained using a high-speed camera together with a mirror system were used to measure displacement fields, from which strains, strain rates, and stresses on the superior surface of the model vocal folds were computed. The dependence of these variables on flow rate was established. A Hertzian impact model was used to estimate the contact pressure on the medial surface from superior surface strains. A tensile stress dominated state was observed on the superior surface, including during collision between the model folds. Collision between the model vocal folds limits the medial-lateral stress levels on the superior surface, in conjunction with compressive stress or contact pressure on the medial surface.     

A model for vocal fold vibratory motion, contact area, and the electroglottogram

The electroglottogram (EGG) has been conjectured to be related to the area of contact between the vocal folds. This hypothesis has been substantiated only partially via direct and indirect observations. In this paper, a simple model of vocal fold vibratory motion is used to estimate the vocal fold contact area as a function of time. This model employs a limited number of vocal fold vibratory features extracted from ultra high-speed laryngeal films. These characteristics include the opening and closing vocal fold angles and the lag (phase difference) between the upper and lower vocal fold margins. The electroglottogram is simulated using the contact area, and the EGG waveforms are compared to measured EGGs for normal male voices producing both modal and pulse register tones. The model also predicts EGG waveforms for vocal fold vibration associated with a nodule or polyp.     

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 contribution to simulating a three-dimensional larynx model using the finite element method

A three-dimensional model is presented to simulate the larynx during vocalization. The finite element method is used to calculate the airflow velocity and pressure along the larynx as well as tissue displacement. It is assumed that the larynx tissue is transversely isotropic and divided into three tissues: cover, ligament, and body. A contact-impact algorithm is incorporated to deal with the physics of the collision between both true vocal folds. The results show that the simulated larynx can reproduce the vertical and horizontal phase difference in the tissue movements and that the false vocal folds affect the pressure distribution over the larynx surfaces. The effects of exciting the larynx with different pressure drops are also investigated.     

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.     

On the difference between negative damping and eigenmode synchronization as two phonation onset mechanisms

Negative damping and eigenmode synchronization as two different mechanisms of phonation onset are distinguished. Although both mechanisms lead to a favorable phase relationship between the flow pressure and the vocal fold motion as required for a net energy transfer into the vocal folds, the underlying mechanisms for this favorable phase relationship are different. The negative damping mechanism relies on glottal aerodynamics or acoustics to establish before onset and maintain after onset the favorable phase relationship, and therefore has minimum requirements on vocal fold geometry and biomechanics. A single degree-of-freedom vocal fold model is all that is needed for self-oscillation in the presence of a negative damping mechanism. In contrast, the mechanism of eigenmode synchronization critically depends on the geometrical and biomechanical properties of the vocal folds (at least 2-degrees-of-freedom are required), and has little requirement on the glottal aerodynamics other than flow separation. The favorable phase relation is established once synchronization occurs, regardless of the phase relationship imposed by glottal aerodynamics before onset. Unlike that of the negative damping mechanism, initiation of eigenmode synchronization requires neither a velocity-dependent flow pressure nor an alternating convergent-divergent glottis. The clinical implications of the distinctions between these two mechanisms are discussed.     

Sensitivity of a continuum vocal fold model to geometric parameters, constraints, and boundary conditions

The influence of key dimensional parameters, motion constraints, and boundary conditions on the modal properties of an idealized, continuum model of the vocal folds was investigated. The Ritz method and the finite element method were used for the analysis. The model's vibratory modes were determined to be most sensitive to changes in the anterior-posterior length of the vocal fold model, due to the influence of three-dimensional stress components acting in the transverse plane. Anterior/ posterior boundary conditions were found to have a significant influence on the vibratory response. Overestimation of modal frequencies resulted when vibration of the structure was restricted to the transverse plane. The overestimation of each modal frequency was proportional to the ratio of longitudinal to transverse Young's modulus, and was significant for ratio values less than 20.     

Flow separation in a computational oscillating vocal fold model

A finite-volume computational model that solves the time-dependent glottal airflow within a forced-oscillation model of the glottis was employed to study glottal flow separation. Tracheal input velocity was independently controlled with a sinusoidally varying parabolic velocity profile. Control parameters included flow rate (Reynolds number), oscillation frequency and amplitude of the vocal folds, and the phase difference between the superior and inferior glottal margins. Results for static divergent glottal shapes suggest that velocity increase caused glottal separation to move downstream, but reduction in velocity increase and velocity decrease moved the separation upstream. At the fixed frequency, an increase of amplitude of the glottal walls moved the separation further downstream during glottal closing. Increase of Reynolds number caused the flow separation to move upstream in the glottis. The flow separation cross-sectional ratio ranged from approximately 1.1 to 1.9 (average of 1.47) for the divergent shapes. Results suggest that there may be a strong interaction of rate of change of airflow, inertia, and wall movement. Flow separation appeared to be "delayed" during the vibratory cycle, leading to movement of the separation point upstream of the glottal end only after a significant divergent angle was reached, and to persist upstream into the convergent phase of the cycle.     

Voice source characteristics in Mongolian "throat singing" studied with high-speed imaging technique, acoustic spectra, and inverse filtering.

Mongolian "throat singing" can be performed in different modes. In Mongolia, the bass-type is called Kargyraa. The voice source in bass-type throat singing was studied in one male singer. The subject alternated between modal voice and the throat singing mode. Vocal fold vibrations were observed with high-speed photography, using a computerized recording system. The spectral characteristics of the sound signal were analyzed. Kymographic image data were compared to the sound signal and flow inverse filtering data from the same singer were obtained on a separate occasion. It was found that the vocal folds vibrated at the same frequency throughout both modes of singing. During throat singing the ventricular folds vibrated with complete but short closures at half the frequency of the true vocal folds, covering every second vocal fold closure. Kymographic data confirmed the findings. The spectrum contained added subharmonics compared to modal voice. In the inverse filtered signal the amplitude of every second airflow pulse was considerably lowered. The ventricular folds appeared to modulate the sound by reducing the glottal flow of every other vocal fold vibratory cycle.     

Voice Source Characteristics in Mongolian “Throat Singing” Studied with High-Speed Imaging Technique, Acoustic Spectra, and Inverse Filtering

Mongolian “throat singing” can be performed in different modes. In Mongolia, the bass-type is called Kargyraa. The voice source in bass-type throat singing was studied in one male singer. The subject alternated between modal voice and the throat singing mode. Vocal fold vibrations were observed with high-speed photography, using a computerized recording system. The spectral characteristics of the sound signal were analyzed. Kymographic image data were compared to the sound signal and flow inverse filtering data from the same singer were obtained on a separate occasion. It was found that the vocal folds vibrated at the same frequency throughout both modes of singing. During throat singing the ventricular folds vibrated with complete but short closures at half the frequency of the true vocal folds, covering every second vocal fold closure. Kymographic data confirmed the findings. The spectrum contained added subharmonics compared to modal voice. In the inverse filtered signal the amplitude of every second airflow pulse was considerably lowered. The ventricular folds appeared to modulate the sound by reducing the glottal flow of every other vocal fold vibratory cycle.     

Effect of inferior surface angle on the self-oscillation of a computational vocal fold model

Geometry of the human vocal folds strongly influences their oscillatory motion. While the effect of intraglottal geometry on phonation has been widely investigated, the study of the geometry of the inferior surface of the vocal folds has been limited. In this study the way in which the inferior vocal fold surface angle affects vocal fold vibration was explored using a two-dimensional, self-oscillating finite element vocal fold model. The geometry was parameterized to create models with five different inferior surface angles. Four of the five models exhibited self-sustained oscillations. Comparisons of model motion showed increased vertical displacement and decreased glottal width amplitude with decreasing inferior surface angle. In addition, glottal width and air flow rate waveforms changed as the inferior surface angle was varied. Structural, rather than aerodynamic, effects are shown to be the cause of the changes in model response as the inferior surface angle was varied. Supporting data including glottal pressure distribution, average intraglottal pressure, energy transfer, and flow separation point locations are discussed, and suggestions for future research are given.     

The effect of subglottal resonance upon vocal fold vibration

An experiment with excised larynges was undertaken to investigate theinteraction between acoustic pressures in a pseudotrachea and the amplitude of vibration of the vocal folds. Pressure was measured beneath the vocal folds at three specific moments of the vibratory cycle: (a) when the superior margin of the vocal folds began to separate, (b) when the vocal folds were maximally apart, and (c) when the inferior margin of the vocal folds began to touch. Results indicate that in half the larynges investigated, the maximum amplitude of vibration increased as a function of: (a) increased positive pressure at the moment of opening, and (b) reduced subglottal pressure when the vocal folds were maximally apart. The implications of these experiments regarding involuntary register transitions related to trachea resonance are discussed in light of a previously proposed register theory.     

Treatment of sulcus vocalis: Auditory perceptual and acoustical analysis of the slicing mucosa surgical technique

Sulcus vocalis is a furrow in the mucosa of the vocal folds that causes glottic insufficiency, irregular vibrations, and a distinctive dysphonia of variable severity. Treatment of sulcus vocalis needs to achieve anatomical and functional improvements that satisfy the behavior of the larynx and vocal quality. We suggest a combined approach, which includes vocal fold mucosal surgery through a technique called slicing mucosa, followed by intensive vocal rehabilitation. The aim of the surgery is to detach the mucosa of the sulcus and to interrupt the longitudinal fibrotic tension lines. The results of this procedure in 10 patients were analyzed through visual, auditory, and acoustical methods, and the data lead us to conclude that the treatment proposed offered a good anatomical result, a better vibratory pattern, and an improvement in overall vocal quality.     

Chaotic vibrations of a vocal fold model with a unilateral polyp

A nonlinear model was proposed to study chaotic vibrations of vocal folds with a unilateral vocal polyp. The model study found that the vocal polyp affected glottal closure and caused aperiodic vocal fold vibrations. Using nonlinear dynamic methods, aperiodic vibrations of the vocal fold model with a polyp were attributed to low-dimensional chaos. Bifurcation diagrams showed that vocal polyp size, stiffness, and damping had important effects on vocal fold vibrations. An increase in polyp size tended to induce subharmonic patterns and chaos. This study provides a theoretical basis to model aperiodic vibrations of vocal folds with a laryngeal mass.     

In Vivo Quantification of the Intraglottal Pressure: Modal Phonation and Voice Onset

Intraglottal pressure is the driving force of vocal fold vibration. Its time course during the open phase of the vibratory cycle is essential in the mechanics of phonation, but measuring it directly is difficult and may hinder spontaneous voicing. However, it can be computed from the in vivo measured transglottal flow and glottal area (hence the air particle velocity) on the basis of the Bernoulli energy law and the interaction with the inertance of the vocal tract. As to sustained modal phonation, calculations are presented for the two possible shapes of glottal duct: convergent and divergent, including absolute calibration in order to obtain quantitative physical values. Whatever the glottal duct configuration, the calculations based on measured values of glottal area and air flow show that the integrated intraglottal pressure during the opening phase systematically exceeds that during the closing phase, which is the basic condition for sustaining vocal fold oscillation. The key point is that the airflow curve is skewed to the right relative to the glottal area curve. The skewing results from air compressibility and vocal tract inertance. The intraglottal pressure becomes negative during the closing phase. As to the soft (or physiological) voice onset, a similar approach shows that the integrated pressure differences (opening phase − closing phase) actually increase as the onset progresses, and this applies to the results based on Bernoulli's energy law as well as to those based on the interaction with the inertance of the vocal tract. Furthermore and similarly, the phase lead of the pressure wave with respect to the glottal opening progressively increases. The underlying explanation lies in the progressively increasing skewing of the airflow curve to the right with respect to the glottal area curve.     

On the application of the lattice Boltzmann method to the investigation of glottal flow

The production of voice is directly related to the vibration of the vocal folds, which is generated by the interaction between the glottal flow and the tissue of the vocal folds. In the current study, the aerodynamics of the symmetric glottis is investigated numerically for a number of static configurations. The numerical investigation is based on the lattice Boltzmann method (LBM), which is an alternative approach within computational fluid dynamics. Compared to the traditional Navier-Stokes computational fluid dynamics methods, the LBM is relatively easy to implement and can deal with complex geometries without requiring a dedicated grid generator. The multiple relaxation time model was used to improve the numerical stability. The results obtained with LBM were compared to the results provided by a traditional Navier-Stokes solver and experimental data. It was shown that LBM results are satisfactory for all the investigated cases.     

Vocal fold and ventricular fold vibration in period-doubling phonation: physiological description and aerodynamic modeling

Occurrences of period-doubling are found in human phonation, in particular for pathological and some singing phonations such as Sardinian A Tenore Bassu vocal performance. The combined vibration of the vocal folds and the ventricular folds has been observed during the production of such low pitch bass-type sound. The present study aims to characterize the physiological correlates of this acoustical production and to provide a better understanding of the physical interaction between ventricular fold vibration and vocal fold self-sustained oscillation. The vibratory properties of the vocal folds and the ventricular folds during phonation produced by a professional singer are analyzed by means of acoustical and electroglottographic signals and by synchronized glottal images obtained by high-speed cinematography. The periodic variation in glottal cycle duration and the effect of ventricular fold closing on glottal closing time are demonstrated. Using the detected glottal and ventricular areas, the aerodynamic behavior of the laryngeal system is simulated using a simplified physical modeling previously validated in vitro using a larynx replica. An estimate of the ventricular aperture extracted from the in vivo data allows a theoretical prediction of the glottal aperture. The in vivo measurements of the glottal aperture are then compared to the simulated estimations.     

Voice production model integrating boundary-layer analysis of glottal flow and source-filter coupling

A voice production model is created in this work by considering essential aerodynamic and acoustic phenomena in human voice production. A precise flow analysis is performed based on a boundary-layer approximation and the viscous-inviscid interaction between the boundary layer and the core flow. This flow analysis can supply information on the separation point of the glottal flow and the thickness of the boundary layer, both of which strongly depend on the glottal configuration and yield an effective prediction of the flow behavior. When the flow analysis is combined with the modified two-mass model of the vocal fold [Pelorson et al. (1994). J. Acoust. Soc. Am. 96, 3416-3431], the resulting acoustic wave travels through the vocal tract and a pressure change develops in the vicinity of the glottis. This change can affect the glottal flow and the motion of the vocal folds, causing source-filter coupling. The property of the acoustic feedback is explicitly expressed in the frequency domain by using an acoustic tube model, allowing a clear interpretation of the coupling. Numerical experiments show that the vocal-tract input impedance and frequency responses representing the source-filter coupling have dominant peaks corresponding to the fourth and fifth formants. Results of time-domain simulations also suggest the importance of these high-frequency peaks in voice production.