Modeling viscous dissipation during vocal fold contact: the influence of tissue viscosity and thickness with implications for hydration

Estimating the Pathophysiology of Phonotraumatic Vocal Hyperfunction Using Ambulatory Data and a Computational Model

PURPOSE

This study uses a voice production model to estimate muscle activation levels and subglottal pressure (PS) in patients with phonotraumatic vocal hyperfunction (PVH), based on ambulatory measurements of sound pressure level (SPL) and spectral tilt (H1-H2). In addition, variations in these physiological parameters are evaluated with respect to different values of the Daily Phonotrauma Index (DPI).

METHOD

The study obtained ambulatory voice data from patients diagnosed with PVH and a matched control group. To infer physiological parameters, ambulatory data were mapped onto synthetic data generated by a physiologically relevant voice production model. Inverse mapping strategies involved selecting model simulations that represented ambulatory distributions using stochastic (random) sampling weighted by probability with which different vowels occur in English. A categorical approach assessed the relationship between different values of DPI and changes in estimated physiological parameters.

RESULTS

Results showed significant differences between the PVH and control groups in key parameters, including statistical moments of H1-H2, SPL, PS, and muscle activity of lateral cricoarytenoid (LCA) and cricothyroid (CT) muscles. Higher DPI values, reflecting more severe PVH, were associated with increased mean LCA activation and decreased LCA variability, along with decreased mean CT activation and increased median PS. These findings highlight the relationship between muscle activation patterns, PS, and the severity of vocal pathology as indicated by the DPI. It is hypothesized that a major driver of muscle activation and PS changes is the variation in maladaptive adjustments (vocal effort) when compensating for the presence of vocal pathology.

CONCLUSIONS

This study demonstrated that noninvasive ambulatory voice data could be used to drive a voice production modeling process, providing valuable insights into underlying physiological parameters associated with PVH. Future research will focus on refining the predictive power of the modeling process and exploring the implications of these findings in further delineating the etiology and pathophysiology of PVH, with the ultimate goal to develop improved methods for the prevention, diagnosis, and treatment of PVH.

SUPPLEMENTAL MATERIAL

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

Biophysical aspects of mechanotransduction in cells and their physiological/biological implications in vocal fold vibration: a narrative review

Mechanotransduction is a crucial property in all organisms, modulating cellular behaviors in response to external mechanical stimuli. Given the high mobility of vocal folds, it is hypothesized that mechanotransduction significantly contributes to their tissue homeostasis. Recent studies have identified mechanosensitive proteins in vocal fold epithelia, supporting this hypothesis. Voice therapy, which, involves the mobilization of vocal folds, aims to rehabilitate vocal function and restore homeostasis. However, establishing a direct causal link between specific mechanical stimuli and therapeutic benefits is challenging due to the variability in voice therapy techniques. This challenge is further compounded when investigating biological benefits in humans. Vocal fold tissue cannot be biopsied without significant impairment of the vibratory characteristics of the vocal folds. Conversely, studies using vocal fold mimetic bioreactors have demonstrated that mechanical stimulation of vocal fold fibroblasts can lead to highly heterogeneous responses, depending on the nature and parameters of the induced vibration. These responses can either aid or impede vocal fold vibration at the physiological level. Future research is needed to determine the specific mechanical parameters that are biologically beneficial for vocal fold function.

Synchronous 3D Imaging of the Medial and Superior Vocal Fold Surface in an Excised Human Hemilarynx

OBJECTIVE

This study investigates relationships between the oscillation behavior of the medial and superior vocal fold (VF) surfaces during sustained phonation in a human cadaver hemilarynx.

METHODS

An experimental test stand synchronously captured the medial and superior VF surfaces of a human ex vivo hemilarynx during sustained phonation using two high-speed camera setups in 24 experimental settings. The 3D coordinates of the medial VF surface were reconstructed by triangulation of sewn-in marker points, while laser-based reconstruction was used for the superior VF surface. Correlation analysis and linear regression were used to quantify the connections of the mean and maximal vertical and lateral VF displacements and the VF velocities. Additionally, stepwise linear regression was used to analyze the impact of the measurement variables mean flow rate, adduction and elongation.

