Dynamic vocal fold parameters with changing adduction in ex-vivo hemilarynx experiments

Three-Dimensional Analysis of Vocal Fold Oscillations: Correlating Superior and Medial Surface Dynamics Using Ex Vivo Human Hemilarynges

The primary acoustic signal of the voice is generated by the complex oscillation of the vocal folds (VFs), whereby physicians can barely examine the medial VF surface due to its anatomical inaccessibility. In this study, we investigated possibilities to infer medial surface dynamics by analyzing correlations in the oscillatory behavior of the superior and medial VF surfaces of four human hemilarynges, each in 24 different combinations of flow rate, VF adduction, and elongation. The two surfaces were recorded synchronously during sustained phonation using two high-speed camera setups and were subsequently 3D-reconstructed. The 3D surface parameters of mean and maximum velocities and displacements and general phonation parameters were calculated. The VF oscillations were also analyzed using empirical eigenfunctions (EEFs) and mucosal wave propagation, calculated from medial surface trajectories. Strong linear correlations were found between the 3D parameters of the superior and medial VF surfaces, ranging from 0.8 to 0.95. The linear regressions showed similar values for the maximum velocities at all hemilarynges (0.69-0.9), indicating the most promising parameter for predicting the medial surface. Since excessive VF velocities are suspected to cause phono-trauma and VF polyps, this parameter could provide added value to laryngeal diagnostics in the future.

Asymmetric triangular body-cover model of the vocal folds with bilateral intrinsic muscle activation

Many voice disorders are linked to imbalanced muscle activity and known to exhibit asymmetric vocal fold vibration. However, the relation between imbalanced muscle activation and asymmetric vocal fold vibration is not well understood. This study introduces an asymmetric triangular body-cover model of the vocal folds, controlled by the activation of bilateral intrinsic laryngeal muscles, to investigate the effects of muscle imbalance on vocal fold oscillation. Various scenarios were considered, encompassing imbalance in individual muscles and muscle pairs, as well as accounting for asymmetry in lumped element parameters. Measurements of amplitude and phase asymmetries were employed to match the oscillatory behavior of two pathological cases: unilateral paralysis and muscle tension dysphonia. The resulting simulations exhibit muscle imbalance consistent with expectations in the composition of these voice disorders, yielding asymmetries exceeding 30% for paralysis and below 5% for dysphonia. This underscores the relevance of muscle imbalance in representing phonatory scenarios and its potential for characterizing asymmetry in vocal fold vibration.

Reconstruction of Vocal Fold Medial Surface 3D Trajectories: Effects of Neuromuscular Stimulation and Airflow

INTRODUCTION

Analysis of medial surface dynamics of the vocal folds (VF) is critical to understanding voice production and treatment of voice disorders. We analyzed VF medial surface vibratory dynamics, evaluating the effects of airflow and nerve stimulation using 3D reconstruction and empirical eigenfunctions (EEF).

STUDY DESIGN

In vivo canine hemilarynx phonation.

METHODS

An in vivo canine hemilarynx was phonated while graded stimulation of the recurrent and superior laryngeal nerves (RLN and SLN) was performed. For each phonatory condition, vibratory cycles were 3D reconstructed from tattooed landmarks on the VF medial surface at low, medium, and high airflows. Parameters describing medial surface trajectory shape were calculated, and underlying patterns were emphasized using EEFs. Fundamental frequency and smoothed cepstral peak prominence (CPPS) were calculated from acoustic data.

RESULTS

Convex-hull area of landmark trajectories increased with increasing flow and decreasing nerve activation level. Trajectory shapes observed included circular, ellipsoid, bent, and figure-eight. They were more circular on the superior and anterior VF, and more elliptical and line-like on the inferior and posterior VF. The EEFs capturing synchronal opening and closing (EEF1) and alternating convergent/divergent (EEF2) glottis shapes were mostly unaffected by flow and nerve stimulation levels. CPPS increased with higher airflow except for low RLN activation and very dominant SLN stimulation.

CONCLUSION

We analyzed VF vibration as a function of neuromuscular stimulation and airflow levels. Oscillation patterns such as figure-eight and bent trajectories were linked to high nerve activation and flow. Further studies investigating longer sections of 3D reconstructed oscillations are needed.

LEVEL OF EVIDENCE

N/A, Basic Science Laryngoscope, 134:1249-1257, 2024.

Effects of Thyroarytenoid Activation Induced Vibratory Asymmetry on Voice Acoustics and Perception

INTRODUCTION

