Simulation of vocal fold impact pressures with a self-oscillating finite-element model

Vocal fold impact pressures were studied using a self-oscillating finite-element model capable of simulating vocal fold vibration and airflow. The calculated airflow pressure is applied on the vocal fold as the driving force. The airflow region is then adjusted according to the calculated vocal fold displacement. The interaction between airflow and the vocal folds produces a self-oscillating solution. Lung pressures between 0.2 and 2.5 kPa were used to drive this self-oscillating model. The spatial distribution of the impact pressure was studied. Studies revealed that the tissue collision during phonation produces a very large impact pressure which correlates with the lung pressure and glottal width. Larger lung pressure and a narrower glottal width increase the impact pressure. The impact pressure was found to be roughly the square root of lung pressure. In the inferior-superior direction, the maximum impact pressure is related to the narrowest glottis. In the anterior-posteriorfirection, the greatest impact pressure appears at the midpoint of the vocal fold. The match between our numerical simulations and clinical observations suggests that this self-oscillating finite-element model might be valuable for predicting mechanical trauma of the vocal folds.     

Histologic and Rheologic Characterization of Vocal Fold Scarring

Scarring of the vocal fold causes considerable dysphonia and presents significant treatment challenges. A rabbit model was developed to investigate the histologic ultrastructure and rheologic properties of the scarred vocal fold lamina propria. Eleven rabbit larynges were scarred by means of forcep biopsy. Sixty days postoperatively, the rabbits were sacrificed and their vocal folds were harvested. Histological analysis of the scarred and normal lamina propria was completed for collagen, procollagen, elastin, and hyaluronic acid. Linear viscoelastic shear properties of the tissues were also measured, including elastic shear modulus and dynamic viscosity. Compared to normal vocal fold lamina propria, scarred tissues demonstrated significantly less collagen, an increase in procollagen, and a decrease in elastin. Rheologically, both elastic shear modulus and dynamic viscosity were significantly higher for the scarred tissues. Increased stiffness and viscosity do not appear to result from an increase in collagen, but rather appear to be related to the presence of new, disorganized collagen scaffolding. Results are interpreted in terms of the possible role of interstitial proteins in the etiology of increased stiffness and viscosity, which requires further investigation. This animal model should allow for systematic future investigations of vocal fold scarring and its treatment.     

Influence of supraglottal structures on the glottal jet exiting a two-layer synthetic, self-oscillating vocal fold model

A synthetic two-layer, self-oscillating, life-size vocal fold model was used to study the influence of the vocal tract and false folds on the glottal jet. The model vibrated at frequencies, pressures, flow rates, and amplitudes consistent with human phonation, although some differences in behavior between the model and the human vocal folds are noted. High-speed images of model motion and flow visualization were acquired. Phase-locked ensemble-averaged glottal jet velocity measurements using particle image velocimetry (PIV) were acquired with and without an idealized vocal tract, with and without false folds. PIV data were obtained with varying degrees of lateral asymmetric model positioning. Glottal jet velocity magnitudes were consistent with those measured using excised larynges. A starting vortex was observed in all test cases. The false folds interfered with the starting vortex, and in some cases vortex shedding from the false folds was observed. In asymmetric cases without false folds, the glottal jet tended to skew toward the nearest wall; with the false folds, the opposite trend was observed. rms velocity calculations showed the jet shear layer and laminar core. The rms velocities were higher in the vocal tract cases compared to the open jet and false fold cases.     

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.     

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.     

Characterization of chronic vocal fold scarring in a rabbit model

The purpose of the current study was to assess the histologic and rheologic properties of the scarred vocal fold lamina propria during a chronic phase of wound repair in a rabbit model. Eighteen rabbit larynges were scarred using a procedure that involved stripping the vocal fold lamina propria down to the thyroarytenoid muscle, using 3-mm microforceps. The approximate dimension of injury to the vocal fold was 3 x 1.5 x 0.5 mm [length x width x depth]. At 6 months postoperatively, histologic analysis of the scarred and control lamina propria in eight of these rabbits was completed for collagen, procollagen, elastin, and hyaluronic acid. Compared with control samples, scarred tissue samples revealed fragmented and disorganized elastin fibers. Additionally, collagen was significantly increased, organized, and formed thick bundles in the scarred vocal fold lamina propria. Measurements of the viscoelastic shear properties of the scarred and control lamina propria in the remaining 10 rabbits revealed increased elastic shear modulus (G') in 8 of 10 scarred samples and increased dynamic viscosity (eta') in 9 of 10 scarred samples. Although rheologic differences were not statistically significant, they revealed that on average, scarred samples were stiffer and more viscous than the normal controls. Histologic data are interpreted as indicating that by 6 months postinjury, the scarred rabbit vocal fold has reached a mature phase of wound repair, characterized by an increased, organized, and thick bundle collagen matrix. Rheologic data are interpreted as providing support for the potential role of increased, thick bundle collagen, and a disorganized elastin network on shear stiffness and dynamic viscosity in the chronic vocal fold scar. Based on these results, a 6-month postoperative time frame is proposed for future studies of chronic vocal fold scarring using the rabbit animal model.     

