Comparison of biomechanical modeling of register transitions and voice instabilities with excised larynx experiments

Voice instabilities were studied using excised human larynx experiments and biomechanical modeling. With a controlled elongation of the vocal folds, the experiments showed registers with chest-like and falsetto-like vibrations. Observed instabilities included abrupt jumps between the two registers exhibiting hysteresis, aphonic episodes, subharmonics, and chaos near the register transitions. In order to model these phenomena, a three-mass model was constructed by adding a third mass on top of the simplified two-mass model. Simulation studies showed that the three-mass model can vibrate in both chest-like and falsetto-like patterns. Variation of tension parameters which mimic activities of laryngeal muscles could induce transitions between both registers. For reduced prephonatory areas and damping constants, extended coexistence of chest and falsetto registers was found, in agreement with experimental data. Subharmonics and deterministic chaos were observed close to transitions between the registers. It is concluded that the abrupt changes between chest and falsetto registers can be understood as shifts in dominance of eigenmodes of the vocal folds.     

Medial surface dynamics of an in vivo canine vocal fold during phonation

Quantitative measurement of the medial surface dynamics of the vocal folds is important for understanding how sound is generated within the larynx. Building upon previous excised hemilarynx studies, the present study extended the hemilarynx methodology to the in vivo canine larynx. Through use of an in vivo model, the medial surface dynamics of the vocal fold were examined as a function of active thyroarytenoid muscle contraction. Data were collected using high-speed digital imaging at a sampling frequency of 2000 Hz, and a spatial resolution of 1024 x 1024 pixels. Chest-like and fry-like vibrations were observed, but could not be distinguished based on the input stimulation current to the recurrent laryngeal nerve. The subglottal pressure did distinguish the registers, as did an estimate of the thyroarytenoid muscle activity. Upon quantification of the three-dimensional motion, the method of Empirical Eigenfunctions was used to extract the underlying modes of vibration, and to investigate mechanisms of sustained oscillation. Results were compared with previous findings from excised larynx experiments and theoretical models.     

Nonlinear dynamics of the voice: Signal analysis and biomechanical modeling

Irregularities in voiced speech are often observed as a consequence of vocal fold lesions, paralyses, and other pathological conditions. Many of these instabilities are related to the intrinsic nonlinearities in the vibrations of the vocal folds. In this paper, bifurcations in voice signals are analyzed using narrow-band spectrograms. We study sustained phonation of patients with laryngeal paralysis and data from an excised larynx experiment. These spectrograms are compared with computer simulations of an asymmetric 2-mass model of the vocal folds. (c) 1995 American Institute of Physics.     

Differences Among Mixed,Chest, and Falsetto Registers: A Multiparametric Study

IntroductionTypical singing registers are the chest and falsetto; however, trained singers have an additional register, namely, the mixed register. The mixed register, which is also called “mixed voice” or “mix,” is an important technique for singers, as it can help bridge from the chest voice to falsetto without noticeable voice breaks.ObjectiveThe present study aims to reveal the nature of the voice-production mechanism of the different registers (chest, mix, and falsetto) using high-speed digital imaging (HSDI), electroglottography (EGG), and acoustic and aerodynamic measurements.Study DesignCross-sectional study.MethodsAerodynamic measurements were acquired for twelve healthy singers (six men and women) during the phonation of a variety of pitches using three registers. HSDI and EGG devices were simultaneously used on three healthy singers (two men and one woman) from which an open quotient (OQ) and speed quotient (SQ) were detected. Audio signals were recorded for five sustained vowels, and a spectral analysis was conducted to determine the amplitude of each harmonic component. Furthermore, the absolute (not relative) value of the glottal volume flow was estimated by integrating data obtained from the HSDI and aerodynamic studies.ResultsFor all singers, the subglottal pressure (P_Sub) was the highest for the chest in the three registers, and the mean flow rate (MFR) was the highest for the falsetto. Conversely, the P_Sub of the mix was as low as the falsetto, and the MFR of the mix was as low as the chest. The HSDI analysis showed that the OQ differed significantly among the registers, even when the fundamental frequency was the same; the OQ of the mix was higher than that of the chest but lower than that of the falsetto. The acoustic analysis showed that, for the mix, the harmonic structure was intermediate between the chest and falsetto. The results of the glottal volume-flow analysis revealed that the maximum volume velocity was the least for the mix register at every fundamental frequency. The first and second harmonic (H1-H2) difference of the voice source spectrum was the greatest for the falsetto, then the mix, and finally, the chest.ConclusionsWe found differences in the registers in terms of the aeromechanical mechanisms and vibration patterns of the vocal folds. The mixed register proved to have a distinct voice-production mechanism, which can be differentiated from those of the chest or falsetto registers.     

