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.     

Hyaluronic acid-based microgels and microgel networks for vocal fold regeneration

Vocal fold scarring disrupts the viscoelastic properties of the lamina propria that are critical for normal phonation. There is a clinical need for the development of advanced biomaterials that approximate the mechanical properties of the lamina propria for in vivo vocal fold regeneration. We have developed hyaluronic acid (HA)-based microgels and cross-linked microgel networks with tunable degradation and mechanical properties. HA microgels were prepared by cross-linking HA derivatives carrying hydrazide (HAADH) and aldehyde (HAALD) functionalities within the inverse emulsion droplets. Alternatively, poly(ethylene glycol) dialdehyde (PEGDiALD) was employed in place of HAALD. Microgels based on HAADH/HAALD are more resistant to enzymatic degradation than those generated from HAADH/PEGDiALD. In vitro cytotoxicity studies using vocal fold fibroblasts indicate that microgels synthesized from HAADH/HAALD are essentially nontoxic, whereas microgels derived from HAADH/PEGDiALD exhibit certain adverse effects on the cultured cells at high concentration (> or =2 mg/mL). These microgels exhibit residual functional groups that can be used as reactive handles for covalent conjugation of therapeutic molecules. The presence of residual functional groups also allows for subsequent cross-linking of the microgels with other reactive polymers, giving rise to doubly cross-linked networks (DXNs) with tunable viscoelasticity. Mechanical measurements using a torsional wave apparatus indicate that HA-based DXNs exhibit elastic moduli that are similar to those of vocal fold lamina propria at frequencies close to the range of human phonation. These HA-based microgel systems are promising candidates for the treatment of vocal fold scarring, not just as biocompatible filler materials, but as smart entities that can repair focal defects, smooth the vocal fold margin, and potentially soften and dissolve scar tissue.     

Vocal fold tissue failure: preliminary data and constitutive modeling

In human voice production (phonation), linear small-amplitude vocal fold oscillation occurs only under restricted conditions. Physiologically, phonation more often involves large-amplitude oscillation associated with tissue stresses and strains beyond their linear viscoelastic limits, particularly in the lamina propria extracellular matrix (ECM). This study reports some preliminary measurements of tissue deformation and failure response of the vocal fold ECM under large-strain shear The primary goal was to formulate and test a novel constitutive model for vocal fold tissue failure, based on a standard-linear cohesive-zone (SL-CZ) approach. Tissue specimens of the sheep vocal fold mucosa were subjected to torsional deformation in vitro, at constant strain rates corresponding to twist rates of 0.01, 0.1, and 1.0 rad/s. The vocal fold ECM demonstrated nonlinear stress-strain and rate-dependent failure response with a failure strain as low as 0.40 rad. A finite-element implementation of the SL-CZ model was capable of capturing the rate dependence in these preliminary data, demonstrating the model's potential for describing tissue failure. Further studies with additional tissue specimens and model improvements are needed to better understand vocal fold tissue failure.     

Regeneration of Vocal Fold Mucosa Using Tissue-Engineered Structures with Oral Mucosal Cells

Objectives

Scarred vocal folds result in irregular vibrations during phonation due to stiffness of the vocal fold mucosa. To date, a completely satisfactory corrective procedure has yet to be achieved. We hypothesize that a potential treatment option for this disease is to replace scarred vocal folds with organotypic mucosa. The purpose of this study is to regenerate vocal fold mucosa using a tissue-engineered structure with autologous oral mucosal cells.

Study Design

Animal experiment using eight beagles (including three controls).

Methods

A 3 mm by 3 mm specimen of canine oral mucosa was surgically excised and divided into epithelial and subepithelial tissues. Epithelial cells and fibroblasts were isolated and cultured separately. The proliferated epithelial cells were co-cultured on oriented collagen gels containing the proliferated fibroblasts for an additional two weeks. The organotypic cultured tissues were transplanted to the mucosa-deficient vocal folds. Two months after transplantation, vocal fold vibrations and morphological characteristics were observed.

Results

A tissue-engineered vocal fold mucosa, consisting of stratified epithelium and lamina propria, was successfully fabricated to closely resemble the normal layered vocal fold mucosa. Laryngeal stroboscopy revealed regular but slightly small mucosal waves at the transplanted site. Immunohistochemically, stratified epithelium expressed cytokeratin, and the distributed cells in the lamina propria expressed vimentin. Elastic Van Gieson staining revealed a decreased number of elastic fibers in the lamina propria of the transplanted site.