RESULTS

Strong linear relationships between all of the tested corresponding parameter pairs of the superior and medial VF surfaces were found (p<.001). Mean and maximum vertical displacements of the medial surface were both approximately 50% of the superior surface. The mean lateral displacements for the medial surface were 12% below the superior surface but 12% higher for the maximum values. The mean and maximum VF velocities were 32% and 36% lower for the medial surface.

CONCLUSION

The suggested multi-modal test stand allows efficient, comprehensive analysis of human hemilarynges and provides promising information about the interaction of the different VF areas and opens up the systematic analysis of multiple hemilarynges.

SIGNIFICANCE

In future, our results could integrate into ENT diagnostics using 3D laryngoscopy where the hidden medial VF surface dynamics may be predicted from the observable superior surface.

Voice Problems Among School Teachers employing the Tele-teaching Modality

OBJECTIVES

To assess the prevalence of voice problems among teachers in Riyadh during tele-teaching and examine the relationship between the Voice Handicap Index 10 (VHI10) scores and a variety of risk factors believed to be related to voice problems. We also assessed awareness of voice hygiene and therapy among teachers.

STUDY DESIGN

An observational cross-sectional study conducted using a multistage random sampling method among Riyadh school teachers who taught by tele-teaching for a minimum of one year.

METHODS

A self-assessment questionnaire which included demographic information about teachers, factors related to their teaching backgrounds, tele-teaching settings, effects of tele-teaching on the voice, medical and social histories, reports of voice and reflux symptoms, VHI10, and general knowledge about voice hygiene. This was distributed to school teachers using an SMS link through the Ministry of Education's IT department.

RESULTS

A total 495 were included in the study after exclusions. The prevalence of teachers who had significant voice problems during tele-teaching (VHI10>11) was 21.6%. Multiple risk factors significantly increased the risk of voice problems during tele-teaching. These factors included being female, teacher age, the presence of background noise from both teachers and students, loud voices, using an open camera during the teaching, stress and anxiety, allergies, respiratory disease, reflux, hearing problems, and a family history of voice problems. Only 4.6% of respondents were familiar with voice hygiene and voice therapy, but 65% believe that it is important for teachers to be knowledgeable about them.

CONCLUSIONS

Due to the lower prevalence of voice disorders among tele-teaching compared to traditional teaching methods, tele-teaching may be a viable option for teachers who have voice problems. There are still several factors influencing voice problems among tele-teachers. To attenuate potential risks, it is crucial that teachers are aware of the concepts of voice hygiene and voice therapy.

The effect of swelling on vocal fold kinematics and dynamics

Swelling in the vocal folds is caused by the local accumulation of fluid, and has been implicated as a phase in the development of phonotraumatic vocal hyperfunction and related structural pathologies, such as vocal fold nodules. It has been posited that small degrees of swelling may be protective, but large amounts may lead to a vicious cycle wherein the engorged folds lead to conditions that promote further swelling, leading to pathologies. As a first effort to explore the mechanics of vocal fold swelling and its potential role in the etiology of voice disorders, this study employs a finite-element model with swelling confined to the superficial lamina propria, which changes the volume, mass, and stiffness of the cover layer. The impacts of swelling on a number of vocal fold kinematic and damage measures, including von Mises stress, internal viscous dissipation, and collision pressure, are presented. Swelling has small but consistent effects on voice outputs, including a reduction in fundamental frequency with increasing swelling (10 Hz at 30 % swelling). Average von Mises stress decreases slightly for small degrees of swelling but increases at large magnitudes, consistent with expectations for a vicious cycle. Both viscous dissipation and collision pressure consistently increase with the magnitude of swelling. This first effort at modeling the impact of swelling on vocal fold kinematics, kinetics, and damage measures highlights the complexity with which phonotrauma can influence performance metrics. Further identification and exploration of salient candidate measures of damage and refined studies coupling swelling with local phonotrauma are expected to shed further light on the etiological pathways of phonotraumatic vocal hyperfunction.