Material parameter computation for multi-layered vocal fold models

Today, the prevention and treatment of voice disorders is an ever-increasing health concern. Since many occupations rely on verbal communication, vocal health is necessary just to maintain one's livelihood. Commonly applied models to study vocal fold vibrations and air flow distributions are self sustained physical models of the larynx composed of artificial silicone vocal folds. Choosing appropriate mechanical parameters for these vocal fold models while considering simplifications due to manufacturing restrictions is difficult but crucial for achieving realistic behavior. In the present work, a combination of experimental and numerical approaches to compute material parameters for synthetic vocal fold models is presented. The material parameters are derived from deformation behaviors of excised human larynges. The resulting deformations are used as reference displacements for a tracking functional to be optimized. Material optimization was applied to three-dimensional vocal fold models based on isotropic and transverse-isotropic material laws, considering both a layered model with homogeneous material properties on each layer and an inhomogeneous model. The best results exhibited a transversal-isotropic inhomogeneous (i.e., not producible) model. For the homogeneous model (three layers), the transversal-isotropic material parameters were also computed for each layer yielding deformations similar to the measured human vocal fold deformations.     

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.     

Nonlinear source-filter coupling in phonation: vocal exercises

Nonlinear source-filter coupling has been demonstrated in computer simulations, in excised larynx experiments, and in physical models, but not in a consistent and unequivocal way in natural human phonations. Eighteen subjects (nine adult males and nine adult females) performed three vocal exercises that represented a combination of various fundamental frequency and formant glides. The goal of this study was to pinpoint the proportion of source instabilities that are due to nonlinear source-tract coupling. It was hypothesized that vocal fold vibration is maximally destabilized when F(0) crosses F(1), where the acoustic load changes dramatically. A companion paper provides the theoretical underpinnings. Expected manifestations of a source-filter interaction were sudden frequency jumps, subharmonic generation, or chaotic vocal fold vibrations that coincide with F(0)-F(1) crossovers. Results indicated that the bifurcations occur more often in phonations with F(0)-F(1) crossovers, suggesting that nonlinear source-filter coupling is partly responsible for source instabilities. Furthermore it was observed that male subjects show more bifurcations in phonations with F(0)-F(1) crossovers, presumably because in normal speech they are less likely to encounter these crossovers as much as females and hence have less practice in suppressing unwanted instabilities.     

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.     

Medial Surface Dynamics as a Function of Subglottal Pressure in a Canine Larynx Model

During vocal fold vibration, there may be a mucosal wave in the superior-inferior (vertical) direction, resulting in a convergent shape during opening and a divergent shape during closing. Most of our understanding of the converging/diverging shape of the glottis has come from studies in a hemilarynx model. Previous work has shown that vibratory patterns in the full excised larynx are different than the hemilarynx. This study characterized the dynamics of the medial glottal wall geometry during vibrations in the full excised canine larynx model. Using particle image velocimetry, the intraglottal geometry was measured at the midmembranous coronal plane in an excised canine larynx model. Measurements of the glottal area were taken simultaneously using high-speed imaging. The results show that skewing of the glottal area waveform occurs without the presence of a vocal tract and that the phase-lag of the superior edge relative to the inferior edge is smaller than reported and depends on the subglottal pressure. In addition, it shows that the glottal divergence angle during closing is proportional to the magnitude of the acoustic intensity and the intraglottal negative pressure. This preliminary data suggests that more studies are needed to determine the important mechanisms determining the relationship between intraglottal flow, intraglottal geometry, and acoustics.     

Bifurcations in excised larynx experiments

Bifurcation analysis was applied to vocal fold vibration in excised larynx experiments. Phonation onset and vocal instabilities were studied in a parameter plane spanned by subglottal pressure and asymmetry of either vocal fold adduction or elongation. Various phonatory regimes were observed, including single vocal fold oscillations. Selected spectra demonstrated correspondence between these regimes and vocal registers noted in the literature. To illustrate the regions spanned by the various phonatory regimes, two-dimensional bifurcation diagrams were generated. Many instabilities or bifurcations were noted in the regions of coexistence, i.e., regions in which the phonatory regimes overlap. Bifurcations were illustrated with spectrograms and fundamental frequency contours. Where possible, results from these studies were related to clinical observations.     