Conclusion

The fabricated mucosa with autologous oral mucosal cells successfully restored the vocal fold mucosa. This reconstruction technique could offer substantial clinical advantages for treating intractable diseases such as scarring of the vocal folds.     

Cervids with different vocal behavior demonstrate different viscoelastic properties of their vocal folds

The authors test the hypothesis that vocal fold morphology and biomechanical properties covary with species‐specific vocal function. They investigate mule deer (Odocoileus hemionus) vocal folds, building on, and extending data on a related cervid, the Rocky Mountain elk (Cervus elaphus nelsoni). The mule deer, in contrast to the elk, is a species with relatively little vocal activity in adult animals. Mule deer and elk vocal folds show the typical three components of the mammalian vocal fold (epithelium, lamina propria and thyroarytenoid muscle). The vocal fold epithelium and the lamina propria were investigated in two sets of tensile tests. First, creep rupture tests demonstrated that ultimate stress in mule deer lamina propria is of the same magnitude as in elk. Second, cyclic loading tests revealed similar elastic moduli for the vocal fold epithelium in mule deer and elk. The elastic modulus of the lamina propria is also similar between the two species in the low‐strain region, but differs at strains larger than 0.3. Sex differences in the stress–strain response, which have been reported for elk and human vocal folds, were not found for mule deer vocal folds. The laminae propriae in mule deer and elk vocal folds are comparatively large. In general, a thick and uniformly stiff lamina propria does not self‐oscillate well, even when high subglottic pressure is applied. If the less stiff vocal fold seen in elk is associated with a differentiated lamina propria it would allow the vocal fold to vibrate at high tension and high subglottic pressure. The results of this study support the hypothesis that viscoelastic properties of vocal folds varies with function and vocal behavior. J. Morphol., 2010. © 2009 Wiley‐Liss, Inc.     

Paired versus two-group experimental design for rheological studies of vocal fold tissues

Vibratory function of the vocal folds is largely determined by the rheological properties or viscoelastic shear properties of the vocal fold lamina propria. To date, investigation of the sample size estimation and statistical experimental design for vocal fold rheological studies is nonexistent. The current work provides the closed-form sample size formulas for two major study designs (i.e. paired and two-group designs) in vocal fold research. Our results demonstrated that the paired design could greatly increase the statistical power compared to the two-group design. By comparing the variance of estimated treatment effect, this study also confirms that ignoring within-subject and within-vocal fold correlations during rheological data analysis will likely increase type I errors. Finally, viscoelastic shear properties of intact and scarred rabbit vocal fold lamina propria were measured and used to illustrate theoretical findings in a realistic scenario and project sample size requirement for future studies.     

Predicting Achievable Fundamental Frequency Ranges in Vocalization Across Species

Vocal folds are used as sound sources in various species, but it is unknown how vocal fold morphologies are optimized for different acoustic objectives. Here we identify two main variables affecting range of vocal fold vibration frequency, namely vocal fold elongation and tissue fiber stress. A simple vibrating string model is used to predict fundamental frequency ranges across species of different vocal fold sizes. While average fundamental frequency is predominantly determined by vocal fold length (larynx size), range of fundamental frequency is facilitated by (1) laryngeal muscles that control elongation and by (2) nonlinearity in tissue fiber tension. One adaptation that would increase fundamental frequency range is greater freedom in joint rotation or gliding of two cartilages (thyroid and cricoid), so that vocal fold length change is maximized. Alternatively, tissue layers can develop to bear a disproportionate fiber tension (i.e., a ligament with high density collagen fibers), increasing the fundamental frequency range and thereby vocal versatility. The range of fundamental frequency across species is thus not simply one-dimensional, but can be conceptualized as the dependent variable in a multi-dimensional morphospace. In humans, this could allow for variations that could be clinically important for voice therapy and vocal fold repair. Alternative solutions could also have importance in vocal training for singing and other highly-skilled vocalizations.     

Synthetic,Multi-Layer,Self-Oscillating Vocal Fold Model Fabrication

