Phonation thresholds as a function of laryngeal size in a two-mass model of the vocal folds

This letter analyzes the oscillation onset-offset conditions of the vocal folds as a function of laryngeal size. A version of the two-mass model of the vocal folds is used, coupled to a two-tube approximation of the vocal tract in configuration for the vowel /a/. The standard male configurations of the laryngeal and vocal tract models are used as reference, and their dimensions are scaled using a single factor. Simulations of the vocal fold oscillation and oral output are produced for varying values of the scaling factor. The results show that the oscillation threshold conditions become more restricted for smaller laryngeal sizes, such as those appropriate for females and children.     

On the relation between the phonation threshold lung pressure and the oscillation frequency of the vocal folds

This Letter presents an extension of a previous equation for the phonation threshold pressure by Titze [I. R. Titze, J. Acoust. Soc. Am. 83, 1536-1552 (1988)]. The extended equation contains the vocal-fold oscillation frequency as an explicit factor. It is derived from the mucosal wave model of the vocal folds by considering the general case of an arbitrary time delay for the mucosal wave to travel the glottal height. The results are illustrated with a numerical example, which shows good qualitative agreement with experimental measures.     

Phonation threshold pressure: comparison of calculations and measurements taken with physical models of the vocal fold mucosa

In an important paper on the physics of small amplitude oscillations, Titze showed that the essence of the vertical phase difference, which allows energy to be transferred from the flowing air to the motion of the vocal folds, could be captured in a surface wave model, and he derived a formula for the phonation threshold pressure with an explicit dependence on the geometrical and biomechanical properties of the vocal folds. The formula inspired a series of experiments [e.g., R. Chan and I. Titze, J. Acoust. Soc. Am 119, 2351-2362 (2006)]. Although the experiments support many aspects of Titze's formula, including a linear dependence on the glottal half-width, the behavior of the experiments at the smallest values of this parameter is not consistent with the formula. It is shown that a key element for removing this discrepancy lies in a careful examination of the properties of the entrance loss coefficient. In particular, measurements of the entrance loss coefficient at small widths done with a physical model of the glottis (M5) show that this coefficient varies inversely with the glottal width. A numerical solution of the time-dependent equations of the surface wave model shows that adding a supraglottal vocal tract lowers the phonation threshold pressure by an amount approximately consistent with Chan and Titze's experiments.     

Mechanisms of irregular vibration in a physical model of the vocal folds

Previous investigations have shown that one mechanism of irregular vocal fold vibration may be a desynchronization of two or more vibratory modes of the vocal fold tissues. In the current investigation, mechanisms of irregular vibration were further examined using a self-oscillating, physical model of vocal fold vibration, a hemi-model methodology, and high-speed, stereoscopic, digital imaging. Using the method of empirical eigen-functions, a spatiotemporal analysis revealed mechanisms of irregular vibration in subharmonic phonation and biphonation, which were not disclosed in a standard acoustic spectrum.     

Coherent structures of the near field flow in a self-oscillating physical model of the vocal folds

Current theories of voice production depend critically upon knowledge of the near field flow which emanates from the glottis. While most modern theories predict complex, three-dimensional structures in the near field flow, few investigations have attempted to quantify such structures. Using methods of flow visualization and digital particle image velocimetry, this study measured the near field flow structures immediately downstream of a self-oscillating, physical model of the vocal folds, with a vocal tract attached. A spatio-temporal analysis of the structures was performed using the method of empirical orthogonal eigenfunctions. Some of the observed flow structures included vortex generation, vortex convection, and jet flapping. The utility of such data in the future development of more accurate, low-dimensional models of voice production is discussed.     

Phonation threshold pressure: a missing link in glottal aerodynamics.