Asymmetry of vocal fold (VF) vibration is common in patients with voice complaints and also observed in 10% of normophonic individuals. Although thyroarytenoid (TA) muscle activation plays a crucial role in regulating VF vibration, how TA activation asymmetry relates to voice acoustics and perception is unclear. We evaluated the relationship between TA activation asymmetry and the resulting acoustics and perception.

METHODS

An in vivo canine model of phonation was used to create symmetric and increasingly asymmetric VF vibratory conditions via graded stimulation of bilateral TA muscles. Naïve listeners (n = 89) rated the perceptual quality of 100 unique voice samples using a visual sort-and-rate task. For each phonatory condition, cepstral peak prominence (CPP), harmonic amplitude (H1-H2), and root-mean-square (RMS) energy of the voice were measured. The relationships between these metrics, vibratory asymmetry, and perceptual ratings were evaluated.

RESULTS

Increasing levels of TA asymmetry resulted in declining listener preference. Furthermore, only severely asymmetric audio samples were perceptually distinguishable from symmetric and mildly asymmetric conditions. CPP was negatively correlated with TA asymmetry: voices produced with larger degrees of asymmetry were associated with lower CPP values. Listeners preferred audio samples with higher values of CPP, high RMS energy, and lower H1-H2 (less breathy).

CONCLUSION

Listeners are sensitive to changes in voice acoustics related to vibratory asymmetry. Although increasing vibratory asymmetry is correlated with decreased perceptual ratings, mild asymmetries are perceptually tolerated. This study contributes to our understanding of voice production and quality by identifying perceptually salient and clinically meaningful asymmetry.

LEVEL OF EVIDENCE

N/A (Basic Science Study) Laryngoscope, 134:1327-1332, 2024.

What kind of phonation causes the strongest vocal fold collision? - A hemi-larynx phonation contact pressure study

This paper presents a measurement setup which is able to measure the distribution of small scale pressure on an area of 15.2 mm × 30.4 mm with a sample rate up to 1.2 kHz. It was used to investigate the contact pressures of vocal folds during phonation. This was performed in ex vivo experiments of 11 porcine larynges. The contact pressure at the medial surface and other phonation parameters, as the glottal resistance and the closing velocity of the vocal fold, were measured at different adduction and elongation levels and air flow rates. A statistical analysis was carried out. It could be shown that the contact pressure rises, when the vocal fold is manipulated or when the flow rate is increased.

Neural network-based estimation of biomechanical vocal fold parameters

Vocal fold (VF) vibrations are the primary source of human phonation. High-speed video (HSV) endoscopy enables the computation of descriptive VF parameters for assessment of physiological properties of laryngeal dynamics, i.e., the vibration of the VFs. However, underlying biomechanical factors responsible for physiological and disordered VF vibrations cannot be accessed. In contrast, physically based numerical VF models reveal insights into the organ's oscillations, which remain inaccessible through endoscopy. To estimate biomechanical properties, previous research has fitted subglottal pressure-driven mass-spring-damper systems, as inverse problem to the HSV-recorded VF trajectories, by global optimization of the numerical model. A neural network trained on the numerical model may be used as a substitute for computationally expensive optimization, yielding a fast evaluating surrogate of the biomechanical inverse problem. This paper proposes a convolutional recurrent neural network (CRNN)-based architecture trained on regression of a physiological-based biomechanical six-mass model (6 MM). To compare with previous research, the underlying biomechanical factor "subglottal pressure" prediction was tested against 288 HSV ex vivo porcine recordings. The contributions of this work are two-fold: first, the presented CRNN with the 6 MM handles multiple trajectories along the VFs, which allows for investigations on local changes in VF characteristics. Second, the network was trained to reproduce further important biomechanical model parameters like VF mass and stiffness on synthetic data. Unlike in a previous work, the network in this study is therefore an entire surrogate of the inverse problem, which allowed for explicit computation of the fitted model using our approach. The presented approach achieves a best-case mean absolute error (MAE) of 133 Pa (13.9%) in subglottal pressure prediction with 76.6% correlation on experimental data and a re-estimated fundamental frequency MAE of 15.9 Hz (9.9%). In-detail training analysis revealed subglottal pressure as the most learnable parameter. With the physiological-based model design and advances in fast parameter prediction, this work is a next step in biomechanical VF model fitting and the estimation of laryngeal kinematics.

Flow-induced oscillations of vocal-fold replicas with tuned extensibility and material properties

