Optical measurements of vocal fold tensile properties: implications for phonatory mechanics

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.

Mechanics of human vocal folds layers during finite strains in tension, compression and shear

During phonation, human vocal fold tissues are subjected to combined tension, compression and shear loading modes from small to large finite strains. Their mechanical behaviour is however still not well understood. Herein, we complete the existing mechanical database of these soft tissues, by characterising, for the first time, the cyclic and finite strains behaviour of the lamina propria and vocalis layers under these loading modes. To minimise the inter or intra-individual variability, particular attention was paid to subject each tissue sample successively to the three loadings. A non-linear mechanical behaviour is observed for all loading modes: a J-shape strain stiffening in longitudinal tension and transverse compression, albeit far less pronounced in shear, stress accommodation and stress hysteresis whatever the loading mode. In addition, recorded stress levels during longitudinal tension are much higher for the lamina propria than for the vocalis. Conversely, the responses of the lamina propria and the vocalis in transverse compression as well as transverse and longitudinal shears are of the same orders of magnitude. We also highlight the strain rate sensitivity of the tissues, as well as their anisotropic properties.

Collagen fibrillar structures in vocal fold scarring and repair using stem cell therapy: a detailed histological, immunohistochemical and atomic force microscopy study

Regenerative medicine opens new opportunities in the repair of cicatricial lesions of the vocal folds. Here, we present a thorough morphological study, with the focus on the collagen structures in the mucosa of the vocal folds, dedicated to the effects of stem cells on the vocal folds repair after cicatricial lesions. We used a conventional experimental model of a mature scar of the rabbit vocal folds, which was surgically excised with a simultaneous implantation of autologous bone marrow-derived mesenchymal stem cells (MSC) into the defect. The restoration of the vocal folds was studied 3 months postimplantation of stem cells and 6 months after the first surgery. The collagen structure assessment included histology, immunohistochemistry and atomic force microscopy (AFM) studies. According to the data of optical microscopy and AFM, as well as to immunohistochemical analysis, MSC implantation into the vocal fold defect leads not only to the general reduction of scarring, normal ratio of collagens type I and type III, but also to a more complete restoration of architecture and ultrastructure of collagen fibres in the mucosa, as compared to the control. The collagen structures in the scar tissue in the vocal folds with implanted MSC are more similar to those in the normal mucosa of the vocal folds than to those of the untreated scars. AFM has proven to be an instrumental technique in the assessment of the ultrastructure restoration in such studies. LAY DESCRIPTION: Regenerative medicine opens new opportunities in the repair of the vocal fold scars. Because collagen is a main component in the vocal fold mucosa responsible for the scar formation and repair, we focus on the collagen structures in the mucosa of the vocal folds, using a thorough morphological study based on histology and atomic force microscopy (AFM). Atomic force microscopy is a scanning microscopic technique which allows revealing the internal structure of a tissue with a resolution up to nanometres. We used a conventional experimental model of a mature scar of the rabbit vocal folds, surgically excised and treated with a mesenchymal stem cells transplant. Our morphological study, primarily AFM, explicitly shows that the collagen structures in the scarred vocal folds almost completely restore after the stem cell treatment. Thus, the modern microscopic methods, and especially AFM are instrumental tools for monitoring the repair of the vocal folds scars.

In situ vocal fold properties and pitch prediction by dynamic actuation of the songbird syrinx

The biomechanics of sound production forms an integral part of the neuromechanical control loop of avian vocal motor control. However, we critically lack quantification of basic biomechanical parameters describing the vocal organ, the syrinx, such as material properties of syringeal elements, forces and torques exerted on, and motion of the syringeal skeleton during song. Here, we present a novel marker-based 3D stereoscopic imaging technique to reconstruct 3D motion of servo-controlled actuation of syringeal muscle insertions sites in vitro and focus on two muscles controlling sound pitch. We furthermore combine kinematic analysis with force measurements to quantify elastic properties of sound producing medial labia (ML). The elastic modulus of the zebra finch ML is 18 kPa at 5% strain, which is comparable to elastic moduli of mammalian vocal folds. Additionally ML lengthening due to musculus syringealis ventralis (VS) shortening is intrinsically constraint at maximally 12% strain. Using these values we predict sound pitch to range from 350–800 Hz by VS modulation, corresponding well to previous observations. The presented methodology allows for quantification of syringeal skeleton motion and forces, acoustic effects of muscle recruitment, and calibration of computational birdsong models, enabling experimental access to the entire neuromechanical control loop of vocal motor control.

Computational modeling of vibration-induced systemic hydration of vocal folds over a range of phonation conditions

