Geometry of human vocal folds and glottal channel for mathematical and biomechanical modeling of voice production

Current models of the vocal folds derive their shape from approximate information rather than from exactly measured data. The objective of this study was to obtain detailed measurements on the geometry of human vocal folds and the glottal channel in phonatory position. A non-destructive casting methodology was developed to capture the vocal fold shape from excised human larynges on both medial and superior surfaces. Two female larynges, each in two different phonatory configurations corresponding to low and high fundamental frequency of the vocal fold vibrations, were measured. A coordinate measuring machine was used to digitize the casts yielding 3D computer models of the vocal fold shape. The coronal sections were located in the models, extracted and fitted by piecewise-defined cubic functions allowing a mathematical expression of the 2D shape of the glottal channel. Left-right differences between the cross-sectional shapes of the vocal folds were found in both the larynges.     

The effect of three-dimensional glottal geometry on intraglottal quasi-steady flow distributions and their relationship with phonation

Vocal fold geometry plays an important role in human phonation. The intraglottal quasi-steady pressure and velocity distributions depend upon the shape, size, and diameter of the glottis. This study reports the effects of the variation of glottal shapes on intraglottal pressures and velocities using a Plexiglas model with a glottis having nine symmetric glottal angles (uniform, as well as convergent and divergent 5°, 10°, 20° and 40°), while the minimal glottal diameter was held constant at 0.06 cm. The empirical data were supported by penalty finite element computational results. The results suggest that larger convergent glottal angles correspond to increased pressures and decreased velocities in the glottis upstream of the minimum glottal location, with a reversal of this pattern at the minimal glottal diameter location. The pressure dip near the glottal entrance for divergent glottal angles was greatest for the 10° divergence angle condition, and was sequentially less for 5°, 20°, and 40°. Flow resistance was greater for a convergent angle than a divergent angle of the same value, and least for the 10° divergent condition. Pressure recovery in the glottis suggested that the optimal glottal diffuser angle was near 10°. Results suggest that the glottal geometry has a critical relationship with phonation (especially for vocal efficiency), and therefore important significance to understanding artistic voice and clinical voice management.     

Numerical analysis and comparison of flow fields in normal larynx and larynx with unilateral vocal fold paralysis

In this study, laryngeal flow fields are investigated and compared in normal larynx and models of larynx with unilateral vocal fold paralysis (UVFP). In paralytic models, three fixed initial glottal gaps are considered to understand the positive or probable negative impacts of surgical operation on unilaterally paralytic larynx, by which the paralyzed vocal fold is brought closer to the mid-plane. Various features of the flow fields have been discussed in detail including glottal gap width, glottal flow rate, glottal exit pressure pattern and glottal jet evolution. The numerical solution of fluid-structure interaction is carried out using ANSYS, and the results confirm some of the favorable effects of surgery on the patient’s larynx. It is also shown that by tightening the glottal gap, some of the problems caused by the presence of a motionless vocal fold, such as leakage through glottal gap in the closure phase resulting in breathy voice can be moderated, although some of the symptoms of this disorder remain relatively unchanged.     

The effect of three-dimensional glottal geometry on intraglottal quasi-steady flow distributions and their relationship with phonation

In order to model laryngeal aerodynamics from a quasi-steady point of view[1], both the dynamic distri-bution of intraglottal air pressures that act upon the vocal folds and the tissue properties of the vocal folds are required[2]. Concerning the first po…     

Comparative analysis of normal voice characteristics using simultaneous electroglottography and high speed digital imaging

Vocal fold (VF) activity in sustained phonation captured through simultaneous visual observations via direct method (high speed digital imaging: HSDI) and via indirect electroglottographic (EGG) method in normal population were submitted to contrastive analysis. The quantitative parameters used to describe the dynamic properties of the VF vibration included open/closed quotients, and indices that we defined as harmonic ratios. Our analysis also included a Nyquist plot approach that we introduced previously. Results indicate that EGG tracings delivered useful results for the analysis of normal voices and derived a reasonably accurate open quotient value in comparison to the visual analysis. Moreover, recognizable Nyquist plot patterns were obtained for EGG tracings although not identical to that obtained from the HSDI-derived glottal area waveform. Therefore ipso-facto we evaluate the usefulness of EGG compared to visual data in describing specific voice characteristics of normal speakers.     