Human vocal folds are highly deformable non-linear oscillators. During phonation, they stretch up to 50% under the complex action of laryngeal muscles. Exploring the fluid/structure/acoustic interactions on a human-scale replica to study the role of the laryngeal muscles remains a challenge. For that purpose, we designed a novel in vitro testbed to control vocal-folds pre-phonatory deformation. The testbed was used to study the vibration and the sound production of vocal-fold replicas made of (i) silicone elastomers commonly used in voice research and (ii) a gelatin-based hydrogel we recently optimized to approximate the mechanics of vocal folds during finite strains under tension, compression and shear loadings. The geometrical and mechanical parameters measured during the experiments emphasized the effect of the vocal-fold material and pre-stretch on the vibration patterns and sounds. In particular, increasing the material stiffness increases glottal flow resistance, subglottal pressure required to sustain oscillations and vibratory fundamental frequency. In addition, although the hydrogel vocal folds only oscillate at low frequencies (close to 60 Hz), the subglottal pressure they require for that purpose is realistic (within the range 0.5-2 kPa), as well as their glottal opening and contact during a vibration cycle. The results also evidence the effect of adhesion forces on vibration and sound production.

Effect of Ligament Fibers on Dynamics of Synthetic, Self-Oscillating Vocal Folds in a Biomimetic Larynx Model

Synthetic silicone larynx models are essential for understanding the biomechanics of physiological and pathological vocal fold vibrations. The aim of this study is to investigate the effects of artificial ligament fibers on vocal fold vibrations in a synthetic larynx model, which is capable of replicating physiological laryngeal functions such as elongation, abduction, and adduction. A multi-layer silicone model with different mechanical properties for the musculus vocalis and the lamina propria consisting of ligament and mucosa was used. Ligament fibers of various diameters and break resistances were cast into the vocal folds and tested at different tension levels. An electromechanical setup was developed to mimic laryngeal physiology. The measurements included high-speed video recordings of vocal fold vibrations, subglottal pressure and acoustic. For the evaluation of the vibration characteristics, all measured values were evaluated and compared with parameters from ex and in vivo studies. The fundamental frequency of the synthetic larynx model was found to be approximately 200-520 Hz depending on integrated fiber types and tension levels. This range of the fundamental frequency corresponds to the reproduction of a female normal and singing voice range. The investigated voice parameters from vocal fold vibration, acoustics, and subglottal pressure were within normal value ranges from ex and in vivo studies. The integration of ligament fibers leads to an increase in the fundamental frequency with increasing airflow, while the tensioning of the ligament fibers remains constant. In addition, a tension increase in the fibers also generates a rise in the fundamental frequency delivering the physiological expectation of the dynamic behavior of vocal folds.

Predicting 3D soft tissue dynamics from 2D imaging using physics informed neural networks

Tissue dynamics play critical roles in many physiological functions and provide important metrics for clinical diagnosis. Capturing real-time high-resolution 3D images of tissue dynamics, however, remains a challenge. This study presents a hybrid physics-informed neural network algorithm that infers 3D flow-induced tissue dynamics and other physical quantities from sparse 2D images. The algorithm combines a recurrent neural network model of soft tissue with a differentiable fluid solver, leveraging prior knowledge in solid mechanics to project the governing equation on a discrete eigen space. The algorithm uses a Long-short-term memory-based recurrent encoder-decoder connected with a fully connected neural network to capture the temporal dependence of flow-structure-interaction. The effectiveness and merit of the proposed algorithm is demonstrated on synthetic data from a canine vocal fold model and experimental data from excised pigeon syringes. The results showed that the algorithm accurately reconstructs 3D vocal dynamics, aerodynamics, and acoustics from sparse 2D vibration profiles.

Collision Pressure and Dissipated Power Dose in a Self-Oscillating Silicone Vocal Fold Model With a Posterior Glottal Opening

PURPOSE

The goal of this study was to experimentally evaluate how compensating for the adverse acoustic effects of a posterior glottal opening (PGO) by increasing subglottal pressure and changing supraglottal compression, as have been associated with vocal hyperfunction, influences the risk of vocal fold (VF) trauma.

METHOD

A self-oscillating synthetic silicone model of the VFs with an airflow bypass that modeled a PGO was investigated in a hemilaryngeal flow facility. The influence of compensatory mechanisms on collision pressure and dissipated collision power was investigated for different PGO areas and supraglottal compression. Compensatory behaviors were mimicked by increasing the subglottal pressure to achieve a target sound pressure level (SPL).

RESULTS

