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.     

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.     

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.     

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 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.     

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.     

Influence of plastic deformation on fracture of an aluminum single crystal under shock-wave loading

The plastic deformation and the onset of fracture of single-crystal metals under shock-wave loading have been studied using aluminum as an example by the molecular dynamics method. The mechanisms of plastic deformation under compression in a shock wave and under tension in rarefaction waves have been investigated. The influence of the defect structure formed in the compression wave on the spall strength and the fracture mechanism has been analyzed. The dependence of the spall strength on the strain rate has been obtained.     

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.     

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.     

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.     

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.     

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.     

The dynamics of length change in canine vocal folds

The time courses of vocal fold elongation and contraction have beenmeasured as a function of intrinsic laryngeal muscle activity. The superior and recurrent laryngeal nerves of anesthetized canines were stimulated supramaximally (on-off in all combinations) while the vocal folds were surgically exposed and illuminated for conventional and higher speed (300 frames per second) video recording. Microsutures were placed on various points on the vocal folds to measure elongation and contraction. Vocal fold strain, defined as elongation divided by rest length, ranged from −17% to +45%. The typical time constant for exponential increase or decrease in strain was about 30 ms. This reflects primarily the intrinsic muscle activation times rather than a passive (inertial or viscoelastic) response of cricothyroid joint rotation or translation.     

Mathematical model of acoustic speech production with mobile walls of the vocal tract

A mathematical speech production model is considered that describes acoustic oscillation propagation in a vocal tract with mobile walls. The wave field function satisfies the Helmholtz equation with boundary conditions of the third kind (impedance type). The impedance mode corresponds to a threeparameter pendulum oscillation model. The experimental research demonstrates the nonlinear character of how the mobility of the vocal tract walls influence the spectral envelope of a speech signal.     

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.     

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.