Sound for the human voice is produced via flow-induced vocal fold vibration. The vocal folds consist of several layers of tissue, each with differing material properties ¹. Normal voice production relies on healthy tissue and vocal folds, and occurs as a result of complex coupling between aerodynamic, structural dynamic, and acoustic physical phenomena. Voice disorders affect up to 7.5 million annually in the United States alone ² and often result in significant financial, social, and other quality-of-life difficulties. Understanding the physics of voice production has the potential to significantly benefit voice care, including clinical prevention, diagnosis, and treatment of voice disorders.Existing methods for studying voice production include in vivo experimentation using human and animal subjects, in vitro experimentation using excised larynges and synthetic models, and computational modeling. Owing to hazardous and difficult instrument access, in vivo experiments are severely limited in scope. Excised larynx experiments have the benefit of anatomical and some physiological realism, but parametric studies involving geometric and material property variables are limited. Further, they are typically only able to be vibrated for relatively short periods of time (typically on the order of minutes).Overcoming some of the limitations of excised larynx experiments, synthetic vocal fold models are emerging as a complementary tool for studying voice production. Synthetic models can be fabricated with systematic changes to geometry and material properties, allowing for the study of healthy and unhealthy human phonatory aerodynamics, structural dynamics, and acoustics. For example, they have been used to study left-right vocal fold asymmetry ^3,4, clinical instrument development ⁵, laryngeal aerodynamics ^6-9, vocal fold contact pressure ¹⁰, and subglottal acoustics ¹¹ (a more comprehensive list can be found in Kniesburges et al. ¹²)Existing synthetic vocal fold models, however, have either been homogenous (one-layer models) or have been fabricated using two materials of differing stiffness (two-layer models). This approach does not allow for representation of the actual multi-layer structure of the human vocal folds ¹ that plays a central role in governing vocal fold flow-induced vibratory response. Consequently, one- and two-layer synthetic vocal fold models have exhibited disadvantages ^3,6,8 such as higher onset pressures than what are typical for human phonation (onset pressure is the minimum lung pressure required to initiate vibration), unnaturally large inferior-superior motion, and lack of a "mucosal wave" (a vertically-traveling wave that is characteristic of healthy human vocal fold vibration).In this paper, fabrication of a model with multiple layers of differing material properties is described. The model layers simulate the multi-layer structure of the human vocal folds, including epithelium, superficial lamina propria (SLP), intermediate and deep lamina propria (i.e., ligament; a fiber is included for anterior-posterior stiffness), and muscle (i.e., body) layers ¹. Results are included that show that the model exhibits improved vibratory characteristics over prior one- and two-layer synthetic models, including onset pressure closer to human onset pressure, reduced inferior-superior motion, and evidence of a mucosal wave.     

声带组织切片方法的比较

目的:探讨犬声带冠状位切片与水平位切片各自的特点,为声带实验提供合适的切片方法。方法:家犬4只,2只取材后行冠状位石蜡切片,2只取材后行水平位石蜡切片。通过HE染色观察声带固有层的一般组织结构,Masson三色染色观察固有层中胶原的排列情况。结果:HE染色示冠状位、水平位切片均可见声带表面被覆复层鳞状上皮,固有层内有大量排列紧密的纤维组织,纤维组织中夹杂少量腺体,固有层下方为肌层。冠状位切片可观察声带某一点冠状面固有层的情况,若观察整个声带的情况需声带连续切片;水平位切片可在一张切片中观察到前联合、声带膜部及声带突部位的固有层情况,解剖标志明显,利于定位。Masson三色染色示冠状位、水平位切片均可见固有层浅层有较细的胶原纤维束,中层有较粗的纤维束与较细的纤维束交织排列,深层纤维束排列更紧密。结论:冠状位切片可观察声带某一点冠状面固有层的整体情况,水平位切片可在一张切片中观察到前联合、膜部及声带突部位的固有层情况。     

Mechanotransduction of vocal fold fibroblasts and mesenchymal stromal cells in the context of the vocal fold mechanome