Increasing the subglottal pressure to compensate for decreased SPL due to a PGO produced higher values for both collision pressure and dissipated collision power. Whereas a 10-mm² PGO area produced a 12% increase in the peak collision pressure, the dissipated collision power increased by 122%, mainly due to an increase in the magnitude of the collision velocity. This suggests that the value of peak collision pressure may not fully capture the mechanisms by which phonotrauma occurs. It was also found that an optimal value of supraglottal compression exists that maximizes the radiated SPL, indicating the potential utility of supraglottal compression as a compensatory mechanism.

CONCLUSIONS

Larger PGO areas are expected to increase the risk of phonotrauma due to the concomitant increase in dissipated collision power associated with maintaining SPL. Furthermore, the risk of VF damage may not be fully characterized by only the peak collision pressure.

Simulated Laryngeal High-Speed Videos for the Study of Normal and Dysphonic Vocal Fold Vibration

PURPOSE

Laryngeal high-speed videoendoscopy (LHSV) has been recognized as a highly valuable modality for the scientific investigations of vocal fold (VF) vibrations. In contrast to stroboscopic imaging, LHSV enables visualizing aperiodic VF vibrations. However, the technique is less well established in the clinical care of disordered voices, partly because the properties of aperiodic vibration patterns are not yet described comprehensively. To address this, a computer model for simulation of VF vibration patterns observed in a variety of different phonation types is proposed.

METHOD

A previously published kinematic model of mucosal wave phenomena is generalized to be capable of left-right asymmetry and to simulate endoscopic videos instead of only kymograms of VF vibrations at single sagittal positions. The most influential control parameters are the glottal halfwidths, the oscillation frequencies, the amplitudes, and the phase delays.

RESULTS

The presented videos demonstrate zipper-like vibration, pressed voice, voice onset, constant and time-varying left-right and anterior-posterior phase differences, as well as left-right frequency differences of the VF vibration. Video frames, videokymograms, phonovibrograms, glottal area waveforms, and waveforms of VF contact area relating to electroglottograms are shown, as well as selected kinematic parameters.

CONCLUSION

The presented videos demonstrate the ability to produce vibration patterns that are similar to those typically seen in endoscopic videos obtained from vocally healthy and dysphonic speakers.

SUPPLEMENTAL MATERIAL

https://doi.org/10.23641/asha.20151833.

Computational Modeling of Voice Production Using Excised Canine Larynx

A combined experimental-numerical work was conducted to comprehensively validate a subject-specific continuum model of voice production in larynx using excised canine laryngeal experiments. The computational model is a coupling of the Navier-Stokes equations for glottal flow dynamics and a finite element model of vocal fold dynamics. The numerical simulations employed a cover-body vocal fold structure with the geometry reconstructed from magnetic resonance imaging scans and the material properties determined through an optimization-based inverse process of experimental indentation measurement. The results showed that the simulations predicted key features of the dynamics observed in the experiments, including the skewing of the glottal flow waveform, mucosal wave propagation, continuous increase of the divergent angle and intraglottal swirl strength during glottal closing, and flow recirculation between glottal jet and vocal fold. The simulations also predicted the increase of the divergent angle, glottal jet speed, and intraglottal flow swirl strength with the subglottal pressure, same as in the experiments. Quantitatively, the simulations over-predicted the frequency and jet speed and under-predicted the flow rate and divergent angle for the larynx under study. The limitations of the model and their implications were discussed.

3D-FV-FE Aeroacoustic Larynx Model for Investigation of Functional Based Voice Disorders

For the clinical analysis of underlying mechanisms of voice disorders, we developed a numerical aeroacoustic larynx model, called simVoice, that mimics commonly observed functional laryngeal disorders as glottal insufficiency and vibrational left-right asymmetries. The model is a combination of the Finite Volume (FV) CFD solver Star-CCM+ and the Finite Element (FE) aeroacoustic solver CFS++. simVoice models turbulence using Large Eddy Simulations (LES) and the acoustic wave propagation with the perturbed convective wave equation (PCWE). Its geometry corresponds to a simplified larynx and a vocal tract model representing the vowel /a/. The oscillations of the vocal folds are externally driven. In total, 10 configurations with different degrees of functional-based disorders were simulated and analyzed. The energy transfer between the glottal airflow and the vocal folds decreases with an increasing glottal insufficiency and potentially reflects the higher effort during speech for patients being concerned. This loss of energy transfer may also have an essential influence on the quality of the sound signal as expressed by decreasing sound pressure level (SPL), Cepstral Peak Prominence (CPP), and Vocal Efficiency (VE). Asymmetry in the vocal fold oscillations also reduces the quality of the sound signal. However, simVoice confirmed previous clinical and experimental observations that a high level of glottal insufficiency worsens the acoustic signal quality more than oscillatory left-right asymmetry. Both symptoms in combination will further reduce the quality of the sound signal. In summary, simVoice allows for detailed analysis of the origins of disordered voice production and hence fosters the further understanding of laryngeal physiology, including occurring dependencies. A current walltime of 10 h/cycle is, with a prospective increase in computing power, auspicious for a future clinical use of simVoice.

