Finite element modeling of the effects of velopharyngeal opening on vocal tract reactance in female voice.

Classical singers use nasal consonants as "resonance exercises," and experimental results have shown that singers may use some velopharyngeal opening (VPO), most often in [a:] and more seldom in [i:] or [u:]. In particular, male singers have been found to increase VPO as pitch rises toward register change frequencies (passaggi). Laryngoscopic findings have shown that some VPO stabilizes vocal fold vibration; the effect is related to positive reactance. This study investigates the effects of VPO on vocal tract (VT) reactance over the range of fundamental frequencies (f0) used in singing using a computerized tomography-based finite element model of the VT and nose of a female. According to the results, by raising the lowest VT resonances, the VPO increased the VT reactance in the frequency ranges 207-359 Hz for [i:], 265-411 Hz for [u:], and 500-611 Hz for [a:], depending on the VPO size (full or half VPO). These frequency ranges are close to the first and second passaggio of a female singer. The change may have an especially practical stabilizing effect for [a:], which is otherwise characterized by very large changes in VT reactance and negative reactance at the second passaggio.

Influence of nasal cavities on voice quality: Computer simulations and experiments.

Nasal cavities are known to introduce antiresonances (dips) in the sound spectrum reducing the acoustic power of the voice. In this study, a three-dimensional (3D) finite element (FE) model of the vocal tract (VT) of one female subject was created for vowels [a:] and [i:] without and with a detailed model of nasal cavities based on CT (Computer Tomography) images. The 3D FE models were then used for analyzing the resonances, antiresonances and the acoustic pressure response spectra of the VT. The computed results were compared with the measurements of a VT model for the vowel [a:], obtained from the FE model by 3D printing. The nasality affects mainly the lowest formant frequency and decreases its peak level. The results confirm the main effect of nasalization, i.e., that sound pressure level decreases in the frequency region of the formants F1-F2 and emphasizes the frequency region of the formants F3-F5 around the singer's formant cluster. Additionally, many internal local resonances in the nasal and paranasal cavities were found in the 3D FE model. Their effect on the acoustic output was found to be minimal, but accelerometer measurements on the walls of the 3D-printed model suggested they could contribute to structure vibrations.

Glottal Imaging Study Comparing Vowel Phonation with Semioccluded Vocal Tract Exercises (SOVTEs) Produced Loudly.

OBJECTIVES: This study aimed to estimate vocal loading in loud phonation of a vowel and two widely used semiocclusion voice exercises (SOVTEs). Impact stress (IS) was estimated from glottal closing speed, inertial forces from the second derivative of glottal opening and closing. STUDY DESIGN: Experimental study in vivo. METHODS: A vocally healthy male sustained the [o:] vowel with habitual loudness and loudly: (1) without a tube, (2) into a silicone "Lax Vox" type tube (35 cm in length, 10 mm in diameter) outer end submerged 10 cm in water, and (3) into a straw (length 12.6 cm, diameter 2.5 mm) the outer end in air. He tried to use equal effort in all loud samples. High-speed video-laryngo-endoscopy was performed with a rigid scope. Oral air pressure (Poral) was registered in a mouthpiece through which the endoscope was inserted into the larynx and to which the tubes were attached air-tightly. RESULTS: Compared with vowel phonation at habitual loudness, mean of maximal glottal width (max GW) increased by 44.1% for loud tube phonation and decreased by 1.8% for loud straw phonation, and mean absolute value of minimum GW time derivative dmin increased by 57.1% for tube and by 29.5% for straw suggesting faster glottal closing. Compared with loud vowel phonation, max GW increased by 22.6% for loud tube phonation, while it decreased by 16.6% for loud straw phonation. For the tube, dmindecreased by 7.6% and for the straw by 23.8%. Maximal acceleration (ACC) and deceleration (DC) values were larger for the tube and smaller for the straw than the values for both vowel phonations. CONCLUSIONS: IS, deduced from dmin, increased in loud SOVTEs compared to vowel phonation at a conversational loudness, but remained lower in loud SOVTEs than in loud vowel phonation, particularly with a narrow straw, which also reduced inertial forces, as suggested by the reduced ACC and DC.