The design of cell-based therapies for vocal fold tissue engineering requires an understanding of how cells adapt to the dynamic mechanical forces found in the larynx. Our objective was to compare mechanotransductive processes in therapeutic cell candidates (mesenchymal stromal cells from adipose tissue and bone marrow, AT-MSC and BM-MSC) to native cells (vocal fold fibroblasts-VFF) in the context of vibratory strain. A bioreactor was used to expose VFF, AT-MSC, and BM-MSC to axial tensile strain and vibration at human physiological levels. Microarray, an empirical Bayes statistical approach, and geneset enrichment analysis were used to identify significant mechanotransductive pathways associated with the three cell types and three mechanical conditions. Two databases (Gene Ontology, Kyoto Encyclopedia of Genes and Genomes) were used for enrichment analyses. VFF shared more mechanotransductive pathways with BM-MSC than with AT-MSC. Gene expression that appeared to distinguish the vibratory strain condition from polystyrene condition for these two cells types related to integrin activation, focal adhesions, and lamellipodia activity, suggesting that vibratory strain may be associated with cytoarchitectural rearrangement, cell reorientation, and extracellular matrix remodeling. In response to vibration and tensile stress, BM-MSC better mimicked VFF mechanotransduction than AT-MSC, providing support for the consideration of BM-MSC as a cell therapy for vocal fold tissue engineering. Future research is needed to better understand the sorts of physical adaptations that are afforded to vocal fold tissue as a result of focal adhesions, integrins, and lamellipodia, and how these adaptations could be exploited for tissue engineering.     

Influence of asymmetric stiffness on the structural and aerodynamic response of synthetic vocal fold models

The influence of asymmetric vocal fold stiffness on voice production was evaluated using life-sized, self-oscillating vocal fold models with an idealized geometry based on the human vocal folds. The models were fabricated using flexible, materially-linear silicone compounds with Young's modulus values comparable to that of vocal fold tissue. The models included a two-layer design to simulate the vocal fold layered structure. The respective Young's moduli of elasticity of the “left” and “right” vocal fold models were varied to create asymmetric conditions. High-speed videokymography was used to measure maximum vocal fold excursion, vibration frequency, and left–right phase shift, all of which were significantly influenced by asymmetry. Onset pressure, a measure of vocal effort, increased with asymmetry. Particle image velocimetry (PIV) analysis showed significantly greater skewing of the glottal jet in the direction of the stiffer vocal fold model. Potential applications to various clinical conditions are mentioned, and suggestions for future related studies are presented.     

Reinke's Edema: investigations on the role of MIB-1 and hepatocyte growth factor

Reinke''s edema is a benign disease of the human vocal fold, which mainly affects the sub-epithelial layer of the vocal fold. Microscopic observations show a strongly oedematous epithelium with loosened intercellular junctions, a disruption of the extracellular connections between mucosal epithelium and connective tissue, closely adherent to the thyroarytenoid muscle. Thickening of the basal layer of epithelium, known as Reinke''s space, high deposition of fibronectin and chronic inflammatory infiltration it is also visible. We analyzed, together with the hepatocyte growth factor (HGF), the expression level of MIB-1 in samples harvested from patients affected by Reinke''s edema, in order to define its biological role and consider it as a possible prognostic factor in the follow-up after surgical treatment. We observed a moderate expression of HGF in the lamina propria of the human vocal fold and in the basal membrane of the mucosal epithelium. Our finding suggests that this growth factor acts as an anti - fibrotic agent in Reinke''s space and affects the fibronectin deposition in the lamina propria. MIB-1, on the contrary, showed a weak expression in the basement membrane of the mucosal epithelium and a total absence in the lamina propria deep layer, thus suggesting that only the superficial layer is actively involved in the reparatory process with a high regenerative capacity, together with a high deposition of fibronectin. The latter is necessary for the cellular connections reconstruction, after the inflammatory infiltration.Key words: Reinke''s edema, Ki67, HGF, fibronectin, immunohistochemistry, extracellular matrix.     

Cross-sample validation provides enhanced proteome coverage in rat vocal fold mucosa

The vocal fold mucosa is a biomechanically unique tissue comprised of a densely cellular epithelium, superficial to an extracellular matrix (ECM)-rich lamina propria. Such ECM-rich tissues are challenging to analyze using proteomic assays, primarily due to extensive crosslinking and glycosylation of the majority of high M(r) ECM proteins. In this study, we implemented an LC-MS/MS-based strategy to characterize the rat vocal fold mucosa proteome. Our sample preparation protocol successfully solubilized both proteins and certain high M(r) glycoconjugates and resulted in the identification of hundreds of mucosal proteins. A straightforward approach to the treatment of protein identifications attributed to single peptide hits allowed the retention of potentially important low abundance identifications (validated by a cross-sample match and de novo interpretation of relevant spectra) while still eliminating potentially spurious identifications (global single peptide hits with no cross-sample match). The resulting vocal fold mucosa proteome was characterized by a wide range of cellular and extracellular proteins spanning 12 functional categories.     