Measurements of vocal fold tissue viscoelasticity: approaching the male phonatory frequency range

Viscoelastic shear properties of human vocal fold tissues have been reported previously. However, data have only been obtained at very low frequencies (< or = 15 Hz). This necessitates data extrapolation to the frequency range of phonation based on constitutive modeling and time-temperature superposition. This study attempted to obtain empirical measurements at higher frequencies with the use of a controlled strain torsional rheometer, with a design of directly controlling input strain that introduced significantly smaller system inertial errors compared to controlled stress rheometry. Linear viscoelastic shear properties of the vocal fold mucosa (cover) from 17 canine larynges were quantified at frequencies of up to 50 Hz. Consistent with previous data, results showed that the elastic shear modulus (G'), viscous shear modulus (G"), and damping ratio (zeta) of the vocal fold mucosa were relatively constant across 0.016-50 Hz, whereas the dynamic viscosity (eta') decreased monotonically with frequency. Constitutive characterization of the empirical data by a quasilinear viscoelastic model and a statistical network model demonstrated trends of viscoelastic behavior at higher frequencies generally following those observed at lower frequencies. These findings supported the use of controlled strain rheometry for future investigations of the viscoelasticity of vocal fold tissues and phonosurgical biomaterials at phonatory frequencies.     

Voice source characteristics in Mongolian "throat singing" studied with high-speed imaging technique, acoustic spectra, and inverse filtering.

Mongolian "throat singing" can be performed in different modes. In Mongolia, the bass-type is called Kargyraa. The voice source in bass-type throat singing was studied in one male singer. The subject alternated between modal voice and the throat singing mode. Vocal fold vibrations were observed with high-speed photography, using a computerized recording system. The spectral characteristics of the sound signal were analyzed. Kymographic image data were compared to the sound signal and flow inverse filtering data from the same singer were obtained on a separate occasion. It was found that the vocal folds vibrated at the same frequency throughout both modes of singing. During throat singing the ventricular folds vibrated with complete but short closures at half the frequency of the true vocal folds, covering every second vocal fold closure. Kymographic data confirmed the findings. The spectrum contained added subharmonics compared to modal voice. In the inverse filtered signal the amplitude of every second airflow pulse was considerably lowered. The ventricular folds appeared to modulate the sound by reducing the glottal flow of every other vocal fold vibratory cycle.     

Laser Doppler Measurements of the Vocal Fold Blood Micro-Circulation

The authors have developed a laser Doppler method, referred to as LDF+LS (laser Doppler flowmetry by laryngostereometry), for measuring the vocal fold micro-circulation. The vocal folds of 103 patients were examined during general anesthesia. A laser Doppler probe was placed on defined positions at the vocal fold edge: midmembranous position (MP), 2 mm and 4 mm behind, and 1 mm anterior to MP. Three parameters were derived, ie, the concentration (CMBC) and velocity (V) of moving blood cells, and the product of both resulting in perfusion: P = CMBCxV. The results of 53 subjects with normal vocal fold status showed that the V at MP was significantly lower than 2 mm behind MP (P < 0.05). Men had significantly higher velocity than women (P < 0.05). Subjects older than 65 years had higher perfusion and CMBC as compared with younger subjects (P < 0.01). Smokers with normal vocal folds had higher velocity than non-smokers (P < 0.05). Measurements of healthy vocal folds were compared with benign vocal fold pathology. Vocal fold polyps had significantly higher perfusion than normal vocal folds (P < 0.01). In conclusion, the laser Doppler measurement is a sensitive tool estimating superficial blood flow in both normal and pathological vocal folds. The blood flow in normals differs with respect to gender and age probably due to variations in micro-circulation of the mucosa and lamina propria. Vocal fold pathology and external factors such as smoking result in blood flow changes.     

Vocal Fold Impact Stress Analysis