High-Speed Image Analysis Comparing Loading of Vocal Folds During Coughing and Phonation: A Case Study.

PURPOSE: Coughing is related to voice problems, since it involves firm glottal closure, fast glottal opening, and high subglottic pressure and flow rate. In this study, the glottal area variation and movements of laryngeal structures during coughing and phonation are compared. METHODS: High-speed laryngoscopy recordings were made of a normophonic male participant with a healthy larynx producing a neutral vowel and coughing. Oral air pressure was registered in a mouthpiece, through which an endoscope was inserted into the pharynx. Electroglottography, acoustic, and pressure signals were recorded simultaneously. The glottal width variation at the membranous and cartilaginous parts of the glottis was derived from the high-speed images, and the strong vibration of the false vocal folds was also registered. RESULTS: In coughing, compared to ordinary vowel phonation in nearly the same sound pressure level (93-94 dB6cm), the glottal width was 25% larger at the middle of the vocal folds, the maximum glottal opening velocity was 39% higher, and the maximum glottal width declination rate during glottal closing was up to three times higher. The maximum acceleration was 40% higher, and the maximum deceleration was 47% higher. Fundamental frequency f0 was the highest (ca. 400 Hz) at the beginning of the first phase of a typical coughing process. During the last part of the coughing process, f0 decreased from ca. 250 Hz to ca. 85 Hz at the phonation offset. CONCLUSIONS: The remarkable increase in maximum glottal width declination rate implies much higher vocal fold loading in coughing compared to phonation.

Vocal tract changes caused by phonation into a tube: a case study using computer tomography and finite-element modeling.

Phonation into a glass tube is a voice training and therapy method that leads to beneficial effects in voice production. It has not been known, however, what changes occur in the vocal tract during and after the phonation into a tube. This pilot study examined the vocal tract shape in a female subject before, during, and after phonation into a tube using computer tomography (CT). Three-dimensional finite-element models (FEMs) of the vocal tract were derived from the CT images and used to study changes in vocal tract input impedance. When phonating on vowel [a:] the data showed tightened velopharyngeal closure and enlarged cross-sectional areas of the oropharyngeal and oral cavities during and after the tube-phonation. FEM calculations revealed an increased input inertance of the vocal tract and an increased acoustic energy radiated out of the vocal tract after the tube-phonation. The results indicate that the phonation into a tube causes changes in the vocal tract which remain also when the tube is removed. These effects may help improving voice production in patients and voice professionals.

How Much Loading Does Water Resistance Voice Therapy Impose on the Vocal Folds? An Experimental Human Study.

OBJECTIVES: Water resistance voice therapy applies phonation into water through a tube. This study investigates how strenuous this therapy can be for the vocal folds in terms of impact stress (IS). It further examines whether it is possible to estimate the IS using the contact quotient (CQ) and maximum derivative from an electroglottogram (EGG). STUDY DESIGN: Experimental study. METHODS: A male participant sustained a rounded back vowel [u:] or [o:] at a comfortable speaking pitch and loudness, and phonated into a silicone "Lax Vox" tube submerged 2 cm in water. High-speed videolaryngoscopy was performed with a rigid scope. Oral air pressure (Poral) was registered in a mouthpiece through which an endoscope was inserted into the larynx. An EGG was recorded. RESULTS: The CQEGG from the EGG and the closed quotient from the glottal width (CQarea) increased, while the maximum glottal amplitude and absolute value of derivative minimum (dmin) and also the derivative maximum from the EGG decreased for phonation into water. Normalized amplitude quotient from the glottal width variation also decreased but the change was not significant. CONCLUSIONS: Based on the glottal area findings, water resistance therapy does not seem to increase vocal fold loading (in terms of increased IS) even if the increase of CQarea, and CQEGG suggest so. CQEGG may qualitatively correspond to that of area, but the reliability of CQ (from the glottal area or the EGG) and the maximum derivative from the EGG as estimates of IS in semiocclusion exercises warrant further studies.