Phonation threshold pressure has previously been defined as the minimum lung pressure required to initiate phonation. By modeling the dependence of this pressure on fundamental frequency, it is shown that relatively simple aerodynamic relations for time-varying flow in the glottis are obtained. Lung pressure and peak glottal flow are nearly linearly related, but not proportional. For this reason, traditional power law relations between vocal power and lung pressure may not hold. Glottal impendance for time-varying flow should be defined differentially rather than as a simple ratio between lung pressure and peak flow. It is shown that the peak flow, the peak flow derivative, the open quotient, and the speed quotient of inverse-filtered glottal flow waveforms all depend explicitly on phonation threshold pressure. Data from singers are compared with those from nonsingers. The primary difference is that singers obtain two to three times greater peak flow for a given lung pressure, suggesting that they adjust their glottal or vocal tract impedance for optimal flow transfer between the source and the resonantor.     

Unsteady behavior of flow in a scaled-up vocal folds model

Measurements of the fluid flow through a scaled-up model of the human glottis are presented to determine whether glottal flow may be approximated as unsteady. Time- and space-resolved velocity vector fields from digital particle image velocimetry (DPIV) measurements of the flow through the gap between two moving, rigid walls are presented in four cases, over a range of Strouhal numbers: 0.010, 0.018, 0.035, 0.040, corresponding to life-scale f(0) of 30, 58, 109, and 126 Hz, respectively, at a Reynolds number of 8000. It is observed that (1) glottal flow onset is delayed after glottal opening and (2) glottal flow shutoff occurs prior to closure. A comparison between flow through a fully open, nonmoving glottis and that through the moving vocal folds shows a marked difference in spatial structure of the glottal jet. The following features of the flow are seen to exhibit strong dependence on cycle frequency: (a) glottal exit plane velocity, (b) volume flow, (c) vortex shedding rates, and (d) vortex amplitude. Vortex shedding appears to be a factor both in controlling flow resistance and in cycle-to-cycle volume flow variations. All these observations strongly suggest that glottal flow is inherently unsteady.     

A theoretical model of the pressure field arising from asymmetric intraglottal flows applied to a two-mass model of the vocal folds

A theoretical flow solution is presented for predicting the pressure distribution along the vocal fold walls arising from asymmetric flow that forms during the closing phases of speech. The resultant wall jet was analyzed using boundary layer methods in a non-inertial reference frame attached to the moving wall. A solution for the near-wall velocity profiles on the flow wall was developed based on a Falkner-Skan similarity solution and it was demonstrated that the pressure distribution along the flow wall is imposed by the velocity in the inviscid core of the wall jet. The method was validated with experimental velocity data from 7.5 times life-size vocal fold models, acquired for varying flow rates and glottal divergence angles. The solution for the asymmetric pressures was incorporated into a widely used two-mass model of vocal fold oscillation with a coupled acoustical model of sound propagation. Asymmetric pressure loading was found to facilitate glottal closure, which yielded only slightly higher values of maximum flow declination rate and radiated sound, and a small decrease in the slope of the spectral tilt. While the impact on symmetrically tensioned vocal folds was small, results indicate the effect becomes more significant for asymmetrically tensioned vocal folds.     

Two-mass model of the vocal folds: negative differential resistance oscillation

The vocal folds and glottis are analyzed as a single system rather than as two separate but interacting systems, i.e., an aerodynamic one (the glottis) and a mechanical one (the vocal folds). Simplified steady flow calculations based on the two-mass model, and similar to those of Ishizaka and Matsudaira [SCRL Monograph No. 8, Santa Barbara, CA (1972)], are made except that flexible walls are assumed for both dc and ac flows. A negative differential resistance is found for steady flow when the coupling spring is weak compared to that of the lower mass. Dynamic transverse motion of the masses is represented by two transverse series resonant circuits in parallel within the glottis. The vocal tract is represented by a lumped resistance and inertance in series. Sustained, self-excited, small-amplitude oscillations can be obtained when the magnitude of the negative differential resistance is equal to the real part of the impedance of the rest of the circuit. The oscillation frequency depends only on the elasticity and mass of the vocal folds. The present analysis differs from Ishizaka and Matsudaira's analysis because their oscillation frequency decreases as dc volume velocity increases.     

Dependence of phonation threshold pressure on vocal tract acoustics and vocal fold tissue mechanics

