Triangular body-cover model of the vocal folds with coordinated activation of the five intrinsic laryngeal muscles

A computational study of the influence of thyroarytenoid and cricothyroid muscle interaction on vocal fold dynamics in an MRI-based human laryngeal model

A human laryngeal model, incorporating all the cartilages and the intrinsic muscles, was reconstructed based on MRI data. The vocal fold was represented as a multilayer structure with detailed inner components. The activation levels of the thyroarytenoid (TA) and cricothyroid (CT) muscles were systematically varied from zero to full activation allowing for the analysis of their interaction and influence on vocal fold dynamics and glottal flow. The finite element method was employed to calculate the vocal fold dynamics, while the one-dimensional Bernoulli equation was utilized to calculate the glottal flow. The analysis was focused on the muscle influence on the fundamental frequency (fo). We found that while CT and TA  activation increased the fo in most of the conditions, TA activation resulted in a frequency drop when it was moderately activated. We show that this frequency drop was associated with the sudden increase of the vertical motion when the vibration transited from involving the whole tissue to mainly in the cover layer. The transition of the vibration pattern was caused by the increased body-cover stiffness ratio that resulted from TA activation.

Asymmetric triangular body-cover model of the VFs 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 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. The results highlight the antagonistic effect between the thyroarytenoid and cricothyroid muscles on the elastic and mass components of the vocal folds, as well as the impact on the vocal process from the imbalance in the lateral cricoarytenoid and interarytenoid adductor muscles. Measurements of amplitude and phase asymmetry were employed to emulate 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 versatility of muscle imbalance in representing phonatory scenarios and its potential for characterizing asymmetry in vocal fold vibration.

An Euler-Bernoulli-type beam model of the vocal folds for describing curved and incomplete glottal closure patterns

Incomplete glottal closure is a laryngeal configuration wherein the glottis is not fully obstructed prior to phonation. It has been linked to inefficient voice production and voice disorders. Various incomplete glottal closure patterns can arise and the mechanisms driving them are not well understood. In this work, we introduce an Euler-Bernoulli composite beam vocal fold (VF) model that produces qualitatively similar incomplete glottal closure patterns as those observed in experimental and high-fidelity numerical studies, thus offering insights into the potential underlying physical mechanisms. Refined physiological insights are pursued by incorporating the beam model into a VF posturing model that embeds the five intrinsic laryngeal muscles. Analysis of the combined model shows that co-activating the lateral cricoarytenoid (LCA) and interarytenoid (IA) muscles without activating the thyroarytenoid (TA) muscle results in a bowed (convex) VF geometry with closure at the posterior margin only; this is primarily attributed to the reactive moments at the anterior VF margin. This bowed pattern can also arise during VF compression (due to extrinsic laryngeal muscle activation for example), wherein the internal moment induced passively by the TA muscle tissue is the predominant mechanism. On the other hand, activating the TA muscle without incorporating other adductory muscles results in anterior and mid-membranous glottal closure, a concave VF geometry, and a posterior glottal opening driven by internal moments induced by TA muscle activation. In the case of initial full glottal closure, the posterior cricoarytenoid (PCA) muscle activation cancels the adductory effects of the LCA and IA muscles, resulting in a concave VF geometry and posterior glottal opening. Furthermore, certain maneuvers involving co-activation of all adductory muscles result in an hourglass glottal shape due to a reactive moment at the anterior VF margin and moderate internal moment induced by TA muscle activation. These findings have implications regarding potential laryngeal maneuvers in patients with voice disorders involving imbalances or excessive tension in the laryngeal muscles such as muscle tension dysphonia.

Modeling the influence of the extrinsic musculature on phonation

Neck muscles play important roles in various physiological tasks, including swallowing, head stabilization, and phonation. The mechanisms by which neck muscles influence phonation are not well understood, with conflicting reports on the change in fundamental frequency for ostensibly the same neck muscle activation scenarios. In this work, we introduce a reduced-order muscle-controlled vocal fold model, comprising both intrinsic muscle control and extrinsic muscle effects. The model predicts that when the neck muscles pull the thyroid cartilage in the superior-anterior direction (with a sufficiently large anterior component), inferior direction, or inferior-anterior direction, tension in the vocal folds increases, leading to fundamental frequency rise during sustained phonation. On the other hand, pulling in the superior direction, superior-posterior direction, or inferior-posterior direction (with a sufficiently large posterior component) tends to decrease vocal fold tension and phonation fundamental frequency. Varying the pulling force location alters the posture and phonation biomechanics, depending on the force direction. These findings suggest potential roles of particular neck muscles in modulating phonation fundamental frequency, with implications for vocal hyperfunction.

An Euler-Bernoulli-Type Beam Model of the Vocal Folds for Describing Curved and Incomplete Glottal Closure Patterns

Incomplete glottal closure is a laryngeal configuration wherein the glottis is not fully obstructed prior to phonation. In this work, we introduce an Euler-Bernoulli composite beam vocal fold (VF) model that produces qualitatively similar incomplete glottal closure patterns as those observed in experimental and high-fidelity numerical studies, thus offering insights in to the potential underlying physical mechanisms. Refined physiological insights are pursued by incorporating the beam model into a VF posturing model that embeds the five intrinsic laryngeal muscles. Analysis of the combined model shows that co-activating the lateral cricoarytenoid (LCA) and interarytenoid (IA) muscles without activating the thyroarytenoid (TA) muscle results in a bowed (convex) VF geometry with closure at the posterior margin only; this is primarily attributed to the reactive moments at the anterior VF margin. This bowed pattern can also arise during VF compression (due to extrinsic laryngeal muscle activation for example), wherein the internal moment induced passively by the TA muscle tissue is the predominant mechanism. On the other hand, activating the TA muscle without incorporating other adductory muscles results in anterior and mid-membranous glottal closure, a concave VF geometry, and a posterior glottal opening driven by internal moments induced by TA muscle activation. In the case of initial full glottal closure, the posterior cricoarytenoid (PCA) muscle activation cancels the adductory effects of the LCA and IA muscles, resulting in a concave VF geometry and posterior glottal opening. Furthermore, certain maneuvers involving co-activation of all adductory muscles result in an hourglass glottal shape due to a reactive moment at the anterior VF margin and moderate internal moment induced by TA muscle activation.

Exploring the mechanics of fundamental frequency variation during phonation onset

Fundamental frequency patterns during phonation onset have received renewed interest due to their promising application in objective classification of normal and pathological voices. However, the associated underlying mechanisms producing the wide array of patterns observed in different phonetic contexts are not yet fully understood. Herein, we employ theoretical and numerical analyses in an effort to elucidate the potential mechanisms driving opposing frequency patterns for initial/isolated vowels versus vowels preceded by voiceless consonants. Utilizing deterministic lumped-mass oscillator models of the vocal folds, we systematically explore the roles of collision and muscle activation in the dynamics of phonation onset. We find that an increasing trend in fundamental frequency, as observed for initial/isolated vowels, arises naturally through a progressive increase in system stiffness as collision intensifies as onset progresses, without the need for time-varying vocal fold tension or changes in aerodynamic loading. In contrast, reduction in cricothyroid muscle activation during onset is required to generate the decrease in fundamental frequency observed for vowels preceded by voiceless consonants. For such phonetic contexts, our analysis shows that the magnitude of reduction in the cricothyroid muscle activation and the activation level of the thyroarytenoid muscle are potential factors underlying observed differences in (relative) fundamental frequency between speakers with healthy and hyperfunctional voices. This work highlights the roles of sometimes competing laryngeal factors in producing the complex array of observed fundamental frequency patterns during phonation onset.