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 (f_o). We found that while CT and TA  activation increased the f_o 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.