A simple-shear rheometer for linear viscoelastic characterization of vocal fold tissues at phonatory frequencies

Previous studies reporting the linear viscoelastic shear properties of the human vocal fold cover or mucosa have been based on torsional rheometry, with measurements limited to low audio frequencies, up to around 80 Hz. This paper describes the design and validation of a custom-built, controlled-strain, linear, simple-shear rheometer system capable of direct empirical measurements of viscoelastic shear properties at phonatory frequencies. A tissue specimen was subjected to simple shear between two parallel, rigid acrylic plates, with a linear motor creating a translational sinusoidal displacement of the specimen via the upper plate, and the lower plate transmitting the harmonic shear force resulting from the viscoelastic response of the specimen. The displacement of the specimen was measured by a linear variable differential transformer whereas the shear force was detected by a piezoelectric transducer. The frequency response characteristics of these system components were assessed by vibration experiments with accelerometers. Measurements of the viscoelastic shear moduli (G' and G") of a standard ANSI S2.21 polyurethane material and those of human vocal fold cover specimens were made, along with estimation of the system signal and noise levels. Preliminary results showed that the rheometer can provide valid and reliable rheometric data of vocal fold lamina propria specimens at frequencies of up to around 250 Hz, well into the phonatory range.     

Resonance properties of the vocal folds: in vivo laryngoscopic investigation of the externally excited laryngeal vibrations

The study presents the first attempt to investigate resonance properties of the living vocal folds by means of laryngoscopy. Laryngeal vibrations were excited via a shaker placed on the neck of a male subject and observed by means of videostroboscopy and videokymography (VKG). When the vocal folds were tuned to the phonation frequency of 110 Hz and sinusoidal vibration with sweeping frequency (in the range 50-400 Hz) was delivered to the larynx, three clearly pronounced resonance peaks at frequencies around 110, 170, and 240 Hz were identified in the vocal fold tissues. Different modes of vibration of the vocal folds, observed as distinct lateral-medial oscillations with one, two, and three half-wavelengths along the glottal length, respectively, were associated with these resonance frequencies. At the external excitation frequencies below 100 Hz, vibrations of the ventricular folds, aryepiglottic folds and arytenoid cartilages were dominant in the larynx.     

Simulation based estimation of dynamic mechanical properties for viscoelastic materials used for vocal fold models

In order to obtain a deeper understanding of the human phonation process and the mechanisms generating sound, realistic setups are built up containing artificial vocal folds. Usually, these vocal folds consist of viscoelastic materials (e.g., polyurethane mixtures). Reliable simulation based studies on the setups require the mechanical properties of the utilized viscoelastic materials. The aim of this work is the identification of mechanical material parameters (Young's modulus, Poisson's ratio, and loss factor) for those materials. Therefore, we suggest a low-cost measurement setup, the so-called vibration transmission analyzer (VTA) enabling to analyze the transfer behavior of viscoelastic materials for propagating mechanical waves. With the aid of a mathematical Inverse Method, the material parameters are adjusted in a convenient way so that the simulation results coincide with the measurement results for the transfer behavior. Contrary to other works, we determine frequency dependent functions for the mechanical properties characterizing the viscoelastic material in the frequency range of human speech (100–250 Hz). The results for three different materials clearly show that the Poisson's ratio is close to 0.5 and that the Young's modulus increases with higher frequencies. For a frequency of 400 Hz, the Young's modulus of the investigated viscoelastic materials is approximately 80% higher than for the static case (0 Hz). We verify the identified mechanical properties with experiments on fabricated vocal fold models. Thereby, only small deviations between measurements and simulations occur.     

Optoreflectometry Determination of the Resonance Properties of a Vocal Fold

A new method of measuring the resonance properties of a vocal fold using electromagnetic excitation and laser optoreflectometry for response monitoring is described. Two resonance peaks were experimentally identified with one magnet stuck on the vocal fold at frequencies F0(1m)=54.7 Hz and F0'(1m)=35.8 Hz. The addition of a second magnet allowed calculation of the actual viscoelastic properties of the vocal fold: F0=71.8 Hz; quality factor Q=8.03; mass m=0.057 g; stiffness k=11.6 Nm; and damping zeta=0.0032 Nm(-1). A numerical simulation of a two-layered model verified the experimental data.     

The shear modulus of the human vocal fold in a transverse direction

The aim of this study was to measure the shear modulus of the vocal fold in a human hemilarynx, such that the data can be related to direction of applied stress and anatomical context. Dynamic spring rate data were collected using a modified linear skin rheometer using human hemilarynges, and converted to estimated shear modulus via application of a simple shear model. The measurement probe was attached to the epithelial layer of the vocal fold cover using suction. A sinusoidal force of 3g was applied to the epithelium, and the resultant displacement logged at a rate of 1kHz. Force measurement accuracy was 20microg and position measurement accuracy was 4microm. The force was applied in a transverse direction at the midmembranous point between the vocal process and the anterior commissure. The shear modulus of the three female vocal folds ranged from 814 to 1232Pa. The shear modulus of the three male vocal folds ranged from 1021 to 1796Pa. These data demonstrate that it is possible to obtain estimates for the shear modulus of the vocal fold while preserving anatomical context. The modulus values reported here are higher than those reported using parallel plate rheometry. This is to be expected as the tissue is attached to surrounding structures, and is under natural tension.     

Histologic and Rheologic Characterization of Vocal Fold Scarring

Scarring of the vocal fold causes considerable dysphonia and presents significant treatment challenges. A rabbit model was developed to investigate the histologic ultrastructure and rheologic properties of the scarred vocal fold lamina propria. Eleven rabbit larynges were scarred by means of forcep biopsy. Sixty days postoperatively, the rabbits were sacrificed and their vocal folds were harvested. Histological analysis of the scarred and normal lamina propria was completed for collagen, procollagen, elastin, and hyaluronic acid. Linear viscoelastic shear properties of the tissues were also measured, including elastic shear modulus and dynamic viscosity. Compared to normal vocal fold lamina propria, scarred tissues demonstrated significantly less collagen, an increase in procollagen, and a decrease in elastin. Rheologically, both elastic shear modulus and dynamic viscosity were significantly higher for the scarred tissues. Increased stiffness and viscosity do not appear to result from an increase in collagen, but rather appear to be related to the presence of new, disorganized collagen scaffolding. Results are interpreted in terms of the possible role of interstitial proteins in the etiology of increased stiffness and viscosity, which requires further investigation. This animal model should allow for systematic future investigations of vocal fold scarring and its treatment.     

Slow nonlinearity kinetics of the viscoelastic properties of oil during shear vibrations

Earlier we discovered the slow evolution of viscoelastic moduli of heavy crude oil. The shear modulus was measured at frequencies of 0.5, 5, and 50 Hz at different temperatures over 72 h. New studies of the dependence of the complex shear modulus on the strain amplitude revealed a logarithmic increase in the nonlinearity parameter as a function of time for this oil sample. It was experimentally established that the complex shear modulus is a linear function of the amplitude of shear perturbations. This is possible in the case of a linear dependence of values of the viscoelastic characteristics on the medium modulus of deformation.     

Dependence of phonation threshold pressure on vocal tract acoustics and vocal fold tissue mechanics

Analytical and computer simulation studies have shown that the acoustic impedance of the vocal tract as well as the viscoelastic properties of vocal fold tissues are critical for determining the dynamics and the energy transfer mechanism of vocal fold oscillation. In the present study, a linear, small-amplitude oscillation theory was revised by taking into account the propagation of a mucosal wave and the inertive reactance (inertance) of the supraglottal vocal tract as the major energy transfer mechanisms for flow-induced self-oscillation of the vocal fold. Specifically, analytical results predicted that phonation threshold pressure (Pth) increases with the viscous shear properties of the vocal fold, but decreases with vocal tract inertance. This theory was empirically tested using a physical model of the larynx, where biological materials (fat, hyaluronic acid, and fibronectin) were implanted into the vocal fold cover to investigate the effect of vocal fold tissue viscoelasticity on Pth. A uniform-tube supraglottal vocal tract was also introduced to examine the effect of vocal tract inertance on Pth. Results showed that Pth decreased with the inertive impedance of the vocal tract and increased with the viscous shear modulus (G") or dynamic viscosity (eta') of the vocal fold cover, consistent with theoretical predictions. These findings supported the potential biomechanical benefits of hyaluronic acid as a surgical bioimplant for repairing voice disorders involving the superficial layer of the lamina propria, such as scarring, sulcus vocalis, atrophy, and Reinke's edema.     

Characterization of chronic vocal fold scarring in a rabbit model

The purpose of the current study was to assess the histologic and rheologic properties of the scarred vocal fold lamina propria during a chronic phase of wound repair in a rabbit model. Eighteen rabbit larynges were scarred using a procedure that involved stripping the vocal fold lamina propria down to the thyroarytenoid muscle, using 3-mm microforceps. The approximate dimension of injury to the vocal fold was 3 x 1.5 x 0.5 mm [length x width x depth]. At 6 months postoperatively, histologic analysis of the scarred and control lamina propria in eight of these rabbits was completed for collagen, procollagen, elastin, and hyaluronic acid. Compared with control samples, scarred tissue samples revealed fragmented and disorganized elastin fibers. Additionally, collagen was significantly increased, organized, and formed thick bundles in the scarred vocal fold lamina propria. Measurements of the viscoelastic shear properties of the scarred and control lamina propria in the remaining 10 rabbits revealed increased elastic shear modulus (G') in 8 of 10 scarred samples and increased dynamic viscosity (eta') in 9 of 10 scarred samples. Although rheologic differences were not statistically significant, they revealed that on average, scarred samples were stiffer and more viscous than the normal controls. Histologic data are interpreted as indicating that by 6 months postinjury, the scarred rabbit vocal fold has reached a mature phase of wound repair, characterized by an increased, organized, and thick bundle collagen matrix. Rheologic data are interpreted as providing support for the potential role of increased, thick bundle collagen, and a disorganized elastin network on shear stiffness and dynamic viscosity in the chronic vocal fold scar. Based on these results, a 6-month postoperative time frame is proposed for future studies of chronic vocal fold scarring using the rabbit animal model.     

Shear waves in viscoelastic wormlike micellar fluids

Low frequency (61 Hz) shear wave speeds have been measured in viscoelastic wormlike micellar (WM) fluids for a concentration range of 20/12-160/96 mM CTAB/NaSAL. The strain induced birefringence of the WM fluids was exploited to optically track the shear pulse using crossed polarizing filters and high speed video. It was found that shear speed increases roughly linearly with concentration at a rate of 3.5 mm s(-1) mM(-1) CTAB. Further, comparison with elastic and loss moduli obtained from rheology data show that shear wave propagation is dominated by the elastic modulus for this frequency range.     

Transient elastography in anisotropic medium: application to the measurement of slow and fast shear wave speeds in muscles

From the measurement of a low frequency (50-150 Hz) shear wave speed, transient elastography evaluates the Young's modulus in isotropic soft tissues. In this paper, it is shown that a rod source can generate a low frequency polarized shear strain waves. Consequently this technique allows to study anisotropic medium such as muscle. The evidence of the polarization of low frequency shear strain waves is supported by both numeric simulations and experiments. The numeric simulations are based on theoretical Green's functions in isotropic and anisotropic media (hexagonal system). The experiments in vitro led on beef muscle proves the pertinent of this simple anisotropic pattern. Results in vivo on man biceps shows the existence of slow and fast shear waves as predicted by theory.     

Resonance method of measuring shear viscoelastic properties of liquid media based on excitation of torsional oscillations in tubes

Possibilities of using torsional oscillations for measuring viscoelastic properties of liquids are discussed. The theory of torsional oscillations of an elastic tube filled with the media to be investigated possessing viscosity and shear elasticity is developed. It is shown that to determine a complex shear modulus it is sufficient to determine the resonance frequency and Q-factor of torsional oscillations. An experimental installation and the results of measurements of viscoelastic modulus of glycerin and oil of one oilfield within the temperature range from −10° to 60°C are given. The experimental installation allows measuring a viscoelastic modulus within the range of acoustic logging frequencies (10–20 kHz). The obtained results are compared with the results of rheometric measurements.     

Methodology for rheological testing of engineered biomaterials at low audio frequencies

A commercial rheometer (Bohlin CVO120) was used to mechanically test materials that approximate vocal-fold tissues. Application is to frequencies in the low audio range (20-150 Hz). Because commercial rheometers are not specifically designed for this frequency range, a primary problem is maintaining accuracy up to (and beyond) the mechanical resonance frequency of the rotating shaft assembly. A standard viscoelastic material (NIST SRM 2490) has been used to calibrate the rheometric system for an expanded frequency range. Mathematically predicted response curves are compared to measured response curves, and an error analysis is conducted to determine the accuracy to which the elastic modulus and the shear modulus can be determined in the 20-150-Hz region. Results indicate that the inertia of the rotating assembly and the gap between the plates need to be known (or determined empirically) to a high precision when the measurement frequency exceeds the resonant frequency. In addition, a phase correction is needed to account for the magnetic inertia (inductance) of the drag cup motor. Uncorrected, the measured phase can go below the theoretical limit of -pi. This can produce large errors in the viscous modulus near and above the resonance frequency. With appropriate inertia and phase corrections, +/- 10% accuracy can be obtained up to twice the resonance frequency.     

Viscoelastic modeling of canine vocalis muscle in relaxation

Analysis of vocal fold vibration requires information on the viscoelastic properties of the vocalis muscle. The force response of two canine vocalis muscles was measured in one-dimensional, stepwise elongation of the tissue as a function of time with a computer-controlled ergometer. The viscoelastic behavior of the muscle in its passive state was demonstrated through the relaxation test. A quasilinear viscoelastic model was used to parametrize the relaxation function, and results are reported for various levels of strain. Furthermore, a model was used to obtain theoretical time-dependent stress-strain curves to compare with experimental data.     

Visualization and Quantification of the Medial Surface Dynamics of an Excised Human Vocal Fold During Phonation

SUMMARY: The purpose of this investigation was to investigate physical mechanisms of vocal fold vibration during normal phonation through quantification of the medial surface dynamics of the fold. An excised hemilarynx setup was used. The dynamics of 30 microsutures mounted on the medial surface of a human vocal fold were analyzed across 18 phonatory conditions. The vibrations were recorded with a digital high-speed camera at a frequency of 4,000 Hz. The positions of the sutures were extracted and converted to three-dimensional coordinates using a linear approximation technique. The data were reduced to principal eigenfuctions, which captured over 90% of the variance of the data, and suggested mechanisms of sustained vocal fold oscillation. The vibrations were imaged as the following phonatory conditions were manipulated: glottal airflow, an adductory force applied to the muscular process, and an elongation force applied to the thyroid cartilage. Over the range of variables studied, only the variation in glottal airflow yielded significant changes in subglottal pressure and fundamental frequency. All recordings showed high correlation for the distribution of the dynamics across the medial surface of the vocal fold. The distribution of the different displacement directions and velocities showed the highest variations around the superior region of the medial surface. Although the computed vibration patterns of the two largest empirical eigenfunctions were consistent with previous experimental observations, the relative prominence of the two eigenfunctions changed as a function of glottal airflow, impacting theories of vocal efficiency and vocal economy.     

On the relation between subglottal pressure and fundamental frequency in phonation

The change in fundamental frequency with subglottal pressure in phonation is quantified on the basis of the ratio between vibrational amplitude and vocal fold length. This ratio is typically very small in stringed instruments, but becomes quite appreciable in vocal fold vibration. Tension in vocal fold tissues is, therefore, not constant over the vibratory cycle, and a dynamic tension gives rise to amplitude-frequency dependence. It is shown that the typical 2-6 Hz/cm H2O rise in fundamental frequency with subglottal pressure observed in human and canine larynges is a direct and predictable consequence of this amplitude-frequency dependence. Results are presently limited to phonation in the chest register.     

Mechanical behavior of the smectic phases near the phase transitions

We present experimental data on the apparent layer complex rigidity of the smectic A, C, B phases of 40 P90B versus temperature at a fixed frequency of 50 Hz. Across the nematic → smectic A and smectic A → smectic C second order transitions, the modulus of the rigidity coefficient becomes very weak, and its phase undergoes a rapid variation of π, in a range of about 100 mK. At the first order smectic C → smectic $\text{B}$ transition, a corresponding discontinuity is observed. This behavior is tentatively explained in term of coupling of the nematic director (or defects) with the shear flow associated with the compression, and an associated viscoelastic relaxation.     

Using the Relaxation Oscillations Principle for Simple Phonation Modeling

Well-known multimass models of vocal folds are useful to describe main behavior observed in human voicing but their principle of functioning, based on harmonic oscillation, may appear complex. This work is designed to show that a simple one-mass model ruled by laws of relaxation oscillation can also depict main behavior of glottis dynamic. Theory of relaxation oscillation is detailed. A relaxation oscillation model is assessed through a numerical simulation using conventional values for tissue characteristics and subglottal pressure. As expected, raising the mass decreases the fundamental frequency and increases the amplitude of vocal fold vibration: for a mass ranging from 0.01 to 0.4 g, F0 decreased from 297.5 to 42.5 Hz and vibrational amplitude increased from 1.26 to 3.25 mm (for stiffness k=10Nm(-1), damping r=0.015 N s m(-1), and subglottal pressure=1 kPa). Stiffness value has the opposite effect. The subglottal pressure controls the fundamental frequency with a rate ranging from 20 to 50 Hz/kPa. The vibrational amplitude is also controlled linearly by subglottal pressure from 0.22 to 0.26 mm/kPa. The range of phonation threshold pressure (PTP) is close to the values currently proposed, that is, 0.1 to 1 kPa and varies with the fundamental frequency. The relaxation oscillator is a simple and useful tool for modeling vocal fold vibration.