The influence of subglottal acoustics on laboratory models of phonation

Many previous laboratory investigations of phonation involving physical models, excised larynges, and in vivo canine larynges have failed to fully specify the subglottal system. Many of these same studies have reported a variety of nonlinear phenomena, including bifurcations (e.g., various classes of phonation onset and offset, register changes, frequency jumps), subharmonics, and chaos, and attributed such phenomena to the biomechanical properties of the larynx. However, such nonlinear phenomena may also be indicative of strong coupling between the voice source and the subglottal tract. Consequently, in such studies, it has not been clear whether the underlying mechanisms of such nonlinear phenomena were acoustical, biomechanical, or a coupling of the acoustical and biomechanical systems. Using a physical model of vocal fold vibration, and tracheal tube lengths which have been commonly reported in the literature, it is hypothesized and subsequently shown that such nonlinear phenomena may be replicated solely on the basis of laryngeal interactions with the acoustical resonances of the subglottal system. Recommendations are given for ruling out acoustical resonances as the source of nonlinear phenomena in future laboratory studies of phonation.     

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.     

Using the Relaxation Oscillations Principle for Simple Phonation Modeling

Well-known multimass models of vocal folds are useful to describe main behavior observed in human voicing but their principle of functioning, based on harmonic oscillation, may appear complex. This work is designed to show that a simple one-mass model ruled by laws of relaxation oscillation can also depict main behavior of glottis dynamic. Theory of relaxation oscillation is detailed. A relaxation oscillation model is assessed through a numerical simulation using conventional values for tissue characteristics and subglottal pressure. As expected, raising the mass decreases the fundamental frequency and increases the amplitude of vocal fold vibration: for a mass ranging from 0.01 to 0.4 g, F0 decreased from 297.5 to 42.5 Hz and vibrational amplitude increased from 1.26 to 3.25 mm (for stiffness k=10Nm(-1), damping r=0.015 N s m(-1), and subglottal pressure=1 kPa). Stiffness value has the opposite effect. The subglottal pressure controls the fundamental frequency with a rate ranging from 20 to 50 Hz/kPa. The vibrational amplitude is also controlled linearly by subglottal pressure from 0.22 to 0.26 mm/kPa. The range of phonation threshold pressure (PTP) is close to the values currently proposed, that is, 0.1 to 1 kPa and varies with the fundamental frequency. The relaxation oscillator is a simple and useful tool for modeling vocal fold vibration.     

Glottal flow through a two-mass model: comparison of Navier-Stokes solutions with simplified models

A new numerical model of the vocal folds is presented based on the well-known two-mass models of the vocal folds. The two-mass model is coupled to a model of glottal airflow based on the incompressible Navier-Stokes equations. Glottal waves are produced using different initial glottal gaps and different subglottal pressures. Fundamental frequency, glottal peak flow, and closed phase of the glottal waves have been compared with values known from the literature. The phonation threshold pressure was determined for different initial glottal gaps. The phonation threshold pressure obtained using the flow model with Navier-Stokes equations corresponds better to values determined in normal phonation than the phonation threshold pressure obtained using the flow model based on the Bernoulli equation. Using the Navier-Stokes equations, an increase of the subglottal pressure causes the fundamental frequency and the glottal peak flow to increase, whereas the fundamental frequency in the Bernoulli-based model does not change with increasing pressure.     

An anatomically based, time-domain acoustic model of the subglottal system for speech production

A time-domain model of sound wave propagation in the branching airways of the subglottal system is presented. The model is formulated as an extension to an acoustic transmission-line modeling scheme originally developed for simulating the supraglottal system in the time-domain during speech production [Maeda (1982). Speech Commun. 1, 199-229; Mokhtari et al. (2008). Speech Commun. 50, 179-190]. The approach allows for predictions of time-varying acoustic pressure and volume velocity at any point along the various generations of subglottal airways from trachea to alveoli. In addition, the model can be configured so that its overall structure simulates different geometric forms, including airways that branch in a symmetric or asymmetric pattern. Three subglottal configurations, two symmetric and one asymmetric, were represented based on reported anatomical dimensions of the subglottal airways. Estimates of the acoustic input impedances of these subglottal configurations revealed resonant characteristics similar to those found in the previous studies. Simulations of voiced sound propagation into the subglottal airways, achieved by coupling the subglottal model to a two-mass vocal fold model and a supraglottal tract configured for different vowels, yielded predictions of time-domain sound pressure waveforms below the vocal folds that compare favorably to previous measurements in human subjects.     

Nonlinear source-filter coupling in phonation: theory