Comparison of acceleration and impact stress as possible loading factors in phonation: a computer modeling study.

Impact stress (the impact force divided by the contact area of the vocal folds) has been suspected to be the main traumatizing mechanism in voice production, and the main cause of vocal fold nodules. However, there are also other factors, such as the repetitive acceleration and deceleration, which may traumatize the vocal fold tissues. Using an aeroelastic model of voice production, the present study quantifies the acceleration and impact stress values in relation to lung pressure, fundamental frequency (F0) and prephonatory glottal half-width. Both impact stress and acceleration were found to increase with lung pressure. Compared to impact stress, acceleration was less dependent on prephonatory glottal width and, thus, on voice production type. Maximum acceleration values were about 5-10 times greater for high F0 (approx. 400 Hz) compared to low F0 (approx. 100 Hz), whereas maximum impact stress remained nearly unchanged. This suggests that acceleration, i.e. the inertia forces, may present at high F0 a greater load for the vocal folds, and in addition to the collision forces may contribute to the fact that females develop vocal fold nodules and other vocal fold traumas more frequently than males.

Impact Stress in Water Resistance Voice Therapy: A Physical Modeling Study.

OBJECTIVES: Phonation through a tube in water is used in voice therapy. This study investigates whether this exercise may increase mechanical loading on the vocal folds. STUDY DESIGN: This is an experimental modeling study. METHODS: A model with three-layer silicone vocal fold replica and a plexiglass, MK Plexi, Prague vocal tract set for the articulation of vowel [u:] was used. Impact stress (IS) was measured in three conditions: for [u:] (1) without a tube, (2) with a silicon Lax Vox tube (35 cm in length, 1 cm in inner diameter) immersed 2 cm in water, and (3) with the tube immersed 10 cm in water. Subglottic pressure and airflow ranges were selected to correspond to those reported in normal human phonation. RESULTS: Phonation threshold pressure was lower for phonation into water compared with [u:] without a tube. IS increased with the airflow rate. IS measured in the range of subglottic pressure, which corresponds to measurements in humans, was highest for vowel [u:] without a tube and lower with the tube in water. CONCLUSIONS: Even though the model and humans cannot be directly compared, for instance due to differences in vocal tract wall properties, the results suggest that IS is not likely to increase harmfully in water resistance therapy. However, there may be other effects related to it, possibly causing symptoms of vocal fatigue (eg, increased activity in the adductors or high amplitudes of oral pressure variation probably capable of increasing stress in the vocal fold). These need to be studied further, especially for cases where the water bubbling frequency is close to the acoustical-mechanical resonance and at the same time the fundamental phonation frequency is near the first formant frequency of the system.

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.

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.

Estimation of impact stress using an aeroelastic model of voice production.

The maximum impact stress at the contact of the vocal folds achieved during the oscillation cycle was estimated in phonation using an aeroelastic model of voice production. Relations of impact stress to the lung pressure, fundamental frequency of self-oscillations, prephonatory glottal width, sound pressure level generated at the end of the glottis and vibration amplitude of the vocal folds were studied. Using the fundamental frequency, prephonatory glottal width, lung pressure and airflow rate values found in normal speech, maximum impact stress values of 2-3 kPa were obtained. The results fall well within the limits reported for excised canine larynges and human subjects. Impact stress increased with lung pressure almost linearly after the phonation threshold but reached a plateau when the limit of maximum glottal opening was reached. When the fundamental frequency and lung pressure were kept constant, impact stress appears to fit closely with a parabolic function of prephonatory glottal half-width.

Resonance Tube or Lax Vox?