Numerical and experimental investigations on vocal fold approximation in healthy and simulated unilateral vocal fold paralysis

We have developed a novel surgical/computational model for the investigation of unilateral vocal fold paralysis (UVFP) which will be used to inform future in silico approaches to improve surgical outcomes in type I thyroplasty. Healthy phonation (HP) was achieved using cricothyroid suture approximation on both sides of the larynx to generate symmetrical vocal fold closure. Following high-speed videoendoscopy (HSV) capture, sutures on the right side of the larynx were removed, partially releasing tension unilaterally and generating asymmetric vocal fold closure characteristic of UVFP (sUVFP condition). HSV revealed symmetric vibration in HP, while in sUVFP the sutured side demonstrated a higher frequency (10 - 11%). For the computational model, ex vivo magnetic resonance imaging (MRI) scans were captured at three configurations: non-approximated (NA), HP, and sUVFP. A finite-element method (FEM) model was built, in which cartilage displacements from the MRI images were used to prescribe the adduction and the vocal fold deformation was simulated before the eigenmode calculation. The results showed that the frequency comparison between the two sides were consistent with observations from HSV. This alignment between the surgical and computational models supports the future application of these methods for the investigation of treatment for UVFP.

Analysis of the tonal sound generation during phonation with and without glottis closure

The human phonation is characterized by periodical oscillations of the vocal folds with a complete glottis closure. In contrast, a glottal insufficiency (GI) represents an oscillation without glottis closure resulting in a breathy and weak voice. In this study, flow-induced oscillations of silicone vocal folds were modeled with and without glottis closure. The measurements comprised the flow pressure in the model, the generated sound, and the high-speed footage of the vocal fold motion. The analysis revealed that the sound signal for vocal fold oscillations without closure exhibits a lower number of harmonic tones with smaller amplitudes compared to the case with complete closure. The time series of the pressure signals showed small and periodical oscillations occurring less frequently and with smaller amplitude for the GI case. Accordingly, the pressure spectra include fewer harmonics similar to the sound. The analysis of the high-speed videos indicates that the strength of the pressure oscillations correlates with the divergence angle of the glottal duct during the closing motion. Physiologically, large divergence angles typically occur for a pronounced mucosal wave motion with glottis closure. Thus, the results indicate a correlation between the intensity of the mucosal wave and the development of harmonic tones.

Toward Development of a Vocal Fold Contact Pressure Probe: Bench-Top Validation of a Dual-Sensor Probe Using Excised Human Larynx Models

A critical element in understanding voice production mechanisms is the characterization of vocal fold collision, which is widely considered a primary etiological factor in the development of common phonotraumatic lesions such as nodules and polyps. This paper describes the development of a transoral, dual-sensor intraglottal/subglottal pressure probe for the simultaneous measurement of vocal fold collision and subglottal pressures during phonation using two miniature sensors positioned 7.6 mm apart at the distal end of a rigid cannula. Proof-of-concept testing was performed using excised whole-mount and hemilarynx human tissue aerodynamically driven into self-sustained oscillation, with systematic variation of the superior-inferior positioning of the vocal fold collision sensor. In the hemilarynx experiment, signals from the pressure sensors were synchronized with an acoustic microphone, a tracheal-surface accelerometer, and two high-speed video cameras recording at 4000 frames per second for top-down and en face imaging of the superior and medial vocal fold surfaces, respectively. As expected, the intraglottal pressure signal exhibited an impulse-like peak when vocal fold contact occurred, followed by a broader peak associated with intraglottal pressure build-up during the de-contacting phase. As subglottal pressure was increased, the peak amplitude of the collision pressure increased and typically reached a value below that of the average subglottal pressure. Results provide important baseline vocal fold collision pressure data with which computational models of voice production can be developed and in vivo measurements can be referenced.

Juvenile Ovine Ex Vivo Larynges: Phonatory, Histologic, and Micro CT Based Anatomic Analyses