Predicting phonation conditions that are benign to voice health remains a biomechanically relevant problem. Our objective is to provide insight into vocal fold (VF) hydration based on continuum-based VF models that are able to compute VF stresses during phonation and a scheme for the extraction and generalization of such computational data based on the principle of linear superposition. Because VF tissue is poroelastic, spatial gradients of VF hydrostatic stresses computed for a given phonation condition determine VF interstitial fluid flow. The present approach transforms, based on linear superposition principles, the computed interstitial fluid velocities at the particular phonation to those at an arbitrary phonation condition. Intersititial fluid flow characteristics for a range of phonation conditions are compared. For phonation conditions with no or moderate collision, no dehydration per vibration cycle is predicted throughout the VF. For more severe collision conditions, tissue dehydration is restricted to a region close to the glottal surface. Interstitial fluid displacement in the VF is found to be heterogeneous and strongly dependent on the phonation condition. A phonation condition is found to exist for which dehydration peaks. The proposed method significantly expands the scope and relevance of conducting isolated numerical simulations of VF vibration.

Mechanical properties vary for different regions of the finger extensor apparatus

The extensor apparatus, an aponeurosis that covers the dorsal side of each finger, transmits force from a number of musculotendons to the phalanges. Multiple tendons integrate directly into the structure at different sites and the extensor apparatus attaches to the phalanges at multiple points. Thus, prediction of the force distribution within the extensor apparatus, or hood, and the transmission to the phalanges is challenging, especially as knowledge of the underlying mechanical properties of the tissue is limited. We undertook quantification of some of these properties through material testing of cadaver specimens. We punched samples at specified locations from 19 extensor hood specimens. Material testing was performed to failure for each sample with a custom material testing device. Testing revealed significant differences in ultimate load, ultimate strain, thickness, and tangent modulus along the length of the extensor hood. Specifically, thickness, ultimate load, and ultimate strain were greater in the more proximal sections of the extensor hood, while the tangent modulus was greater in the more distal sections. The variations in mechanical properties within the hood may impact prediction of force transmission and, thus, should be considered when modeling the action of the extensor apparatus. Across the extensor hood, tangent modulus values were substantially smaller than values reported for other soft tissues, such as the Achilles tendon and knee ligaments, while ultimate strains were much greater. Thus, the tissue in the extensor apparatus seems to have greater elasticity, which should be modeled accordingly.

Material and shape optimization for multi-layered vocal fold models using transient loadings

Commonly applied models to study vocal fold vibrations in combination with air flow distributions are self-sustained physical models of the larynx consisting of artificial silicone vocal folds. Choosing appropriate mechanical parameters and layer geometries for these vocal fold models while considering simplifications due to manufacturing restrictions is difficult but crucial for achieving realistic behavior. In earlier work by Schmidt et al. [J. Acoust. Soc. Am. 129, 2168-2180 (2011)], the authors presented an approach in which material parameters of a static numerical vocal fold model were optimized to achieve an agreement of the displacement field with data retrieved from hemilarynx experiments. This method is now generalized to a fully transient setting. Moreover in addition to the material parameters, the extended approach is capable of finding optimized layer geometries. Depending on chosen material restriction, significant modifications of the reference geometry are predicted. The additional flexibility in the design space leads to a significantly more realistic deformation behavior. At the same time, the predicted biomechanical and geometrical results are still feasible for manufacturing physical vocal fold models consisting of several silicone layers. As a consequence, the proposed combined experimental and numerical method is suited to guide the construction of physical vocal fold models.

The anisotropic hyperelastic biomechanical response of the vocal ligament and implications for frequency regulation: a case study

One of the primary mechanisms to vary one's vocal frequency is through vocal fold length changes. As stress and deformation are linked to each other, it is hypothesized that the anisotropy in the biomechanical properties of the vocal fold tissue would affect the phonation characteristics. A biomechanical model of vibrational frequency rise during vocal fold elongation is developed which combines an advanced biomechanical characterization protocol of the vocal fold tissue with continuum beam models. Biomechanical response of the tissue is related to a microstructurally informed, anisotropic, nonlinear hyperelastic constitutive model. A microstructural characteristic (the dispersion of collagen) was represented through a statistical orientation function acquired from a second harmonic generation image of the vocal ligament. Continuum models of vibration were constructed based upon Euler-Bernoulli and Timoshenko beam theories, and applied to the study of the vibration of a vocal ligament specimen. From the natural frequency predictions in dependence of elongation, two competing processes in frequency control emerged, i.e., the applied tension raises the frequency while simultaneously shear deformation lowers the frequency. Shear becomes much more substantial at higher modes of vibration and for highly anisotropic tissues. The analysis was developed as a case study based on a human vocal ligament specimen.

A canonical biomechanical vocal fold model