Vocal fold impact stress (force/area) has been implicated as a factor possibly contributing to the formation of nodules and polyps. The force of impact of a moving body is related to its acceleration. Since the mass of the folds is relatively constant, one expects impact force to be directly proportional to acceleration. A measure that reflects the relative displacement of the vocal folds is photoglottography (PGG). The velocity and acceleration of the folds are easily obtained by calculating the first and second derivatives of the PGG displacement waveform. This study, therefore, compared the second derivative of the PGG signal with simultaneously measured impact stress in an excised canine larynx model. Glottal transillumination (PGG) was measured with a subglottic transducer. A miniature force transducer placed in the midline between the vocal folds measured impact stress at the midglottal position. For nine different larynges, there was a positive and linear relationship between the second derivative of PGG and impact stress. The statistically significant results support the hypothesis that the second derivative of PGG m ay provide a use fulnoninvasive way to estimate relative vocal fold impact stress.     

Assessment of local vocal fold deformation characteristics in an in vitro static tensile test

Voice quality is strongly dependent on vocal fold dynamics, which in turn are dependent on lung pressure and vocal fold biomechanics. Numerical and physical models are often used to investigate the interactions of these different subsystems. However, the utility of numerical and physical models is limited unless appropriately validated with data from physiological models. Hence a method that enables analysis of local vocal fold deformations along the entire surface is presented. In static tensile tests, forces are applied to distinctive working points being located in cover and muscle, respectively, so that specific layer properties can be investigated. The forces are directed vertically upward and are applied along or above the vocal fold edge. The resulting deformations are analyzed using multiple perspectives and three-dimensional reconstruction. Deformation characteristics of four human vocal folds were investigated. Preliminary results showed two phases of deformation: a range with a small slope for small deformations fading into a significant nonlinear deformation trend with a high slope. An increase of tissue stiffness from posterior to anterior was detected. This trend is more significant for muscle and in the mid-anterior half of the vocal fold.     

Voice Source Characteristics in Mongolian “Throat Singing” Studied with High-Speed Imaging Technique, Acoustic Spectra, and Inverse Filtering

Mongolian “throat singing” can be performed in different modes. In Mongolia, the bass-type is called Kargyraa. The voice source in bass-type throat singing was studied in one male singer. The subject alternated between modal voice and the throat singing mode. Vocal fold vibrations were observed with high-speed photography, using a computerized recording system. The spectral characteristics of the sound signal were analyzed. Kymographic image data were compared to the sound signal and flow inverse filtering data from the same singer were obtained on a separate occasion. It was found that the vocal folds vibrated at the same frequency throughout both modes of singing. During throat singing the ventricular folds vibrated with complete but short closures at half the frequency of the true vocal folds, covering every second vocal fold closure. Kymographic data confirmed the findings. The spectrum contained added subharmonics compared to modal voice. In the inverse filtered signal the amplitude of every second airflow pulse was considerably lowered. The ventricular folds appeared to modulate the sound by reducing the glottal flow of every other vocal fold vibratory cycle.     

Treatment of sulcus vocalis: Auditory perceptual and acoustical analysis of the slicing mucosa surgical technique

Sulcus vocalis is a furrow in the mucosa of the vocal folds that causes glottic insufficiency, irregular vibrations, and a distinctive dysphonia of variable severity. Treatment of sulcus vocalis needs to achieve anatomical and functional improvements that satisfy the behavior of the larynx and vocal quality. We suggest a combined approach, which includes vocal fold mucosal surgery through a technique called slicing mucosa, followed by intensive vocal rehabilitation. The aim of the surgery is to detach the mucosa of the sulcus and to interrupt the longitudinal fibrotic tension lines. The results of this procedure in 10 patients were analyzed through visual, auditory, and acoustical methods, and the data lead us to conclude that the treatment proposed offered a good anatomical result, a better vibratory pattern, and an improvement in overall vocal quality.     

Computation of physiological human vocal fold parameters by mathematical optimization of a biomechanical model

With the use of an endoscopic, high-speed camera, vocal fold dynamics may be observed clinically during phonation. However, observation and subjective judgment alone may be insufficient for clinical diagnosis and documentation of improved vocal function, especially when the laryngeal disease lacks any clear morphological presentation. In this study, biomechanical parameters of the vocal folds are computed by adjusting the corresponding parameters of a three-dimensional model until the dynamics of both systems are similar. First, a mathematical optimization method is presented. Next, model parameters (such as pressure, tension and masses) are adjusted to reproduce vocal fold dynamics, and the deduced parameters are physiologically interpreted. Various combinations of global and local optimization techniques are attempted. Evaluation of the optimization procedure is performed using 50 synthetically generated data sets. The results show sufficient reliability, including 0.07 normalized error, 96% correlation, and 91% accuracy. The technique is also demonstrated on data from human hemilarynx experiments, in which a low normalized error (0.16) and high correlation (84%) values were achieved. In the future, this technique may be applied to clinical high-speed images, yielding objective measures with which to document improved vocal function of patients with voice disorders.     