Modeling mechanical stresses as a factor in the etiology of benign vocal fold lesions

Vocal fold tissue lesions such as nodules and polyps are thought to develop in response to mechanical stress that occurs during vocal fold collision. Two computational models of vocal fold collision during voice production are used to investigate this hypothesis. A one-dimensional lumped mass model, whose parameters are derived from vocal fold tissue dimensions and material properties, predicts stress perpendicular to the direction of impact (normal stress). A previously published three-dimensional finite element model that incorporates the same dimensions and properties predicts the entire stress tensor. The hypothesis is supported by predictions from the finite element model that three components of normal stress and one component of shear stress are increased during collision in the typical location of lesions (i.e. the center of the superior medial edge of the vocal fold in the middle of the vibrating and contact region). The lumped mass model predicts that mechanical stress is negatively correlated with mucosal thickness (increased by voice warm-up and hydration), is positively correlated with driving force (proportional to voice intensity), and is affected by voice production method. These relationships are consistent with clinical observations of vocal fold lesion risk factors and have implications for improving prevention and treatment of benign vocal fold lesions.     

Influence of asymmetric stiffness on the structural and aerodynamic response of synthetic vocal fold models

The influence of asymmetric vocal fold stiffness on voice production was evaluated using life-sized, self-oscillating vocal fold models with an idealized geometry based on the human vocal folds. The models were fabricated using flexible, materially-linear silicone compounds with Young's modulus values comparable to that of vocal fold tissue. The models included a two-layer design to simulate the vocal fold layered structure. The respective Young's moduli of elasticity of the “left” and “right” vocal fold models were varied to create asymmetric conditions. High-speed videokymography was used to measure maximum vocal fold excursion, vibration frequency, and left–right phase shift, all of which were significantly influenced by asymmetry. Onset pressure, a measure of vocal effort, increased with asymmetry. Particle image velocimetry (PIV) analysis showed significantly greater skewing of the glottal jet in the direction of the stiffer vocal fold model. Potential applications to various clinical conditions are mentioned, and suggestions for future related studies are presented.     

Adapted to roar: functional morphology of tiger and lion vocal folds

Vocal production requires active control of the respiratory system, larynx and vocal tract. Vocal sounds in mammals are produced by flow-induced vocal fold oscillation, which requires vocal fold tissue that can sustain the mechanical stress during phonation. Our understanding of the relationship between morphology and vocal function of vocal folds is very limited. Here we tested the hypothesis that vocal fold morphology and viscoelastic properties allow a prediction of fundamental frequency range of sounds that can be produced, and minimal lung pressure necessary to initiate phonation. We tested the hypothesis in lions and tigers who are well-known for producing low frequency and very loud roaring sounds that expose vocal folds to large stresses. In histological sections, we found that the Panthera vocal fold lamina propria consists of a lateral region with adipocytes embedded in a network of collagen and elastin fibers and hyaluronan. There is also a medial region that contains only fibrous proteins and hyaluronan but no fat cells. Young's moduli range between 10 and 2000 kPa for strains up to 60%. Shear moduli ranged between 0.1 and 2 kPa and differed between layers. Biomechanical and morphological data were used to make predictions of fundamental frequency and subglottal pressure ranges. Such predictions agreed well with measurements from natural phonation and phonation of excised larynges, respectively. We assume that fat shapes Panthera vocal folds into an advantageous geometry for phonation and it protects vocal folds. Its primary function is probably not to increase vocal fold mass as suggested previously. The large square-shaped Panthera vocal fold eases phonation onset and thereby extends the dynamic range of the voice.     

Mechanical stress during phonation in a self-oscillating finite-element vocal fold model