Comparison between electroglottography and electromagnetic glottography

Newly developed glottographic sensors, utilizing high-frequency propagating electromagnetic waves, were compared to a well-established electroglottographic device. The comparison was made on four male subjects under different phonation conditions, including three levels of vocal fold adduction (normal, breathy, and pressed), three different registers (falsetto, chest, and fry), and two different pitches. Agreement between the sensors was always found for the glottal closure event, but for the general wave shape the agreement was better for falsetto and breathy voice than for pressed voice and vocal fry. Differences are attributed to the field patterns of the devices. Whereas the electroglottographic device can operate only in a conduction mode, the electromagnetic device can operate in either the forward scattering (diffraction) mode or in the backward scattering (reflection) mode. Results of our tests favor the diffraction mode because a more favorable angle imposed on receiving the scattered (reflected) signal did not improve the signal strength. Several observations are made on the uses of the electromagnetic sensors for operation without skin contact and possibly in an array configuration for improved spatial resolution within the glottis.     

Vocal tract area functions and formant frequencies in opera tenors' modal and falsetto registers

According to recent model investigations, vocal tract resonance is relevant to vocal registers. However, no experimental corroboration of this claim has been published so far. In the present investigation, ten professional tenors' vocal tract configurations were analyzed using MRI volumetry. All subjects produced a sustained tone on the pitch F4 (349 Hz) on the vowel /a/ (1) in modal and (2) in falsetto register. The area functions were estimated from the MRI data and their associated formant frequencies were calculated. In a second condition the same subjects repeated the same tasks in a sound treated room and their formant frequencies were estimated by means of inverse filtering. In both recordings similar formant frequencies were observed. Vocal tract shapes differed between modal and falsetto register. In modal as compared to falsetto the lip opening and the oral cavity were wider and the first formant frequency was higher. In this sense the presented results are in agreement with the claim that the formant frequencies differ between registers.     

Correspondence of electroglottographic closed quotient to vocal fold impact stress in excised canine larynges

The purpose of this study was to explore the possible use of the electroglottographic closed quotient (EGG CQ) as a noninvasive estimate of vocal fold impact stress (SI). Two excised canine larynges were used. Each larynx was mounted and vocal fold oscillation was induced using a humidified air source. Twentyseven experimental trials were conducted for each larynx. Trials involved variations in vocal process gap, vocal fold elongation, and subglottic pressure. Simultaneous measures were made of vocal fold SI at the midpoint of the membranous vocal folds, and EGG CQ (dimensionless ratio). The results indicated that when threshold and saturation effects were excluded, the SI and the CQ were strongly related (linear correlation r = .83 and .96 for the two individual larynges, and .81 for the combined data). Within the region of linear relation, an increase of .15 in the CQ corresponded to about 1 kPa increase in SI for the combined data. Discussion focuses on possible clinical implications and the likely reasons for threshold and saturation phenomena.     

Visualization and Quantification of the Medial Surface Dynamics of an Excised Human Vocal Fold During Phonation

SUMMARY: The purpose of this investigation was to investigate physical mechanisms of vocal fold vibration during normal phonation through quantification of the medial surface dynamics of the fold. An excised hemilarynx setup was used. The dynamics of 30 microsutures mounted on the medial surface of a human vocal fold were analyzed across 18 phonatory conditions. The vibrations were recorded with a digital high-speed camera at a frequency of 4,000 Hz. The positions of the sutures were extracted and converted to three-dimensional coordinates using a linear approximation technique. The data were reduced to principal eigenfuctions, which captured over 90% of the variance of the data, and suggested mechanisms of sustained vocal fold oscillation. The vibrations were imaged as the following phonatory conditions were manipulated: glottal airflow, an adductory force applied to the muscular process, and an elongation force applied to the thyroid cartilage. Over the range of variables studied, only the variation in glottal airflow yielded significant changes in subglottal pressure and fundamental frequency. All recordings showed high correlation for the distribution of the dynamics across the medial surface of the vocal fold. The distribution of the different displacement directions and velocities showed the highest variations around the superior region of the medial surface. Although the computed vibration patterns of the two largest empirical eigenfunctions were consistent with previous experimental observations, the relative prominence of the two eigenfunctions changed as a function of glottal airflow, impacting theories of vocal efficiency and vocal economy.     