An Evaluation of Voice Hygiene in High School Choir Singers

BACKGROUND/OBJECTIVES

Voice hygiene, including proper hydration and limiting phonotraumatic behaviors and refluxogenic foods when appropriate, is essential to voice care for singers. The purpose of this study was to evaluate voice care practices, knowledge, and factors correlating with improved voice hygiene and knowledge in adolescents, an understudied demographic within the singing population.

METHODS

214 high school choir singers, ages 14 through 18, responded to a voice health questionnaire, which included questions about phonotraumatic behaviors and other lifestyle measures regarding voice care, perceived knowledge of the voice, and potential correlational factors such as past singing experience, having a friend or family member when is a professional voice user or possesses an in-depth knowledge of the voice, use of the Internet to sing better or improve understanding of the voice, and other factors. These data were analyzed using descriptive statistics and multivariable linear regressions.

RESULTS

Results showed poor voice hygiene. Phonotraumatic behaviors were frequent and the mean daily water intake was only 5.24 cups. However, 94.4% of the respondents warmed up their voices before rehearsals or performances. The mean total knowledge score, from 0 to 6, was only 2.09. In addition, feeling that one took care of his/her voice well, and use of Internet to sing better or advance knowledge of the voice, were correlated with improved knowledge or measures of voice health.

CONCLUSION

Data suggest that improved instruction regarding voice health in the high school setting is necessary and that the Internet or social media may be valuable resources for reaching adolescent singers.

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.

LaDIVA: A neurocomputational model providing laryngeal motor control for speech acquisition and production

Many voice disorders are the result of intricate neural and/or biomechanical impairments that are poorly understood. The limited knowledge of their etiological and pathophysiological mechanisms hampers effective clinical management. Behavioral studies have been used concurrently with computational models to better understand typical and pathological laryngeal motor control. Thus far, however, a unified computational framework that quantitatively integrates physiologically relevant models of phonation with the neural control of speech has not been developed. Here, we introduce LaDIVA, a novel neurocomputational model with physiologically based laryngeal motor control. We combined the DIVA model (an established neural network model of speech motor control) with the extended body-cover model (a physics-based vocal fold model). The resulting integrated model, LaDIVA, was validated by comparing its model simulations with behavioral responses to perturbations of auditory vocal fundamental frequency (fo) feedback in adults with typical speech. LaDIVA demonstrated capability to simulate different modes of laryngeal motor control, ranging from short-term (i.e., reflexive) and long-term (i.e., adaptive) auditory feedback paradigms, to generating prosodic contours in speech. Simulations showed that LaDIVA's laryngeal motor control displays properties of motor equivalence, i.e., LaDIVA could robustly generate compensatory responses to reflexive vocal fo perturbations with varying initial laryngeal muscle activation levels leading to the same output. The model can also generate prosodic contours for studying laryngeal motor control in running speech. LaDIVA can expand the understanding of the physiology of human phonation to enable, for the first time, the investigation of causal effects of neural motor control in the fine structure of the vocal signal.

Biomechanical Models to Represent Vocal Physiology: A Systematic Review

Biomechanical modeling allows obtaining information on physical phenomena that cannot be directly observed. This study aims to review models that represent voice production. A systematic review of the literature was conducted using PubMed/Medline, SCOPUS, and IEEE Xplore databases. To select the papers, we used the protocol PRISMA Statement. A total of 53 publications were included in this review. This article considers a taxonomic classification of models found in the literature. We propose four categories in the taxonomy: (1) Models representing the Source (Vocal folds); (2) Models representing the Filter (Vocal Tract); (3) Models representing the Source - Filter Interaction; and (4) Models representing the Airflow - Source Interaction. We include a bibliographic analysis with the evolution of the publications per category. We provide an analysis of the number as well of publications in journals per year. Moreover, we present an analysis of the term occurrence and its frequency of usage, as found in the literature. In each category, different types of vocal production models are mentioned and analyzed. The models account for the analysis of evidence about aerodynamic, biomechanical, and acoustic phenomena and their correlation with the physiological processes involved in the production of the human voice. This review gives an insight into the state of the art related to the mathematical modeling of voice production, analyzed from the viewpoint of vocal physiology.