The stress information during phonation in the vocal folds is helpful in understanding the etiologies of vocal trauma and its related vocal diseases, such as nodules. In this paper, a self-oscillating finite-element model, which combines aerodynamic properties, tissue mechanics, airflow-tissue interactions, and vocal fold collisions, was used to simulate the vocal fold vibration during phonation. The spatial and temporal characteristics of mechanical stress in the vocal folds were predicted by this model. Temporally, it was found that mechanical stress periodically undulates with vibration of the vocal folds and that vocal fold impact causes a jump in the normal stress value. Spatially, the normal stress is significantly higher on the vocal fold surface than inside of the vocal folds. At the midpoint of the medial surface, the peak-to-peak amplitude of the normal stress reaches its maximum value. Using different lung pressures (0-1.5kPa) to drive the self-oscillating model, we found that lower lung pressure can effectively decrease the mechanical stress in the vocal folds. This study supports the fatigue damage hypothesis of vocal trauma. With this hypothesis and the numerical simulation in this study, the clinical observations of vocal fold trauma risk can be explained. This implies the mechanical stress predicted by this self-oscillating model could be valuable for predicting, preventing, and treating vocal fold injury.     

Songbirds use pulse tone register in two voices to generate low-frequency sound

The principal physical mechanism of sound generation is similar in songbirds and humans, despite large differences in their vocal organs. Whereas vocal fold dynamics in the human larynx are well characterized, the vibratory behaviour of the sound-generating labia in the songbird vocal organ, the syrinx, is unknown. We present the first high-speed video records of the intact syrinx during induced phonation. The syrinx of anaesthetized crows shows a vibration pattern of the labia similar to that of the human vocal fry register. Acoustic pulses result from short opening of the labia, and pulse generation alternates between the left and right sound sources. Spontaneously calling crows can also generate similar pulse characteristics with only one sound generator. Airflow recordings in zebra finches and starlings show that pulse tone sounds can be generated unilaterally, synchronously or by alternating between the two sides. Vocal fry-like dynamics therefore represent a common production mechanism for low-frequency sounds in songbirds. These results also illustrate that complex vibration patterns can emerge from the mechanical properties of the coupled sound generators in the syrinx. The use of vocal fry-like dynamics in the songbird syrinx extends the similarity to this unusual vocal register with mammalian sound production mechanisms.     

Predicting Achievable Fundamental Frequency Ranges in Vocalization Across Species

Vocal folds are used as sound sources in various species, but it is unknown how vocal fold morphologies are optimized for different acoustic objectives. Here we identify two main variables affecting range of vocal fold vibration frequency, namely vocal fold elongation and tissue fiber stress. A simple vibrating string model is used to predict fundamental frequency ranges across species of different vocal fold sizes. While average fundamental frequency is predominantly determined by vocal fold length (larynx size), range of fundamental frequency is facilitated by (1) laryngeal muscles that control elongation and by (2) nonlinearity in tissue fiber tension. One adaptation that would increase fundamental frequency range is greater freedom in joint rotation or gliding of two cartilages (thyroid and cricoid), so that vocal fold length change is maximized. Alternatively, tissue layers can develop to bear a disproportionate fiber tension (i.e., a ligament with high density collagen fibers), increasing the fundamental frequency range and thereby vocal versatility. The range of fundamental frequency across species is thus not simply one-dimensional, but can be conceptualized as the dependent variable in a multi-dimensional morphospace. In humans, this could allow for variations that could be clinically important for voice therapy and vocal fold repair. Alternative solutions could also have importance in vocal training for singing and other highly-skilled vocalizations.     

A computational study of systemic hydration in vocal fold collision

Mechanical stresses develop within vocal fold (VF) soft tissues due to phonation-associated vibration and collision. These stresses in turn affect the hydration of VF tissue and thus influence voice health. In this paper, high-fidelity numerical computations are described, taking into account fully 3D geometry, realistic tissue and air properties, and high-amplitude vibration and collision. A segregated solver approach is employed, using sophisticated commercial solvers for both the VF tissue and glottal airflow domains. The tissue viscoelastic properties were derived from a biphasic formulation. Two cases were considered, whereby the tissue viscoelastic properties corresponded to two different volume fractions of the fluid phase of the VF tissue. For each case, hydrostatic stresses occurring as a result of vibration and collision were investigated. Assuming the VF tissue to be poroelastic, interstitial fluid movement within VF tissue was estimated from the hydrostatic stress gradient. Computed measures of overall VF dynamics (peak airflow velocity, magnitude of VF deformation, frequency of vibration and contact pressure) were well within the range of experimentally observed values. The VF motion leading to mechanical stresses within the VFs and their effect on the interstitial fluid flux is detailed. It is found that average deformation and vibration of VFs tend to increase the state of hydration of the VF tissue, whereas VF collision works to reduce hydration.     