Dependence of phonation threshold pressure on vocal tract acoustics and vocal fold tissue mechanics

Analytical and computer simulation studies have shown that the acoustic impedance of the vocal tract as well as the viscoelastic properties of vocal fold tissues are critical for determining the dynamics and the energy transfer mechanism of vocal fold oscillation. In the present study, a linear, small-amplitude oscillation theory was revised by taking into account the propagation of a mucosal wave and the inertive reactance (inertance) of the supraglottal vocal tract as the major energy transfer mechanisms for flow-induced self-oscillation of the vocal fold. Specifically, analytical results predicted that phonation threshold pressure (Pth) increases with the viscous shear properties of the vocal fold, but decreases with vocal tract inertance. This theory was empirically tested using a physical model of the larynx, where biological materials (fat, hyaluronic acid, and fibronectin) were implanted into the vocal fold cover to investigate the effect of vocal fold tissue viscoelasticity on Pth. A uniform-tube supraglottal vocal tract was also introduced to examine the effect of vocal tract inertance on Pth. Results showed that Pth decreased with the inertive impedance of the vocal tract and increased with the viscous shear modulus (G") or dynamic viscosity (eta') of the vocal fold cover, consistent with theoretical predictions. These findings supported the potential biomechanical benefits of hyaluronic acid as a surgical bioimplant for repairing voice disorders involving the superficial layer of the lamina propria, such as scarring, sulcus vocalis, atrophy, and Reinke's edema.     

Vocal fold and ventricular fold vibration in period-doubling phonation: physiological description and aerodynamic modeling

Occurrences of period-doubling are found in human phonation, in particular for pathological and some singing phonations such as Sardinian A Tenore Bassu vocal performance. The combined vibration of the vocal folds and the ventricular folds has been observed during the production of such low pitch bass-type sound. The present study aims to characterize the physiological correlates of this acoustical production and to provide a better understanding of the physical interaction between ventricular fold vibration and vocal fold self-sustained oscillation. The vibratory properties of the vocal folds and the ventricular folds during phonation produced by a professional singer are analyzed by means of acoustical and electroglottographic signals and by synchronized glottal images obtained by high-speed cinematography. The periodic variation in glottal cycle duration and the effect of ventricular fold closing on glottal closing time are demonstrated. Using the detected glottal and ventricular areas, the aerodynamic behavior of the laryngeal system is simulated using a simplified physical modeling previously validated in vitro using a larynx replica. An estimate of the ventricular aperture extracted from the in vivo data allows a theoretical prediction of the glottal aperture. The in vivo measurements of the glottal aperture are then compared to the simulated estimations.     

A simple-shear rheometer for linear viscoelastic characterization of vocal fold tissues at phonatory frequencies

Previous studies reporting the linear viscoelastic shear properties of the human vocal fold cover or mucosa have been based on torsional rheometry, with measurements limited to low audio frequencies, up to around 80 Hz. This paper describes the design and validation of a custom-built, controlled-strain, linear, simple-shear rheometer system capable of direct empirical measurements of viscoelastic shear properties at phonatory frequencies. A tissue specimen was subjected to simple shear between two parallel, rigid acrylic plates, with a linear motor creating a translational sinusoidal displacement of the specimen via the upper plate, and the lower plate transmitting the harmonic shear force resulting from the viscoelastic response of the specimen. The displacement of the specimen was measured by a linear variable differential transformer whereas the shear force was detected by a piezoelectric transducer. The frequency response characteristics of these system components were assessed by vibration experiments with accelerometers. Measurements of the viscoelastic shear moduli (G' and G") of a standard ANSI S2.21 polyurethane material and those of human vocal fold cover specimens were made, along with estimation of the system signal and noise levels. Preliminary results showed that the rheometer can provide valid and reliable rheometric data of vocal fold lamina propria specimens at frequencies of up to around 250 Hz, well into the phonatory range.     

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.     

Vestibular vocal fold behavior during phonation in unilateral vocal fold paralysis

Voice quality in patients with vocal fold paralysis can be affected by several factors, such as the position of the paralyzed vocal fold, its degree of atrophy, the configuration of its free edge, and the level differences between both vocal folds. Depending on the related vocal deficiency the patient will attempt to compensate using different maneuvers, such as increment of vocal tract and neck muscle contraction to improve glottal closure. This is probably one of the reasons why ventricular folds are frequently requested. The objective of this study is to analyze the behavior of the homolateral and contralateral vestibular folds to delineate patterns of vestibular motion during sustained phonation, in cases of unilateral vocal fold paralysis.