Study of spatiotemporal liquid dynamics in a vibrating vocal fold by using a self-oscillating poroelastic model

The main purpose of this study is to investigate the spatiotemporal interstitial fluid dynamics in a vibrating vocal fold. A self-oscillating poroelastic model is proposed to study the liquid dynamics in the vibrating vocal folds by treating the vocal fold tissue as a transversally isotropic, fluid-saturated, porous material. Rich spatiotemporal liquid dynamics have been found. Specifically, in the vertical direction, the liquid is transported from the inferior side to the superior side due to the propagation of the mucosal wave. In the longitudinal direction, the liquid accumulates at the anterior-posterior midpoint. However, the contact between the two vocal folds forces the accumulated liquid out laterally in a very short time span. These findings could be helpful for exploring etiology of some laryngeal pathologies, optimizing laryngeal disease treatment, and understanding hemodynamics in the vocal folds.

Laryngeal strategies to minimize vocal fold contact pressure and their effect on voice production

The goal of this study is to identify laryngeal strategies that minimize vocal fold contact pressure while producing a target sound pressure level (SPL) using a three-dimensional voice production model. The results show that while the subglottal pressure and transverse stiffness can be manipulated to reduce the peak contact pressure, such manipulations also reduce the SPL, and are thus less effective in reducing contact pressure in voice tasks targeting a specific SPL level. In contrast, changes in the initial glottal angle and vocal fold vertical thickness that reduce the contact pressure also increase the SPL. Thus, in voice tasks targeting a specific SPL, such changes in the initial glottal angle and vertical thickness also lower the subglottal pressure, which further reduces the peak contact pressure. Overall the results show that for voice tasks with a target SPL level, vocal fold contact pressure can be significantly reduced by adopting a barely abducted glottal configuration or reducing the vocal fold vertical thickness. Aerodynamic measures are effective in identifying voice production with large initial glottal angles, but by themselves alone are not useful in differentiating hyperadducted vocal folds from barely abducted vocal folds, which may be better differentiated by closed quotient and voice type measures.

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.

Toward Development of a Vocal Fold Contact Pressure Probe: Sensor Characterization and Validation Using Synthetic Vocal Fold Models

Excessive vocal fold collision pressures during phonation are considered to play a primary role in the formation of benign vocal fold lesions, such as nodules. The ability to accurately and reliably acquire intraglottal pressure has the potential to provide unique insights into the pathophysiology of phonotrauma. Difficulties arise, however, in directly measuring vocal fold contact pressures due to physical intrusion from the sensor that may disrupt the contact mechanics, as well as difficulty in determining probe/sensor position relative to the contact location. These issues are quantified and addressed through the implementation of a novel approach for identifying the timing and location of vocal fold contact, and measuring intraglottal and vocal fold contact pressures via a pressure probe embedded in the wall of a hemi-laryngeal flow facility. The accuracy and sensitivity of the pressure measurements are validated against ground truth values. Application to in vivo approaches are assessed by acquiring intraglottal and VF contact pressures using a synthetic, self-oscillating vocal fold model in a hemi-laryngeal configuration, where the sensitivity of the measured intraglottal and vocal fold contact pressure relative to the sensor position is explored.

Vocal fold contact pressure in a three-dimensional body-cover phonation model