The effects of the false vocal fold gaps on intralaryngeal pressure distributions and their effects on phonation

Human phonation does not merely depend on the vibration of the vocal folds. Research by clinical and computer simulations has demonstrated that the false vocal fold (FVF) is an important laryngeal con-striction that plays a vital role during human voice production. This study explored the effects of the FVF gaps using both the three-dimensional Plexiglas model and the numerical computation methods. Twelve FVF gaps (ranging from 0.02 to 2.06 cm) were used in this study at three glottal angles (uniform and convergent/divergent 40°), two minimal glottal diameters (Dg) (0.04 cm and 0.06 cm) separately, and the constant subglottal pressure (8 cm H2O). The results suggested that (1) the intralaryngeal pressure was the lowest and the flow was the highest (least flow resistance) when the FVF gap was 1.5-2 times greater than Dg; (2) the divergent glottal angle gave lower pressure and greater flow than the conver-gent and uniform glottal angle as there were no FVF conditions; (3) the presence of the FVF decreased the effects of the glottal angle to a certain extent; and more importantly, (4) the presence of the FVF also moved the separation points downstream, straightened the glottal jet for a longer distance, decreased the overall laryngeal resistance, and reduced the energy dissipation, suggesting the significance of FVF in efficient voice production. These results may be incorporated in the phonatory models (physical or computational) for better understanding of vocal mechanics. The results might also be helpful in exploring the surgical and rehabilitative intervention of related voice problems.     

The anatomy of vocal divergence in North American Elk and European red deer

Loud and frequent vocalizations play an important role in courtship behavior in Cervus species. European red deer (Cervus elaphus) produce low‐pitched calls, whereas North American elk (Cervus canadensis) produce high‐pitched calls, which is remarkable for one of the biggest land mammals. Both species engage their vocal organs in elaborate maneuvers but the precise mechanism is unknown. Vocal organs were compared by macroscopic and microscopic dissection. The larynx is sexually dimorphic in red deer but not in elk. The laryngeal lumen is more constricted in elk, and narrows further during ontogeny. Several elements of the hyoid skeleton and two of four vocal tract segments are longer in red deer than in elk allowing greater vocal tract expansion and elongation. We conclude that elk submit the larynx and vocal tract to much higher tension than red deer, whereby, enormously stressed vocal folds of reduced effective length create a high resistance glottal source. The narrow, high impedance laryngeal vestibulum matches glottal and vocal tract impedance allowing maximum power transfer. In red deer longer and relaxed vocal folds create a less resistant glottal source and a wider vestibulum matches the low glottal impedance to the vocal tract, thereby also ensuring maximum power transfer. J. Morphol., 2013. © 2012 Wiley Periodicals, Inc.     

Investigating the Three-dimensional Flow Separation Induced by a Model Vocal Fold Polyp

The fluid-structure energy exchange process for normal speech has been studied extensively, but it is not well understood for pathological conditions. Polyps and nodules, which are geometric abnormalities that form on the medial surface of the vocal folds, can disrupt vocal fold dynamics and thus can have devastating consequences on a patient''s ability to communicate. Our laboratory has reported particle image velocimetry (PIV) measurements, within an investigation of a model polyp located on the medial surface of an in vitro driven vocal fold model, which show that such a geometric abnormality considerably disrupts the glottal jet behavior. This flow field adjustment is a likely reason for the severe degradation of the vocal quality in patients with polyps. A more complete understanding of the formation and propagation of vortical structures from a geometric protuberance, such as a vocal fold polyp, and the resulting influence on the aerodynamic loadings that drive the vocal fold dynamics, is necessary for advancing the treatment of this pathological condition. The present investigation concerns the three-dimensional flow separation induced by a wall-mounted prolate hemispheroid with a 2:1 aspect ratio in cross flow, i.e. a model vocal fold polyp, using an oil-film visualization technique. Unsteady, three-dimensional flow separation and its impact of the wall pressure loading are examined using skin friction line visualization and wall pressure measurements.     

Effects of Vocal Fold Nodules on Glottal Cycle Measurements Derived from High-Speed Videoendoscopy in Children