The biomechanics of the medialization laryngoplasty(thyroplasty type 1) in an ex vivo canine model

The biomechanics of medialization laryngoplasty are not well understood. An excised canine larynx model was used to test the effects of various sized silicon implants. The vocal fold length, position, and tension were measured. Medialization laryngoplasty did not affect vocal fold length. At the mid-membranous vocal fold, larger shims resulted in greater medialization and tension. Medialization laryngoplasty neither medialized nor stiffened the vocal process to resist lateralizing forces. We conclude that medialization laryngoplasty provides bulk and support for defects of the membranous region of the vocal fold, but does not appear to close a posterior glottal gap. The selection of a surgical procedure to treat glottal incompetence should take into account the unique biomechanical properties of the anterior (membranous vocal folds) and posterior (cartilaginous portion) glottis.     

Three-dimensional biomechanical properties of human vocal folds: parameter optimization of a numerical model to match in vitro dynamics

The human voice signal originates from the vibrations of the two vocal folds within the larynx. The interactions of several intrinsic laryngeal muscles adduct and shape the vocal folds to facilitate vibration in response to airflow. Three-dimensional vocal fold dynamics are extracted from in vitro hemilarynx experiments and fitted by a numerical three-dimensional-multi-mass-model (3DM) using an optimization procedure. In this work, the 3DM dynamics are optimized over 24 experimental data sets to estimate biomechanical vocal fold properties during phonation. Accuracy of the optimization is verified by low normalized error (0.13 ± 0.02), high correlation (83% ± 2%), and reproducible subglottal pressure values. The optimized, 3DM parameters yielded biomechanical variations in tissue properties along the vocal fold surface, including variations in both the local mass and stiffness of vocal folds. That is, both mass and stiffness increased along the superior-to-inferior direction. These variations were statistically analyzed under different experimental conditions (e.g., an increase in tension as a function of vocal fold elongation and an increase in stiffness and a decrease in mass as a function of glottal airflow). The study showed that physiologically relevant vocal fold tissue properties, which cannot be directly measured during in vivo human phonation, can be captured using this 3D-modeling technique.     

Three registers in an untrained female singer analyzed by videokymography, strobolaryngoscopy and sound spectrography

There has been a lack of objective data on the singing voice registers, particularly on the so called "whistle" register, occurring in the top part of the female pitch range, which is accessible only to some singers. This study offers unique strobolaryngoscopic and high-speed (7812.5 imagess) videokymographic data on the vocal fold behavior of an untrained female singer capable of producing three distinct voice qualities, i.e., the chest, head and whistle registers. The sound was documented spectrographically. The transition from chest to head register, accompanied by pitch jumps, occurred around tones B4-C#5 (500-550 Hz) and was found to be associated with a slight decrease in arytenoids adduction, resulting in decrease of the closed quotient. The register shifts from head to whistle, also accompanied by pitch jumps, occurred around tones E5-B5 (670-1000 Hz) without any noticeable changes in arytenoids adduction. Some evidence was found for the vocal tract influence on this transition. The mechanism of the vocal fold vibration in whistle register was found principally similar to that at lower registers: vibrations along the whole glottal length and vertical phase differences (indicated by sharp lateral peaks in videokymography) were seen on the vocal folds up to the highest tone G6 (1590 Hz).     

Quantifying the complexity of excised larynx vibrations from high-speed imaging using spatiotemporal and nonlinear dynamic analyses

In this paper, we investigate the biomechanical applications of spatiotemporal analysis and nonlinear dynamic analysis to quantitatively describe regular and irregular vibrations of twelve excised larynges from high-speed image recordings. Regular vibrations show simple spatial symmetry, temporal periodicity, and discrete frequency spectra, while irregular vibrations show complex spatiotemporal plots, aperiodic time series, and broadband spectra. Furthermore, the global entropy and correlation length from spatiotemporal analysis and the correlation dimension from nonlinear dynamic analysis reveal a statistical difference between regular and irregular vibrations. In comparison with regular vibrations, the global entropy and correlation dimension of irregular vibrations are statistically higher, while the correlation length is significantly lower. These findings show that spatiotemporal analysis and nonlinear dynamic analysis are capable of describing the complex dynamics of vocal fold vibrations from high-speed imaging and may potentially be helpful for understanding disordered behaviors in biomedical laryngeal systems.