A theory of interaction between the source of sound in phonation and the vocal tract filter is developed. The degree of interaction is controlled by the cross-sectional area of the laryngeal vestibule (epilarynx tube), which raises the inertive reactance of the supraglottal vocal tract. Both subglottal and supraglottal reactances can enhance the driving pressures of the vocal folds and the glottal flow, thereby increasing the energy level at the source. The theory predicts that instabilities in vibration modes may occur when harmonics pass through formants during pitch or vowel changes. Unlike in most musical instruments (e.g., woodwinds and brasses), a stable harmonic source spectrum is not obtained by tuning harmonics to vocal tract resonances, but rather by placing harmonics into favorable reactance regions. This allows for positive reinforcement of the harmonics by supraglottal inertive reactance (and to a lesser degree by subglottal compliant reactance) without the risk of instability. The traditional linear source-filter theory is encumbered with possible inconsistencies in the glottal flow spectrum, which is shown to be influenced by interaction. In addition, the linear theory does not predict bifurcations in the dynamical behavior of vocal fold vibration due to acoustic loading by the vocal tract.     

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.     

Measures of spectral slope using an excised larynx model

Spectral measures of the glottal source were investigated using an excised canine larynx (CL) model for various aerodynamic and phonatory conditions. These measures included spectral harmonic difference H1-H2 and spectral slope that are highly correlated with voice quality but not reported in a systematic manner using an excised larynx model. It was hypothesized that the acoustic spectra of the glottal source were significantly influenced by the subglottal pressure, glottal adduction, and vocal fold elongation, as well as the resulting vibration pattern. CLs were prepared, mounted on the bench with and without false vocal folds, and made to oscillate with a flow of heated and humidified air. Major control parameters were subglottal pressure, adduction, and elongation. Electroglottograph, subglottal pressure, flow rate, and audio signals were analyzed using custom software. Results suggest that an increase in subglottal pressure and glottal adduction may change the energy balance between harmonics by increasing the spectral energy of the first few harmonics in an unpredictable manner. It is suggested that changes in the dynamics of vocal fold motion may be responsible for different spectral patterns. The finding that the spectral harmonics do not conform to previous findings was demonstrated through various cases. Results of this study may shed light on phonatory spectral control when the larynx is part of a complete vocal tract system.     

Correlation analysis of the physiological factors controlling fundamental voice frequency

A technique has been developed to obtain a quantitative measure of correlation between electromyographic (EMG) activity of various laryngeal muscles, subglottal air pressure, and the fundamental frequency of vibration of the vocal folds (Fo). Data were collected and analyzed on one subject, a native speaker of American English. The results show that an analysis of this type can provide a useful measure of correlation between the physiological and acoustical events in speech and, furthermore, can yield detailed insights into the organization and nature of the speech production process. In particular, based on these results, a model is suggested of Fo control involving laryngeal state functions that seems to agree with present knowledge of laryngeal control and experimental evidence.     

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.     

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.     

马铁菊头蝠声道声学特性的数值分析

采用有限元数值计算得到了马铁菊头蝠声道内部的声场分布,给出了马铁菊头蝠声道内几种特殊的腔体结构在蝙蝠发声过程中的作用。通过微型CT扫描并经过三维重构得到了马铁菊头蝠声道的三维立体模型用于有限元数值计算,通过在声门处放置单位声源计算得到了整个声道内部以及鼻孔周围的声压分布。结果表明,马铁菊头蝠声道包含了鼻腔结构后声波在声门上方的声压幅度明显大于不含鼻腔结构的情况,从传输曲线来看,声门上方鼻腔的存在使得系统对声波传输在二次谐波频率处呈现低阻抗效果,同时鼻腔的改变还可影响二次谐波的位置。而声门下方的气管空腔主要影响声波的背向转播,声门下方的气管空腔的存在可明显降低蝙蝠发声时声场在声道声门下方的声压幅度,同时抑制声音背向传播时二次谐波成分的强度。      

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.     

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.     

Dynamic glottal pressures in an excised hemilarynx model

During phonation, air pressures act upon the vocal folds to help maintain their oscillation. The air pressures vary dynamically along the medial surface of the vocal folds, although no live human or excised studies have shown how those pressure profiles vary in time. The purpose of this study was to examine time-dependent glottal pressure profiles using a canine hemilarynx approach. The larynx tissue was cut in the midsaggital plane from the top to about 5 mm below the vocal folds. The right half was replaced with a Plexiglas pane with imbedded pressure taps. Simultaneous recordings were made of glottal pressure signals, subglottal pressure, particle velocity, and average airflow at various levels of adduction. The data indicate that the pressures in the glottis (on the Plexiglas) vary both vertically and longitudinally throughout the phonatory cycle. Pressures vary most widely near the location of maximum vibratory amplitude, and can include negative pressures during a portion of the cycle. Pressures anterior and posterior to the maximum amplitude location may have less variation and may remain positive throughout the cycle, giving rise to a new concept called dynamic bidirectional pressure gradients in the glottis. This is an important concept that may relate strongly to tissue health as well as basic oscillatory mechanics.     