It is well known that the phonatory process changes during the life span. However, detailed investigations on potential factors concerned are rare. To deal with this issue, we performed extended biomechanical, macro anatomical, and histological analyses of the contributing laryngeal structures in ex vivo juvenile sheep models. Altogether twelve juvenile sheep larynges were analyzed within the phonatory experiments. Three different elongation levels and 16 different flow levels were applied to achieve a large variety of phonatory conditions. Vocal fold dynamics and acoustical and subglottal signals could be analyzed for 431 experimental runs. Subsequently, for six juvenile larynges microcomputed tomography following virtual 3D reconstruction was performed. The remaining six juvenile larynges as well as six ex vivo larynges from old sheep were histologically and immunohistologically analyzed. Results for juveniles showed more consistent dynamical behavior compared to old sheep larynges due to vocal fold tissue alterations during the life span. The phonatory process in juvenile sheep seems to be more effective going along with a greater dynamic range. These findings are supported by the histologically detected higher amounts of elastin and hyaluronic acid in the lamina propria of the juvenile sheep. The 3D reconstructions of the thyro-arytenoid muscles (TAM) showed a symmetrical shape. Intraindividual volume and surface differences of the TAM were small and comparable to those of aged sheep. However, TAM dimensions were statistically significant smaller for juvenile larynges. Finally, topographical landmarks were introduced for later comparison with other individuals and species. This work resulted in detailed functional, immunohistological, and anatomical information that was not yet reported. This data will also provide reference information for therapeutic strategies regarding aging effects, e.g. laryngeal muscle treatment by functional electrical stimulation.

A computational study of depth of vibration into vocal fold tissues

The effective depth of vocal fold vibration is self-regulated and generally not known a priori in vocalization. In this study, the effective depth was quantified systematically under various phonatory conditions using a fiber-gel finite element vocal fold model. The horizontal and vertical excursions of each finite element nodal point trajectory were recorded to compute trajectory areas. The extent of vibration was then studied based on the variation of trajectory radii as a function of depth in several coronal sections along the anterior-posterior direction. The results suggested that the vocal fold nodal trajectory excursions decrease systematically as a function of depth but are affected by the layered structure of the vocal folds. The effective depth of vibration was found to range between 15 and 55% of the total anatomical depth across all phonatory conditions. The nodal trajectories from the current study were compared qualitatively with the results from excised human hemi-larynx experiments published in Döllinger and Berry [(2006). J. Voice. 20(3), 401-413]. An estimate of the effective mass of a one-mass vocal fold model was also computed based on the effective depth of vibration observed in this study under various phonatory conditions.

Visual and Automatic Evaluation of Vocal Fold Mucosal Waves Through Sharpness of Lateral Peaks in High-Speed Videokymographic Images

INTRODUCTION

The sharpness of lateral peaks is a visually helpful clinical feature in high-speed videokymographic (VKG) images indicating vertical phase differences and mucosal waves on the vibrating vocal folds and giving insights into the health and pliability of vocal fold mucosa. This study aims at investigating parameters that can be helpful in objectively quantifying the lateral peak sharpness from the VKG images.

METHOD

Forty-five clinical VKG images with different degrees of sharpness of lateral peaks were independently evaluated visually by three raters. The ratings were compared to parameters obtained by automatic image analysis of the vocal fold contours: Open Time Percentage Quotients (OTQ) and Plateau Quotients (PQ). The OTQ parameters were derived as fractions of the period during which the vocal fold displacement exceeds a predetermined percentage of the vibratory amplitude. The PQ parameters were derived similarly but as a fraction of the open phase instead of a period.

RESULTS

The best correspondence between the visual ratings and the automatically derived quotients were found for the OTQ and PQ parameters derived at 95% and 80% of the amplitude, named OTQ₉₅, PQ₉₅, OTQ₈₀ and PQ₈₀. Their Spearman's rank correlation coefficients were in the range of 0.73 to 0.77 (P < 0.001) indicating strong relationships with the visual ratings. The strengths of these correlations were similar to those found from inter-rater comparisons of visual evaluations of peak sharpness.

CONCLUSION

The Open time percentage and Plateau quotients at 95% and 80% of the amplitude stood out as the possible candidates for capturing the sharpness of the lateral peaks with their reliability comparable to that of visual ratings.

Development of Excised Larynx