Analytical and computer simulation studies have shown that the acoustic impedance of the vocal tract as well as the viscoelastic properties of vocal fold tissues are critical for determining the dynamics and the energy transfer mechanism of vocal fold oscillation. In the present study, a linear, small-amplitude oscillation theory was revised by taking into account the propagation of a mucosal wave and the inertive reactance (inertance) of the supraglottal vocal tract as the major energy transfer mechanisms for flow-induced self-oscillation of the vocal fold. Specifically, analytical results predicted that phonation threshold pressure (Pth) increases with the viscous shear properties of the vocal fold, but decreases with vocal tract inertance. This theory was empirically tested using a physical model of the larynx, where biological materials (fat, hyaluronic acid, and fibronectin) were implanted into the vocal fold cover to investigate the effect of vocal fold tissue viscoelasticity on Pth. A uniform-tube supraglottal vocal tract was also introduced to examine the effect of vocal tract inertance on Pth. Results showed that Pth decreased with the inertive impedance of the vocal tract and increased with the viscous shear modulus (G") or dynamic viscosity (eta') of the vocal fold cover, consistent with theoretical predictions. These findings supported the potential biomechanical benefits of hyaluronic acid as a surgical bioimplant for repairing voice disorders involving the superficial layer of the lamina propria, such as scarring, sulcus vocalis, atrophy, and Reinke's edema.     

Observation of perturbations in a lumped-element model of the vocal folds with application to some pathological cases

In this paper a mass-spring model is developed that is a hybrid of the two-mass and the longitudinal string models, proposed by Ishizaka and Flanagan [Bell Sys. Tech. J. 51, 1233-1268 (1972)] and Titze [Phonetica 28, 129-170 (1973)], respectively. The model is used to simulate the vibratory motion of both the normal and asymmetric vocal folds. Mouth-output pressure, lateral tissue displacement, phase plots, and energy diagrams are presented to demonstrate the interaction between vocal fold tissue and the aerodynamic flow between the folds. The results of the study suggest that this interaction is necessary for sustained large amplitude oscillation because the flow supplies the energy lost by the tissue damping. Tissue mass and stiffness were varied locally or uniformly. Decreased stress in the longitudinal string tension produced subharmonic and chaotic vibrations in the displacement, velocity and acceleration phase diagrams. Similar vibratory characteristics also appeared in pathological speech data analyzed using time domain jitter and shimmer measures and a harmonics-to-noise ratio metric. The subharmonics create an effect that has been perceptually described as diplophonia.     

Identification of geometric parameters influencing the flow-induced vibration of a two-layer self-oscillating computational vocal fold model

Simplified models have been used to simulate and study the flow-induced vibrations of the human vocal folds. While it is clear that the models' responses are sensitive to geometry, it is not clear how and to what extent specific geometric features influence model motion. In this study geometric features that played significant roles in governing the motion of a two-layer (body-cover), two-dimensional, finite element vocal fold model were identified. The model was defined using a flow solver based on the viscous, unsteady, Navier-Stokes equations and a solid solver that allowed for large strain and deformation. A screening-type design-of-experiments approach was used to identify the relative importance of 13 geometric parameters. Five output measures were analyzed to assess the magnitude of each geometric parameter's effect on the model's motion. The measures related to frequency, glottal width, flow rate, intraglottal angle, and intraglottal phase delay. The most significant geometric parameters were those associated with the cover--primarily the pre-phonatory intraglottal angle--as well as the body inferior angle. Some models exhibited evidence of improved model motion, including mucosal wave-like motion and alternating convergent-divergent glottal profiles, although further improvements are still needed to more closely mimic human vocal fold motion.     

Influence of acoustic loading on an effective single mass model of the vocal folds

Three-way interactions between sound waves in the subglottal and supraglottal tracts, the vibrations of the vocal folds, and laryngeal flow were investigated. Sound wave propagation was modeled using a wave reflection analog method. An effective single-degree-of-freedom model was designed to model vocal-fold vibrations. The effects of orifice geometry changes on the flow were considered by enforcing a time-varying discharge coefficient within a Bernoulli flow model. The resulting single-degree-of-freedom model allowed for energy transfer from flow to structural vibrations, an essential feature usually incorporated through the use of higher order models. The relative importance of acoustic loading and the time-varying flow resistance for fluid-structure energy transfer was established for various configurations. The results showed that acoustic loading contributed more significantly to the net energy transfer than the time-varying flow resistance, especially for less inertive supraglottal loads. The contribution of supraglottal loading was found to be more significant than that of subglottal loading. Subglottal loading was found to reduce the net energy transfer to the vocal-fold oscillation during phonation, balancing the effects of the supraglottal load.     

Influence of supraglottal structures on the glottal jet exiting a two-layer synthetic, self-oscillating vocal fold model

A synthetic two-layer, self-oscillating, life-size vocal fold model was used to study the influence of the vocal tract and false folds on the glottal jet. The model vibrated at frequencies, pressures, flow rates, and amplitudes consistent with human phonation, although some differences in behavior between the model and the human vocal folds are noted. High-speed images of model motion and flow visualization were acquired. Phase-locked ensemble-averaged glottal jet velocity measurements using particle image velocimetry (PIV) were acquired with and without an idealized vocal tract, with and without false folds. PIV data were obtained with varying degrees of lateral asymmetric model positioning. Glottal jet velocity magnitudes were consistent with those measured using excised larynges. A starting vortex was observed in all test cases. The false folds interfered with the starting vortex, and in some cases vortex shedding from the false folds was observed. In asymmetric cases without false folds, the glottal jet tended to skew toward the nearest wall; with the false folds, the opposite trend was observed. rms velocity calculations showed the jet shear layer and laminar core. The rms velocities were higher in the vocal tract cases compared to the open jet and false fold cases.     

Direct measurement of onset and offset phonation threshold pressure in normal subjects

Phonation threshold pressures were directly measured in five normal subjects in a variety of voicing conditions. The effects of fundamental frequency, intensity, closure speed of the vocal folds, and laryngeal airway resistance on phonation threshold pressures were determined. Subglottic air pressures were measured using percutaneous puncture of the cricothyroid membrane. Both onset and offset of phonation were studied to see if a hysteresis effect produced lower offset pressures than onset pressures. Univariate analysis showed that phonation threshold pressure was influenced most strongly by fundamental frequency and intensity. Multiple linear regression showed that these two variables, as well as laryngeal airway resistance, most strongly predicted phonation threshold pressure. Two of the five subjects demonstrated a significant hysteresis effect, but one subject actually had higher offset pressures than onset pressures.     

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.     

The physics of small-amplitude oscillation of the vocal folds

A theory of vocal fold oscillation is developed on the basis of the body-cover hypothesis. The cover is represented by a distributed surface layer that can propagate a mucosal surface wave. Linearization of the surface-wave displacement and velocity, and further small-amplitude approximations, yields closed-form expressions for conditions of oscillation. The theory predicts that the lung pressure required to sustain oscillation, i.e., the oscillation threshold pressure, is reduced by reducing the mucosal wave velocity, by bringing the vocal folds closer together and by reducing the convergence angle in the glottis. The effect of vocal tract acoustic loading is included. It is shown that vocal tract inertance reduces the oscillation threshold pressure, whereas vocal tract resistance increases it. The treatment, which is applicable to falsetto and breathy voice, as well as onset or release of phonation in the absence of vocal fold collision, is harmonized with former treatments based on two-mass models and collapsible tubes.     

Verification of two minimally invasive methods for the estimation of the contact pressure in human vocal folds during phonation

The contact pressure on the vocal fold surface during high pitch or amplitude voice production is believed to be one major source of phonotrauma. Models for the quantitative estimate of the contact pressure may be valuable for prevention and treatment. Various indirect and minimally invasive approaches have been purported to estimate contact pressure. But the accuracy of these methods has not yet been objectively verified in controlled laboratory settings. In the present study, two indirect approaches for the estimation of the contact pressure were investigated. One is based on a Hertzian impact model, and the other on a finite element model. A probe microphone was used for direct measurements of the contact pressure and verifications of the indirect approaches. A silicone replica of human vocal folds was used as a test bed. Consistent contact pressure estimations were obtained using all three methods. The advantages and disadvantages of each approach for eventual clinical applications are described.