A methodological and preliminary study on the acoustic effect of a trumpet player's vocal tract

A methodological study is presented to examine the acoustic role of the vocal tract in playing the trumpet. Preliminary results obtained for one professional player are also shown to demonstrate the effectiveness of the method. Images of the vocal tract with a resolution of 0.5 mm (2 mm in thickness) were recorded with magnetic resonance imaging to observe the tongue posture and estimate the vocal-tract area function during actual performance. The input impedance was then calculated for the player's air column including both the supra- and subglottal tracts using an acoustic tube model including the effect of wall losses. Finally, a time-domain blowing simulation by Adachi and Sato [J. Acoust. Soc. Am. 99, 1200-1209 (1996)] was performed with a model of the lips. In this simulation, the oscillating frequency of the lips was slightly affected by using different shapes of the vocal tract measured for the player. In particular, when the natural frequency of the lips was gradually increased, the transition to the higher mode occurred at different frequencies for different vocal-tract shapes. Furthermore, simulation results showed that the minimum blowing pressure required to attain the lip oscillation can be reduced by adjusting the vocal-tract shape properly.     

An evaluation of the effects of three laryngeal lubricants on phonation threshold pressure (PTP)

Clinicians frequently offer advice to performers and voice-disordered patients aimed ostensibly to manipulate the water content and/or viscosity of the mucus blanket covering the vocal folds. To evaluate the relative effects of three potential laryngeal lubricants on phonatory function (ie, water, Mannitol--an osmotic agent, and Entertainer's Secret Throat Relief (Kli Corp., Carmel, IN)--a glycerin-based product), phonation threshold pressure (PTP) was measured in 18 healthy, vocally normal female participants twice before (baseline) and then four times after 2 ml of each substance were nebulized. PTP is the minimum subglottal pressure required to initiate vocal fold oscillation, and the lowering of PTP is assumed to correspond to physiologically more efficient phonation and reduced phonatory effort. Over a 3-week period, participants were tested on three separate occasions (at 1-week intervals). On each occasion, a different nebulized treatment was administered. PTP for both comfortable and high fundamental frequency productions was measured using an oral pressure-flow system (Perci-Sars, MicroTronics Corp., Chapel Hill, NC). Analysis of the results revealed that Mannitol, an agent that encourages osmotic water flux to the luminal airway surface, lowered PTP immediately after its administration (ie, p = 0.071, for high-pitched productions only). However, the duration of its PTP lowering effect was less than 20 minutes. The other two substances did not demonstrate any significant postadministration effect on PTP.     

Two-dimensional model of vocal fold vibration for sound synthesis of voice and soprano singing

Voiced sounds were simulated with a computer model of the vocal fold composed of a single mass vibrating both parallel and perpendicular to the airflow. Similarities with the two-mass model are found in the amplitudes of the glottal area and the glottal volume flow velocity, the variation in the volume flow waveform with the vocal tract shape, and the dependence of the oscillation amplitude upon the average opening area of the glottis, among other similar features. A few dissimilarities are also found in the more symmetric glottal and volume flow waveforms in the rising and falling phases. The major improvement of the present model over the two-mass model is that it yields a smooth transition between oscillations with an inductive load and a capacitive load of the vocal tract with no sudden jumps in the vibration frequency. Self-excitation is possible both below and above the first formant frequency of the vocal tract. By taking advantage of the wider continuous frequency range, the two-dimensional model can successfully be applied to the sound synthesis of a high-pitched soprano singing, where the fundamental frequency sometimes exceeds the first formant frequency.     

Spectrum warping based on sub-glottal resonances in speaker-independent speech recognition

To reduce degradation in speech recognition due to varied characteristics of different speakers,a method of perceptual frequency warping based on subglottal resonances for speaker normalization is proposed.The warping factor is extracted from the second subglottal resonance using acoustic coupling between subglottis and vocal tract.The second subglottal resonance is independent of the speech content,which reflects the speaker characteristics more than the third formant.The perceptual minimum variation distortionless response(PMVDR) coefficient is normalized,which is more robust and has better anti-noise capability than MFCC. The normalized coefficients are used in the speech-mode training and speech recognition.Experiments show that the word error rate,as compared with MFCC and the spectrum warping by the third formant,decreases by 4%and 3%respectively in clean speech recognition,and by 9%and 5%respectively in a noisy environment.The results indicate that the proposed method can improve the word recognition accuracy in a speaker-independent recognition system.