The present article aimed at constructing a canonical geometry of the human vocal fold (VF) from subject-specific image slice data. A computer-aided design approach automated the model construction. A subject-specific geometry available in literature, three abstractions (which successively diminished in geometric detail) derived from it, and a widely used quasi two-dimensional VF model geometry were used to create computational models. The first three natural frequencies of the models were used to characterize their mechanical response. These frequencies were determined for a representative range of tissue biomechanical properties, accounting for underlying VF histology. Compared with the subject-specific geometry model (baseline), a higher degree of abstraction was found to always correspond to a larger deviation in model frequency (up to 50% in the relevant range of tissue biomechanical properties). The model we deemed canonical was optimally abstracted, in that it significantly simplified the VF geometry compared with the baseline geometry but can be recalibrated in a consistent manner to match the baseline response. Models providing only a marginally higher degree of abstraction were found to have significant deviation in predicted frequency response. The quasi two-dimensional model presented an extreme situation: it could not be recalibrated for its frequency response to match the subject-specific model. This deficiency was attributed to complex support conditions at anterior-posterior extremities of the VFs, accentuated by further issues introduced through the tissue biomechanical properties. In creating canonical models by leveraging advances in clinical imaging techniques, the automated design procedure makes VF modeling based on subject-specific geometry more realizable.

Pipette aspiration applied to the characterization of nonhomogeneous, transversely isotropic materials used for vocal fold modeling

The etiology and treatment of voice disorders are still not completely understood. Since the vibratory characteristics of vocal folds are strongly influenced by both anatomy and mechanical material properties, measurement methods to analyze the material behavior of vocal fold tissue are required. Due to the limited life time of real tissue in the laboratory, synthetic models are often used to study vocal fold vibrations. In this paper we focus on two topics related to synthetic and real vocal fold materials. First, because certain tissues within the human vocal folds are transversely isotropic, a fabrication process for introducing this characteristic in commonly used vocal fold modeling materials is presented. Second, the pipette aspiration technique is applied to the characterization of these materials. By measuring the displacement profiles of stretched specimens that exhibit varying degrees of transverse isotropy, it is shown that local anisotropy can be quantified using a parameter describing the deviation from an axisymmetric profile. The potential for this technique to characterize homogeneous, anisotropic materials, including soft biological tissues such as those found in the human vocal folds, is supplemented by a computational study.

Empirical measurements of biomechanical anisotropy of the human vocal fold lamina propria

The vocal folds are known to be mechanically anisotropic due to the microstructural arrangement of fibrous proteins such as collagen and elastin in the lamina propria. Even though this has been known for many years, the biomechanical anisotropic properties have rarely been experimentally studied. We propose that an indentation procedure can be used with uniaxial tension in order to obtain an estimate of the biomechanical anisotropy within a single specimen. Experiments were performed on the lamina propria of three male and three female human vocal folds dissected from excised larynges. Two experiments were conducted: each specimen was subjected to cyclic uniaxial tensile loading in the longitudinal (i.e., anterior-posterior) direction, and then to cyclic indentation loading in the transverse (i.e., medial-lateral) direction. The indentation experiment was modeled as contact on a transversely isotropic half-space using the Barnett-Lothe tensors. The longitudinal elastic modulus E(L) was computed from the tensile test, and the transverse elastic modulus E(T) and longitudinal shear modulus G(L) were obtained by inverse analysis of the indentation force-displacement response. It was discovered that the average of E(L) /E(T) was 14 for the vocal ligament and 39 for the vocal fold cover specimens. Also, the average of E(L)/G(L), a parameter important for models of phonation, was 28 for the vocal ligament and 54 for the vocal fold cover specimens. These measurements of anisotropy could contribute to more accurate models of fundamental frequency regulation and provide potentially better insights into the mechanics of vocal fold vibration.

Could spatial heterogeneity in human vocal fold elastic properties improve the quality of phonation?

The physical mechanisms leading to the acoustic and perceptual qualities of voice are not well understood. This study examines the spatial distribution of biomechanical properties in human vocal folds and explores the consequences of these properties on phonation. Vocal fold lamina propria specimens isolated from nine excised human male larynges were tested in uniaxial tension (six from non-smokers, three from smokers). An optical method was employed to determine the local stretch, from which the elastic modulus of three segments in the anterior-posterior direction was calculated. Several specimens exhibited a significant heterogeneity in the modulus with the middle segment stiffer than the other segments. It was concluded that such modulus gradients are stronger in specimens from non-smokers than smokers. To understand the functional implications of a modulus gradient, the first eigenmode of vibration was calculated with a finite element model. With a modulus gradient, the vocal fold's eigenmode deflection was spread along the anterior-posterior length, whereas for a homogeneous modulus distribution, the deflection was more focused around the mid-coronal plane. Consequently, the strong modulus gradient may enable more complete glottal closure, which is important for normal phonation, while a more homogeneous modulus may be responsible for poor glottal closure and a perceived "breathy" voice.