The goal of this study is to quantify the effects of vocal fold nodules on vibratory motion in children using high-speed videoendoscopy. Differences in vibratory motion were evaluated in 20 children with vocal fold nodules (5–11 years) and 20 age and gender matched typically developing children (5–11 years) during sustained phonation at typical pitch and loudness. Normalized kinematic features of vocal fold displacements from the mid-membranous vocal fold point were extracted from the steady-state high-speed video. A total of 12 kinematic features representing spatial and temporal characteristics of vibratory motion were calculated. Average values and standard deviations (cycle-to-cycle variability) of the following kinematic features were computed: normalized peak displacement, normalized average opening velocity, normalized average closing velocity, normalized peak closing velocity, speed quotient, and open quotient. Group differences between children with and without vocal fold nodules were statistically investigated. While a moderate effect size was observed for the spatial feature of speed quotient, and the temporal feature of normalized average closing velocity in children with nodules compared to vocally normal children, none of the features were statistically significant between the groups after Bonferroni correction. The kinematic analysis of the mid-membranous vocal fold displacement revealed that children with nodules primarily differ from typically developing children in closing phase kinematics of the glottal cycle, whereas the opening phase kinematics are similar. Higher speed quotients and similar opening phase velocities suggest greater relative forces are acting on vocal fold in the closing phase. These findings suggest that future large-scale studies should focus on spatial and temporal features related to the closing phase of the glottal cycle for differentiating the kinematics of children with and without vocal fold nodules.     

全喉与半喉模型声门区内准稳态流场分布研究

对于非对称声带发声过程的研究,有助于将正常语音的研究拓展到病理状态,从而为嗓音康复工程打下基础。采用具有嵌入式可活动声门结构的喉部物理模型,研究了声门最小直径为0．0402cm．跨声门压分别为100、500、1000和1500Pa时,全喉及半喉内的准稳态流场分布及其与发声参量的关系。同时．运用三维有限元方法预测了上述边界条件对应的流场分布,计算结果有效地支持了实验数据。结果表明,对称结构下存在着非对称压力和速度成分,但是由于它所占的比率有限(不超过10％),对正常发声的影响不大。非对称结构下,由声门入口处极高的上下表面压力差(通常为跨声门压的1-3倍)导致的倾斜流和涡流的出现、气流分离点位置后移、分离区域增大、声门出入口间压差占跨声门压的比率降低(平均30％)、压力速度场的变化程度减弱等因素以及由此带来的能量损耗,是非对称结构下发声效率降低、发音失真的主要原因。研究提示：声门重建方案的设计应尽量满足对称结构原则。     

Effect of the ventricular folds in a synthetic larynx model

Within the human larynx, the ventricular folds serve primarily as a protecting valve during swallowing. They are located directly above the sound-generating vocal folds. During normal phonation, the ventricular folds are passive structures that are not excited to periodical oscillations. However, the impact of the ventricular folds on the phonation process has not yet been finally clarified.An experimental synthetic human larynx model was used to investigate the effect of the ventricular folds on the phonation process. The model includes self-oscillating vocal fold models and allows the comparison of the pressure distribution at multiple locations in the larynx for configurations with and without ventricular folds.The results indicate that the ventricular folds increase the efficiency of the phonation process by reducing the phonation threshold level of the pressure below the vocal folds. Two effects caused by the ventricular folds could be identified as reasons: (1) a decrease in the mean pressure level in the region between vocal and ventricular folds (ventricles) and (2) an increase in the glottal flow resistance.The reason for the first effect is a reduction of the pressure level in the ventricles due to the jet entrainment and the low static pressure in the glottal jet. The second effect results from an increase in the glottal flow resistance that enhances the aerodynamic energy transfer into the vocal folds. This effect reduces the onset threshold of the pressure difference across the glottis.     

VOCAL NODULES—A High-Speed Photographic Analysis: Notes on Treatment

Vocal nodules are benign tumefactions on the vocal cords due to excessive or improper use of the voice. They vary in size, shape, location and histologic composition and are essentially a clinical rather than a pathologic entity.High-speed motion pictures at 5,000 frames per second revealed that they disturb the normal vibratory pattern in the same manner as benign tumors of the vocal cords generally. Treatment consists of vocal rest, vocal reeducation, and surgical operation, singly or in combination.