The larynx is a complex organ which has a role in a variety of functions such as phonation, breathing, and swallowing. To research these functions, it is widely accepted that in vivo studies provide more anatomically and physiologically relevant findings. However, invasive procedures are generally needed to measure variables such a subglottal pressure, vocal fold tension and stiffness, and cricothyroid muscle stretch. Performing studies using excised larynges is a useful technique which makes it possible to not only measure phonation parameters but control them as well. Early studies using excised larynges mainly focused on controlling specific parameters and mathematical modeling simulations. The use of these studies has helped further research in laryngeal anatomy, imaging techniques, as well as aerodynamic, acoustic, and biomechanical properties. Here, we describe the progress of this research over the past 5 years. The number of accepted animal models has increased and ideas from excised larynx studies are starting to be applied to treatment methods for laryngeal disorders. These experiments are only valid for an excised situation and must continue to be combined with animal experimentation and clinical observations.

Investigation of phonatory characteristics using ex vivo rabbit larynges

Quantitative analysis of phonatory characteristics of rabbits has been widely neglected. However, preliminary studies established the rabbit larynx as a potential model of human phonation. This study reports quantitative data on phonation using ex vivo rabbit larynx models to achieve more insight into dependencies of three main components of the phonation process, including airflow, vocal fold dynamics, and the acoustic output. Sustained phonation was induced in 11 ex vivo rabbit larynges. For 414 phonatory conditions, vocal fold vibrations, acoustic, and aerodynamic parameters were analyzed as functions of longitudinal vocal fold pre-stress, applied air flow, and glottal closure insufficiency. Dimensions of the vocal folds were measured and histological data were analyzed. Glottal closure characteristics improved for increasing longitudinal pre-stress and applied airflow. For the subglottal pressure signal only the cepstral peak prominence showed dependency on glottal closure. In contrast, vibrational, acoustic, and aerodynamic parameters were found to be highly dependent on the degree of glottal closure: The more complete the glottal closure during phonation, the better the aerodynamic and acoustic characteristics. Hence, complete or at least partial glottal closure appears to enhance acoustic signal quality. Finally, results validate the ex vivo rabbit larynx as an effective model for analyzing the phonatory process.

Hemi-laryngeal Setup for Studying Vocal Fold Vibration in Three Dimensions

The voice of humans and most non-human mammals is generated in the larynx through self-sustaining oscillation of the vocal folds. Direct visual documentation of vocal fold vibration is challenging, particularly in non-human mammals. As an alternative, excised larynx experiments provide the opportunity to investigate vocal fold vibration under controlled physiological and physical conditions. However, the use of a full larynx merely provides a top view of the vocal folds, excluding crucial portions of the oscillating structures from observation during their interaction with aerodynamic forces. This limitation can be overcome by utilizing a hemi-larynx setup where one half of the larynx is mid-sagittally removed, providing both a superior and a lateral view of the remaining vocal fold during self-sustained oscillation. Here, a step-by-step guide for the anatomical preparation of hemi-laryngeal structures and their mounting on the laboratory bench is given. Exemplary phonation of the hemi-larynx preparation is documented with high-speed video data captured by two synchronized cameras (superior and lateral views), showing three-dimensional vocal fold motion and corresponding time-varying contact area. The documentation of the hemi-larynx setup in this publication will facilitate application and reliable repeatability in experimental research, providing voice scientists with the potential to better understand the biomechanics of voice production.

Biomechanical simulation of vocal fold dynamics in adults based on laryngeal high-speed videoendoscopy

MOTIVATION

Human voice is generated in the larynx by the two oscillating vocal folds. Owing to the limited space and accessibility of the larynx, endoscopic investigation of the actual phonatory process in detail is challenging. Hence the biomechanics of the human phonatory process are still not yet fully understood. Therefore, we adapt a mathematical model of the vocal folds towards vocal fold oscillations to quantify gender and age related differences expressed by computed biomechanical model parameters.

METHODS

The vocal fold dynamics are visualized by laryngeal high-speed videoendoscopy (4000 fps). A total of 33 healthy young subjects (16 females, 17 males) and 11 elderly subjects (5 females, 6 males) were recorded. A numerical two-mass model is adapted to the recorded vocal fold oscillations by varying model masses, stiffness and subglottal pressure. For adapting the model towards the recorded vocal fold dynamics, three different optimization algorithms (Nelder-Mead, Particle Swarm Optimization and Simulated Bee Colony) in combination with three cost functions were considered for applicability. Gender differences and age-related kinematic differences reflected by the model parameters were analyzed.

RESULTS AND CONCLUSION

