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.     

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.     

Nonlinear oscillation of pathological vocal folds during vocalization

The extended two-mass model is adopted to analyze the nonlinear oscillation of pathological vocal folds during vocalization. Redundant tissue or area in laryngeal patients is modeled as a massless rigid connected to the upper mass of the vocal folds, and a parameter Q is introduced to represent the change of glottal configurations and tension imbalance between the left and right sides of vocal folds. Numerical simulations demonstrate that the pathological vocal-fold decreases the threshold of Q to generate nonlinear vocal oscillation, indicating the improvement of the sensitivity of vocal folds to asymmetries and enhancing the coupling between two sides. Furthermore, the pathological vocal-fold can lower the fundamental frequency and eliminate high-order harmonics, For example, the fundamental frequency decreases from 119.94 Hz to 84.95 Hz when Q=0.58 and the sub-glottal pressure 1450 Pa. However, there are no prominent effects on the amplitudes of sub-harmonic and low-order harmonics.     

The effect of viscosity changes in the vocal folds on the range of oscillation

Changes in vocal fold oscillation threshold pressure were induced in excised canine larynges by experimentally causing fluid movement into and out of the vocal folds. The transport was facilitated by exposing the vocal folds to various osmotic solutions, and it was assumed that changes in hydration caused changes in the internal tissue viscosity. A range of oscillation threshold pressures was measured for each condition of hydration by varying length and glottal width. The oscillation threshold pressure shifted as predicted. Decreased hydration (increased viscosity) raised the threshold of oscillation, and increased hydration (decreased viscosity) lowered the threshold of oscillation. This apparently represents the first in vitro model for the study of the effect of viscosity changes of the internal environment of the vocal folds on phonation.     

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.     

Botox for hyperadduction of the false vocal folds: A case report

We present a patient with severe hyperadduction of the false vocal folds (FVF) treated with Botulinum Toxin injections to each FVF. This patient presented with severe dysphonia and was found to demonstrate severe hyperadduction of the FVF's with all phonatory tasks. The patient was treated with extensive speech therapy without improvement in voice quality nor FVF motion pattern. He was then injected with Botox A bilaterally using a peroral approach to the FVFs. Shortly after treatment the patient experienced dramatic improvement in voice quality. Videolaryngoscopy revealed no adduction of the FVFs with phonation and essentially normal true vocal fold motion. He remained with normal voice quality one year after treatment without any further treatment. Possible mechanism of action of this type of treatment are discussed.     

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.     

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.     

Nonlinear behavior of vocal fold vibration: The role of coupling between the vocal folds

Coupling between the vocal folds is one of the nonlinear mechanisms allowing regulation and synchronization of mucosal vibration. The purpose of this study was to establish that modulations such as diplophonia and abnormalities observed in vocal signals that may be observed in some cases of laryngeal pathology can be considered as nonlinear behavior due to the persistence of some physical interaction (coupling). An experimental model using excised porcine larynx was designed to create tension asymmetry between the vocal folds and to obtain vocal signals with modulations. Signals were analyzed by spectral analysis and the phase portrait method. Results were compared with computer-generated synthetic signals corresponding to nonlinear combinations of sinusoid signals. Under these conditions, evidence of nonlinear behavior was detected in 85% of experimental signals. These findings were interpreted as a demonstration of vocal fold interaction. Based on these findings, the authors conclude that (1) coupling must be taken into account in physical models of laryngeal physiology, and that (2) methods of nonlinear dynamics may be used for objective voice analysis.     

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.     

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.     

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.     

A phenomenology of boundary lubrication: the lumped junction model

We propose a phenomenological model of boundary lubricated junctions consisting of a few layers of small molecules which describes the rheological properties of these sytems both in the static, frozen, and sliding, molten, states as well as the dynamical transition between them. Two dynamical regimes can be distinguished, according to the level of internal damping of the junction, which depends on its thickness and on the normal load. In the overdamped regime, under driving at constant velocity v through an external spring, the motion evolves continuously from “atomic stick-slip” to modulated sliding. Underdamped systems exhibit, under given external stress, a range of dynamic bistability where the sheared static state coexists with a steadily sliding one. The frictional dynamics under shear driving is analyzed in detail, it provides a complete account of the qualitative dynamical scenarios observed by Israelashvili et al., and yields semiquantitative agreement with experimental data. A few complementary experimental tests of the model are suggested. Received: 18 December 1997 / Received in final form and accepted: 26 March 1998     

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.     

A Preisach-based hysteresis model for magnetic and ferroelectric hysteresis

In this paper we present a model for hysteretic nonlinearities with non-local memories. This model can be applied to describe hysteretic material behavior. Common applications are ferromagnetic or ferroelectric materials. Our model consists of an analytic function and a Preisach operator. We define a continuous Preisach weight function and introduce a method for the identification of the model parameters. The model parameters are customized to a set of symmetric hysteresis curves. We verify our model for a soft magnetic material, a hard magnetic material and the ferroelectric behavior of some piezoelectric material. After that, non-symmetric curves like the virgin curve are predicted very well by the model. It is especially useful, if forced magnetization or polarization, that appears beyond technical saturation, come into account.     

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.     

Stroboscopic signs associated with benign lesions of the vocal folds

Stroboscopic signs were systematically rated for a group of 80 patients with benign vocal fold lesions, most of whom had either a nodule or a polyp. Each group revealed a characteristic pattern of ranking of signs and exhibited differences of most predominant signs. The results of the ratings were submitted to a multiple discriminant analysis to determine if post hoc stroboscopic ratings could be used to correctly classify patients into one of four diagnostic groups and into one of two treatment groups. All patients except one were correctly classified into the diagnostic groups, and all were correctly classified into the treatment groups. The important signs for classifying patients into the diagnostic groups were roughness of the edge of the affected vocal fold, phase closure pattern, and phase symmetry. The important signs for classifying patients into the treatment groups were roughness of the edge of the affected vocal fold, glottal closure configuration, and vibration characteristics of the affected (or more affected) vocal fold. The results suggest that objective evaluation of stroboscopic examinations can be valuable in correctly diagnosing patients and in selecting the proper treatment regimen for the patient.     

声带上表面振动黏膜波波速的激光测振估计方法

发声过程中声带组织振动的黏膜波测量对于声带组织力学特性和病理机制研究具有重要意义。本研究利用多普勒激光测振(Laser Doppler Vibrometer,LDV)和电声门图时间同步方法,对声带上表面的振动过程进行了重建,并基于该结果估计了黏膜波传播速度。通过离体犬喉发声实验,激光测振方法能够得到与高速光学方法相一致的黏膜波速度估计,并且能够重建不同时刻的声带上表面形态,证明了该方法应用于声带上表面振动特性研究的可行性和有效性。然而,由于单点测振的局限和时间同步的要求,稳态发声是保证该方法测量准确性的重要条件。