The goal of this study is to identify vocal fold geometric and mechanical conditions that are likely to produce large contact pressure and thus high risk of vocal fold injury. Using a three-dimensional computational model of phonation, parametric simulations are performed with co-variations in vocal fold geometry and stiffness, with and without a vocal tract. For each simulation, the peak contact pressure is calculated. The results show that the subglottal pressure and the transverse stiffness of the vocal folds in the coronal plane have the largest and most consistent effect on the peak contact pressure, indicating the importance of maintaining a balance between the subglottal pressure and transverse stiffness to avoiding vocal fold injury. The presence of a vocal tract generally increases the peak contact pressure, particularly for an open-mouth vocal tract configuration. While a low degree of vocal fold approximation significantly reduces vocal fold contact pressure, for conditions of moderate and tight vocal fold approximation changes in vocal fold approximation may increase or decrease the peak contact pressure. The effects of the medial surface thickness and vocal fold stiffness along the anterior-posterior direction are similarly inconsistent and vary depending on other control parameters and the vocal tract configuration.

Influence of vocal fold cover layer thickness on its vibratory dynamics during voice production

The influence of vocal fold cover layer thickness on the flow-induced vibration and voice production was studied by using a continuum-mechanics based computational model. The cover-body thickness ratio of a three-layer vocal fold was systematically varied. The effect on the vocal fold stiffness, eigenfrequencies and eigenmodes, fundamental frequencies, glottal flow rate, vocal fold vibratory dynamics, and synchronization of the eigenmodes were analyzed by using the structure eigen analysis and flow-structure interaction simulations. It was found that the cover-body layer thickness ratio significantly affected the strength and synchronization of the eigenmodes during flow-structure interactions, and ultimately affected the fundamental frequency and vibration pattern. With the increasing cover-body thickness ratio, the strength of the wave-type higher-eigenfrequency modes increased, and that resulted in a nonlinear bifurcation of the system in which the system evolved from a regular periodic vibration to a periodic doubling vibration and then back to a regular periodic vibration with increased fundamental frequencies. During the transition, the system vibrated chaotically. Because of the increased strength of the wave-type modes, the maximum divergent angle of the glottis was also increased with the increasing cover-body thickness ratio.

Estimating Vocal Fold Contact Pressure from Raw Laryngeal High-Speed Videoendoscopy Using a Hertz Contact Model

The development of trauma-induced lesions of the vocal folds (VFs) has been linked to a high collision pressure on the VF surface. However, there are no direct methods for the clinical assessment of VF collision, thus limiting the objective assessment of these disorders. In this study, we develop a video processing technique to directly quantify the mechanical impact of the VFs using solely laryngeal kinematic data. The technique is based on an edge tracking framework that estimates the kinematic sequence of each VF edge with a Kalman filter approach and a Hertzian impact model to predict the contact force during the collision. The proposed formulation overcomes several limitations of prior efforts since it uses a more relevant VF contact geometry, it does not require calibrated physical dimensions, it is normalized by the tissue properties, and it applies a correction factor for using a superior view only. The proposed approach is validated against numerical models, silicone vocal fold models, and prior studies. A case study with high-speed videoendoscopy recordings provides initial insights between the sound pressure level and contact pressure. Thus, the proposed method has a high potential in clinical practice and could also be adapted to operate with laryngeal stroboscopic systems.

A method to directly estimate the contact pressure of the vocal folds using uncalibrated laryngeal kinematic data is presented. The approach is promising in enhancing the objective assessment of vocal function in clinical settings, especially for studying same-subject variations.

Recursive Least Squares Filtering Algorithms for On-Line Viscoelastic Characterization of Biosamples

The mechanical characterization of biological samples is a fundamental issue in biology and related fields, such as tissue and cell mechanics, regenerative medicine and diagnosis of diseases. In this paper, a novel approach for the identification of the stiffness and damping coefficients of biosamples is introduced. According to the proposed method, a MEMS-based microgripper in operational condition is used as a measurement tool. The mechanical model describing the dynamics of the gripper-sample system considers the pseudo-rigid body model for the microgripper, and the Kelvin–Voigt constitutive law of viscoelasticity for the sample. Then, two algorithms based on recursive least square (RLS) methods are implemented for the estimation of the mechanical coefficients, that are the forgetting factor based RLS and the normalised gradient based RLS algorithms. Numerical simulations are performed to verify the effectiveness of the proposed approach. Results confirm the feasibility of the method that enables the ability to perform simultaneously two tasks: sample manipulation and parameters identification.