The response of vocal fold fibroblasts and mesenchymal stromal cells to vibration

Illumination of cellular changes caused by mechanical forces present within the laryngeal microenvironment may well guide strategies for tissue engineering the vocal fold lamina propria. The purpose of this study was to compare the response of human vocal fold fibroblasts (hVFF) and bone marrow mesenchymal stem cells (BM-MSC) to vibratory stimulus. In order to study these effects, a bioreactor capable of vibrating two cell seeded substrates was developed. The cell seeded substrates contact each other as a result of the sinusoidal frequency, producing a motion similar to the movement of true vocal folds. Utilizing this bioreactor, hVFF and BM-MSC were subjected to 200 Hz vibration and 20% strain for 8 hours. Immunohistochemistry (Ki-67 and TUNEL) was performed to examine cell proliferation and apoptosis respectively, while semi-quantitative RT-PCR was used to assess extracellular matrix related gene expression. HVFF significantly proliferated (p = 0.011) when subjected to 200 Hz vibration and 20% strain, while BM-MSC did not (p = 1.0). A statistically significant increase in apoptosis of BM-MSC (p = 0.0402) was observed under the experimental conditions; however high cell viability (96%) was maintained. HVFF did not have significantly altered apoptosis (p = 0.7849) when subjected to vibration and strain. Semi-quantitative RT-PCR results show no significant differences in expression levels of collagen I (BM-MSC p = 0.1951, hVFF p = v0.3629), fibronectin (BM-MSC p = 0.1951, hVFF p = 0.2513), and TGF-β1 (BM-MSC p = 0.2534, hVFF p = 0.6029) between vibratory and static conditions in either cell type. Finally, smooth muscle actin mRNA was not present in either vibrated or static samples, indicating that no myofibroblast differentiation occurred for either cell type. Together, these results demonstrate that BM-MSC may be a suitable alternative to hVFF for vocal fold tissue engineering. Further investigation into a larger number of gene markers, protein levels, increased number of donors and vibratory conditions are warranted.     

Adapted to roar: functional morphology of tiger and lion vocal folds

Vocal production requires active control of the respiratory system, larynx and vocal tract. Vocal sounds in mammals are produced by flow-induced vocal fold oscillation, which requires vocal fold tissue that can sustain the mechanical stress during phonation. Our understanding of the relationship between morphology and vocal function of vocal folds is very limited. Here we tested the hypothesis that vocal fold morphology and viscoelastic properties allow a prediction of fundamental frequency range of sounds that can be produced, and minimal lung pressure necessary to initiate phonation. We tested the hypothesis in lions and tigers who are well-known for producing low frequency and very loud roaring sounds that expose vocal folds to large stresses. In histological sections, we found that the Panthera vocal fold lamina propria consists of a lateral region with adipocytes embedded in a network of collagen and elastin fibers and hyaluronan. There is also a medial region that contains only fibrous proteins and hyaluronan but no fat cells. Young's moduli range between 10 and 2000 kPa for strains up to 60%. Shear moduli ranged between 0.1 and 2 kPa and differed between layers. Biomechanical and morphological data were used to make predictions of fundamental frequency and subglottal pressure ranges. Such predictions agreed well with measurements from natural phonation and phonation of excised larynges, respectively. We assume that fat shapes Panthera vocal folds into an advantageous geometry for phonation and it protects vocal folds. Its primary function is probably not to increase vocal fold mass as suggested previously. The large square-shaped Panthera vocal fold eases phonation onset and thereby extends the dynamic range of the voice.     

Non-invasive in vivo measurement of the shear modulus of human vocal fold tissue

Voice is the essential part of singing and speech communication. Voice disorders significantly affect the quality of life. The viscoelastic mechanical properties of the vocal fold mucosa determine the characteristics of the vocal folds oscillations, and thereby voice quality. In the present study, a non-invasive method was developed to determine the shear modulus of human vocal fold tissue in vivo via measurements of the mucosal wave propagation speed during phonation. Images of four human subjects' vocal folds were captured using high speed digital imaging (HSDI) and magnetic resonance imaging (MRI) for different phonation pitches, specifically fundamental frequencies between 110 and 440 Hz. The MRI images were used to obtain the morphometric dimensions of each subject's vocal folds in order to determine the pixel size in the high-speed images. The mucosal wave propagation speed was determined for each subject and at each pitch value using an automated image processing algorithm. The transverse shear modulus of the vocal fold mucosa was then calculated from a surface (Rayleigh) wave propagation dispersion equation using the measured wave speeds. It was found that the mucosal wave propagation speed and therefore the shear modulus of the vocal fold tissue were generally greater at higher pitches. The results were in good agreement with those from other studies obtained via in vitro measurements, thereby supporting the validity of the proposed measurement method. This method offers the potential for in vivo clinical assessments of vocal folds viscoelasticity from HSDI.     