OBJECTIVE: This study compares the flow resistance of Resonance tube (RT) and Lax Vox tube (LVT) when submerged 2 cm and 10 cm in water, as well as phonation into the tubes in these conditions. METHODS: In the in vitro experiment, the air pressure for flow rates of 60-600 mL/s was measured at the tube inlet, when the outer end of the tube was submerged 2 cm and 10 cm below water surface in 30°, 45°, and 90° angle. In the in vivo experiment, 14 subjects phonated in habitual loudness and loudly into both tubes, with the outer end 2 cm and 10 cm in water. RT was immersed in a 45° angle and LVT in a 90° angle in water. Oral pressure, contact quotient from electroglottographic signal, and sound pressure level were studied. Sensations during phonation were reported in an interview. RESULTS: Flow resistance was slightly lower with LVT than with RT, and slightly lower for smaller immersion angles. In habitual loudness, transglottic pressure and frequency of oral pressure variation were lower for LVT phonation and amplitude of oral pressure variation was higher for LVT 2 cm in water. Some subjects preferred RT, whereas others preferred LVT or reported no differences between them. CONCLUSIONS: The tubes differed slightly in flow resistance. Higher oral pressure oscillation with LVT 2 cm in water may offer stronger massage effect on vocal folds.

Computational modelling of effect of tonsillectomy on voice production.

The aim of this study is to develop mathematical 3D finite element (FE) models for numerical simulations of vibroacoustic properties of the human vocal tract after a tonsillectomy. Similar experimental studies are not easily realizable on living subjects. The FE models of the acoustic spaces corresponding to the human vocal tract for the Czech vowels /a:/ and /i:/ and the acoustic space around the human head were used in numerical simulations of phonation. The acoustic resonant characteristics of the FE models were studied using modal and transient analyses (excitation by a short pulse). Calculated results show that a tonsillectomy causes a frequency shift of the 3rd (down by approximately 180 Hz) and 4th (down by approximately 120 Hz) formants down to the lower frequencies for the vowel /a:/ and similarly for the 2nd, 4th and 5th formants for the vowel /i:/ (all down by approximately 100 Hz). Similar effects and results can be found in experimental studies in literature. The formant changes are dependent on the size of additional acoustic spaces that occur after a tonsillectomy. The verification of the model was performed on the recordings of patients before and after the tonsillectomy operation. Multi-Dimensional Voice Program (MDVP Advanced, KAY Elemetrics Corp.) was used for the comparing of the formant centre frequencies. Very small differences in the results of subjective evaluation of the voice before and after tonsillectomy were found.

Human vocal tract resonances and the corresponding mode shapes investigated by three-dimensional finite-element modelling based on CT measurement.

Resonance frequencies of the vocal tract have traditionally been modelled using one-dimensional models. These cannot accurately represent the events in the frequency region of the formant cluster around 2.5-4.5 kHz, however. Here, the vocal tract resonance frequencies and their mode shapes are studied using a three-dimensional finite element model obtained from computed tomography measurements of a subject phonating on vowel [a:]. Instead of the traditional five, up to eight resonance frequencies of the vocal tract were found below the prominent antiresonance around 4.7 kHz. The three extra resonances were found to correspond to modes which were axially asymmetric and involved the piriform sinuses, valleculae, and transverse vibrations in the oral cavity. The results therefore suggest that the phenomenon of speaker's and singer's formant clustering may be more complex than originally thought.

Vocal Fold Adjustment Caused by Phonation Into a Tube: A Double-Case Study Using Computed Tomography.

