Investigation of Vibrational Frequency of Canine Vocal Folds Using a Two-Way Fluid-Solid Interaction Analysis

Introduction Speech is an integral component of human communication, requiring the coordinated efforts of various organs to produce sound (Titze & Alipour, 2006). The glottis region, a key player in voice production, assumes a crucial role in this intricate process. As air, emanating from the lungs in a confined space, interacts with the vocal folds (VFs) within the human body, it gives rise to the creation of voice (Alipour & Vigmostad, 2012). Understanding the mechanical intricacies of this process is very important. Studying VFs in vivo situations is hard work. However, the orientation, shape and size of VFs fibers have been extracted with synchrotron X-ray microtomography. (Bailly et al., 2018) The investigation of mechanical properties of both human and animal VFs has been carried out through various methodologies in the literature. The mechanical properties of VFs have been studied using the uniaxial extension test (Alipour & Vigmostad, 2012) assuming a linear behavior, while the nonlinearity and anisotropy of VFs has been determined using a multiscale method as in Miri et al. (2013). Pipette aspiration has also been used to extract in vivo elastic properties of VFs (Scheible et al., 2023). Mechanical behavior of VFs layers in tension, compression and shear has been studied. (Cochereau et al., 2020). Fluid-structure interaction (FSI) simulations provide a valuable tool to gain a deeper understanding of voice production (Ghorbani et al. 2022). These simulations allow us to model the dynamic interplay between the VFs and air. Our research focuses on investigating the mechanical properties of canine vocal folds and utilizing these findings in an FSI simulation. Through this simulation, we aim to unravel how these mechanical properties affect voice production.Methods To investigate the mechanical properties of canine VFs, an in vitro study was conducted involving 6 mixedbreed dogs. The samples were harvested from canine cadavers euthanized for reasons unrelated to this study. In the following, the VFs were harvested and tested upon 3-4 hours post-animal sacrifice. Experimental trials were carried out using the STM-1 device (SANTAM Co.), equipped with a 100 kg load cell. Seven uniaxial tensile tests were done on each sample, with displacement rates of 1, 5, 10, 20, 40, 60, and 120 mm/min. The very slow rate of 1 mm/min was chosen to assess only elastic properties eliminating viscosity effects. Various hyperelastic models were used to fit the experimental data. Subsequently, for each model, both the mean and standard deviation (SD) were determined for the hyperelastic model parameters and their residuals. For FSI analysis we used a simplified laryngeal model as a hollow cylinder with a diameter of 50 mm and a thickness of 3 mm. The overall length of the larynx was set at 100 mm. The VFs were modeled as a circular disc with a small elliptical fissure in the midst of the cylinder section. Boundary conditions were established based on pressure differentials, with the inlet gauge pressure set at 1200 Pa and the relative pressure at the outlet set to 0. To account for the turbulent nature of airflow within the larynx, we employed the K-epsilon method to solve the motion differential equations in a two-way fluid-structure interaction simulation using ANSYS FLUENT 2021. This approach enabled us to investigate how the acquired mechanical properties of canine vocal folds affect the FSI simulations during phonation, resulting in a more comprehensive understanding of their impact. To determine the vibrational frequency of VFs, we calculated the time it took to reach maximum displacement and then quadrupled this value to obtain the period of vibration.