Stress wave velocities in bovine patellar tendon.

The velocity of longitudinal stress waves in an elastic body is given by the square root of the ratio of its elastic modulus to its density. In tendinous and ligamentous tissue, the elastic modulus increases with strain and with strain rate. Therefore, it was postulated that stress wave velocity would also increase with increasing strain and strain rate. The purpose of this study was to determine the velocity of stress waves in tendinous tissue as a function of strain and to compare these values to those predicted using the elastic modulus derived from quasi-static testing. Five bovine patellar tendons were harvested and potted as bone-tendon-bone specimens. Quasi-static mechanical properties were determined in tension at a deformation rate of 100 mm/s. Impact loading was employed to determine wave velocity at various strain levels, achieved by preloading the tendon. Following impact, there was a measurable delay in force transmission across the specimen and this delay decreased with increasing tendon strain. The wave velocities at tendon strains of 0.0075, 0.015, and 0.0225 were determined to be 260 +/- 52 m/s, 360 +/- 71 m/s, and 461 +/- 94 m/s, respectively. These velocities were significantly (p < 0.01) faster than those predicted using elastic moduli derived from the quasi-static tests by 52, 45, and 41 percent, respectively. This study has documented that stress wave velocity in patellar tendon increases with increasing strain and is underestimated with a modulus estimated from quasi-static testing.     

Mechanical stress during phonation in a self-oscillating finite-element vocal fold model

The stress information during phonation in the vocal folds is helpful in understanding the etiologies of vocal trauma and its related vocal diseases, such as nodules. In this paper, a self-oscillating finite-element model, which combines aerodynamic properties, tissue mechanics, airflow-tissue interactions, and vocal fold collisions, was used to simulate the vocal fold vibration during phonation. The spatial and temporal characteristics of mechanical stress in the vocal folds were predicted by this model. Temporally, it was found that mechanical stress periodically undulates with vibration of the vocal folds and that vocal fold impact causes a jump in the normal stress value. Spatially, the normal stress is significantly higher on the vocal fold surface than inside of the vocal folds. At the midpoint of the medial surface, the peak-to-peak amplitude of the normal stress reaches its maximum value. Using different lung pressures (0-1.5kPa) to drive the self-oscillating model, we found that lower lung pressure can effectively decrease the mechanical stress in the vocal folds. This study supports the fatigue damage hypothesis of vocal trauma. With this hypothesis and the numerical simulation in this study, the clinical observations of vocal fold trauma risk can be explained. This implies the mechanical stress predicted by this self-oscillating model could be valuable for predicting, preventing, and treating vocal fold injury.     

A computational study of systemic hydration in vocal fold collision

Mechanical stresses develop within vocal fold (VF) soft tissues due to phonation-associated vibration and collision. These stresses in turn affect the hydration of VF tissue and thus influence voice health. In this paper, high-fidelity numerical computations are described, taking into account fully 3D geometry, realistic tissue and air properties, and high-amplitude vibration and collision. A segregated solver approach is employed, using sophisticated commercial solvers for both the VF tissue and glottal airflow domains. The tissue viscoelastic properties were derived from a biphasic formulation. Two cases were considered, whereby the tissue viscoelastic properties corresponded to two different volume fractions of the fluid phase of the VF tissue. For each case, hydrostatic stresses occurring as a result of vibration and collision were investigated. Assuming the VF tissue to be poroelastic, interstitial fluid movement within VF tissue was estimated from the hydrostatic stress gradient. Computed measures of overall VF dynamics (peak airflow velocity, magnitude of VF deformation, frequency of vibration and contact pressure) were well within the range of experimentally observed values. The VF motion leading to mechanical stresses within the VFs and their effect on the interstitial fluid flux is detailed. It is found that average deformation and vibration of VFs tend to increase the state of hydration of the VF tissue, whereas VF collision works to reduce hydration.