OBJECTIVES: Phonation into a tube is a widely used method for vocal training and therapy. Previous studies and practical experience show that the phonation becomes easier and louder after such an exercise. The purpose of this study was to find out whether there are systematic changes in the vocal fold adjustment after the exercise. METHODS: Two volunteer subjects (1 male and 1 female) without voice disorders were examined with computed tomography (CT). Both produced a sustained vowel [a:] at comfortable pitch and loudness before and after the tube phonation and a vowel-like phonation into the tube. Computed tomography (CT) scans were obtained before, during, and after the exercise, twice for each condition. The gathered CT images were used for measurements of vertical vocal fold thickness, bulkiness, length, and glottal width. RESULTS: No prominent trends common to both subjects were found in vocal fold adjustment during and after the phonation into the tube. Variability observed under the same conditions was usually of the same magnitude as the changes before and after the tube phonation. CONCLUSIONS: Changes in vocal tract configuration observed after the resonance tube exercises in previous related studies were more prominent than the changes in vocal fold configuration observed here.

Vocal tract and glottal function during and after vocal exercising with resonance tube and straw.

OBJECTIVE: The present study aimed to investigate the vocal tract and glottal function during and after phonation into a tube and a stirring straw. METHODS: A male classically trained singer was assessed. Computerized tomography (CT) was performed when the subject produced [a:] at comfortable speaking pitch, phonated into the resonance tube and when repeating [a:] after the exercise. Similar procedure was performed with a narrow straw after 15 minutes silence. Anatomic distances and area measures were obtained from CT midsagittal and transversal images. Acoustic, perceptual, electroglottographic (EGG), and subglottic pressure measures were also obtained. RESULTS: During and after phonation into the tube or straw, the velum closed the nasal passage better, the larynx position lowered, and hypopharynx area widened. Moreover, the ratio between the inlet of the lower pharynx and the outlet of the epilaryngeal tube became larger during and after tube/straw phonation. Acoustic results revealed a stronger spectral prominence in the singer/speaker's formant cluster region after exercising. Listening test demonstrated better voice quality after straw/tube than before. Contact quotient derived from EGG decreased during both tube and straw and remained lower after exercising. Subglottic pressure increased during straw and remained somewhat higher after it. CONCLUSION: CT and acoustic results indicated that vocal exercises with increased vocal tract impedance lead to increased vocal efficiency and economy. One of the major changes was the more prominent singer's/speaker's formant cluster. Vocal tract and glottal modifications were more prominent during and after straw exercising compared with tube phonation.

Case-study magnetic resonance imaging and acoustic investigation of the effects of vocal warm-up on two voice professionals.

Vocal warm-up (WU)-related changes were studied in one male musical singer and one female speech trainer. They sustained vowels before and after WU in a magnetic resonance imaging (MRI) device. Acoustic recordings were made in a studio. The vocal tract area increased after WU, a formant cluster appeared between 2 and 4.5 kHz, and SPL increased. Evidence of larynx lowering was only found for the male. The pharyngeal inlet over the epilaryngeal outlet ratio (A(ph)/A(e)) increased by 10%-28%, being 3-4 for the male and 5-7 for the female. The results seem to represent different voice training traditions. A singer's formant cluster may be achievable without a high A(ph)/A(e) (≥ 6), but limitations of the 2D method should be taken into account.

Bifurcations and chaos in register transitions of excised larynx experiments.

Experimental data from an excised larynx are analyzed in the light of nonlinear dynamics. The excised larynx provides an experimental framework that enables artificial control and direct observation of the vocal fold vibrations. Of particular interest in this experiment is the coexistence of two distinct vibration patterns, which closely resemble chest and falsetto registers of the human voice. Abrupt transitions between the two registers are typically accompanied by irregular vibrations. Two approaches are presented for the modeling of the excised larynx experiment; one is the nonlinear predictive modeling of the experimental time series and the other is the biomechanical modeling (three-mass model) that takes into account basic mechanisms of the vocal fold vibrations. The two approaches show that the chest and falsetto vibrations correspond to two coexisting limit cycles, which jump to each other with a change in the bifurcation parameter. Irregular vibrations observed at the register jumps are due to chaos that exists near the two limit cycles. This provides an alternative mechanism to generate chaotic vibrations in excised larynx experiment, which is different from the conventionally known mechanisms such as strong asymmetry between the left and right vocal folds or excessively high subglottal pressure.