The biomechanical model in combination with numerical optimization techniques allowed phonatory behavior to be simulated and laryngeal parameters involved to be quantified. All three optimization algorithms showed promising results. However, only one cost function seems to be suitable for this optimization task. The gained model parameters reflect the phonatory biomechanics for men and women well and show quantitative age- and gender-specific differences. The model parameters for younger females and males showed lower subglottal pressures, lower stiffness and higher masses than the corresponding elderly groups. Females exhibited higher subglottal pressures, smaller oscillation masses and larger stiffness than the corresponding similar aged male groups. Optimizing numerical models towards vocal fold oscillations is useful to identify underlying laryngeal components controlling the phonatory process.

Influence of glottal closure on the phonatory process in ex vivo porcine larynges

Many cases of disturbed voice signals can be attributed to incomplete glottal closure, vocal fold oscillation asymmetries, and aperiodicity. Often these phenomena occur simultaneously and interact with each other, making a systematic, isolated investigation challenging. Therefore, ex vivo porcine experiments were performed which enable direct control of glottal configurations. Different pre-phonatory glottal gap sizes, adduction levels, and flow rates were adjusted. The resulting glottal closure types were identified in a post-processing step. Finally, the acoustic quality, aerodynamic parameters, and the characteristics of vocal fold oscillation were analyzed in reference to the glottal closure types. Results show that complete glottal closure stabilizes the phonation process indicated through a reduced left-right phase asymmetry, increased amplitude and time periodicity, and an increase in the acoustic quality. Although asymmetry and periodicity parameter variation covers only a small range of absolute values, these small variations have a remarkable influence on the acoustic quality. Due to the fact that these parameters cannot be influenced directly, the authors suggest that the (surgical) reduction of the glottal gap seems to be a promising method to stabilize the phonatory process, which has to be confirmed in future studies.

Phonation Analysis Combined with 3D Reconstruction of the Thyroarytenoid Muscle in Aged Ovine Ex Vivo Larynx Models

OBJECTIVE

The aim of the study was to establish a basic data set of combined functional and anatomical measures of aged sheep larynges using ex vivo models. Combining these two approaches in one and the same larynx is an unmet goal so far yet is important as newer treatment strategies aim to preserve the organ structure and new assessment tools are required. Ovine larynges were used as their dimensions, and muscle fiber type distribution highly resemble the human larynx.

STUDY DESIGN

Ex vivo animal study.

METHODS

Larynges of six sheep (~9 years of age) were subjected to ex vivo functional phonatory experiments. Phonatory characteristics were analyzed as a function of longitudinal vocal fold (VF) prestress. Anatomical measurements of the same larynges comprised micro-computed tomography scans followed by three-dimensional (3D) reconstructions. Using specially adapted radiological scan protocols with subsequent 3D reconstruction, muscle volumes, surface areas, and anatomical measurements were computed.

RESULTS

Increasing longitudinal prestress yielded higher subglottal pressure (P_S) for the same airflow. Quantitative differences to previous studies-such as the increased P_S and increased phonation threshold pressure-were detected. We achieved excellent visualization of the laryngeal muscles and framework, resulting in accurate 3D reconstructions for quantitative analysis. We found no significant intraindividual volume differences of the thyroarytenoid muscles.

CONCLUSION

The established protocol allows precise functional and anatomical measures. The data created provide a reference data set for upcoming therapeutic strategies (eg, growth factor therapy, functional electrical stimulation) that target essential structures of the VFs such as the laryngeal muscles and/or the VF mucosa.

What Can be Learned from the Time Course of Changes in Low-Frequency Stimulated Muscle?

Not available.

Automated setup for ex vivo larynx experiments

Ex vivo larynx experiments are limited in time due to degeneration of the laryngeal tissues. In order to acquire a significant and comparable amount of data, automatization of current manual experimental procedures is desirable. A computer controlled, electro-mechanical setup was developed for time-dependent variation of specific physiological parameters, including adduction and elongation level of the vocal folds and glottal flow. The setup offers a standardized method to induce defined forces on the laryngeal cartilages. Furthermore, phonation onset is detected automatically and the subsequent measurement procedure is automated and standardized to improve the efficiency of the experimental process. The setup was validated using four ex vivo porcine larynges, whereas each validation measurement series was executed with one separate larynx. Altogether 31 single measurements were undertaken, which can be summed up to a total experimental time of about 4 min. Vocal fold elongation and adduction lead both to an increase in fundamental frequency and subglottal pressure. Measurement procedures like applying defined subglottal pressure steps and onset-offset detection were reliably executed. The setup allows for a computer-based parameter control, which enables fast experimental execution over a wide range of laryngeal configurations. This maximizes the number of measurements and reduces personal effort compared with manual procedures.