Research Papers

The Role of Finite Displacements in Vocal Fold Modeling

[+] Author and Article Information
Haoxiang Luo

Department of Mechanical Engineering,
Vanderbilt University,
2301 Vanderbilt Place,
Nashville, TN 37235

James F. Doyle

School of Aeronautics and Astronautics,
Purdue University,
West Lafayette, IN 47907-2045

Bernard Rousseau

Departments of Otolaryngology and
Hearing and Speech Sciences,
Vanderbilt University,
1215 21st Avenue South,
Nashville, TN 37232

Contributed by the Bioengineering Division of ASME for publication in the Journal of Biomechanical Engineering. Manuscript received February 26, 2013; final manuscript received August 22, 2013; accepted manuscript posted September 6, 2013; published online October 1, 2013. Assoc. Editor: Ender A. Finol.

J Biomech Eng 135(11), 111008 (Oct 01, 2013) (8 pages) Paper No: BIO-13-1103; doi: 10.1115/1.4025330 History: Received February 26, 2013; Revised August 22, 2013; Accepted September 06, 2013

Human vocal folds experience flow-induced vibrations during phonation. In previous computational models, the vocal fold dynamics has been treated with linear elasticity theory in which both the strain and the displacement of the tissue are assumed to be infinitesimal (referred to as model I). The effect of the nonlinear strain, or geometric nonlinearity, caused by finite displacements is yet not clear. In this work, a two-dimensional model is used to study the effect of geometric nonlinearity (referred to as model II) on the vocal fold and the airflow. The result shows that even though the deformation is under 1 mm, i.e., less than 10% of the size of the vocal fold, the geometric nonlinear effect is still significant. Specifically, model I underpredicts the gap width, the flow rate, and the impact stress on the medial surfaces as compared to model II. The study further shows that the differences are caused by the contact mechanics and, more importantly, the fluid-structure interaction that magnifies the error from the small-displacement assumption. The results suggest that using the large-displacement formulation in a computational model would be more appropriate for accurate simulations of the vocal fold dynamics.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Ishizaka, K., and Flanagan, J. L., 1972, “Synthesis of Voiced Sounds From a Two-Mass Model of the Vocal Cords,” Bell Syst. Tech. J., 51, pp. 1233–1268. [CrossRef]
Jiang, J., and Zhang, Y., 2002, “Chaotic Vibration Induced by Turbulent Noise in a Two-Mass Model of Vocal Folds,” J. Acoust. Soc. Am., 112(5), pp. 2127–2133. [CrossRef] [PubMed]
Alipour, F., Brucker, C., Cook, D., Gommel, A., Kaltenbacher, M., Mattheus, W., Mongeau, L., Nauman, E., Schwarze, R., Tokuda, I., and Zorner, S., 2011, “Mathematical Models and Numerical Schemes for the Simulation of Human Phonation,” Cur. Bioinf., 6(3), pp. 323–343. [CrossRef]
Alipour, F., Berry, D., and Titze, I., 2000, “A Finite-Element Model of Vocal-Fold Vibration,” J. Acoust. Soc. Am., 14, pp. 442–454.
Luo, H., Mittal, R., and Bielamowicz, S. A., 2009, “Analysis of Flow-Structure Interaction in the Larynx During Phonation Using an Immersed-Boundary Method,” J. Acoust. Soc. Am., 126, pp. 816–824. [CrossRef] [PubMed]
Zheng, X., Xue, Q., Mittal, R., and Bielamowicz, S., 2010, “A Coupled Sharp-Interface Immersed-Boundary-Finiteelement Method for Flow-Structure Interaction With Application to Human Phonation,” ASME J. Biomech. Eng., 132(11), p. 111003. [CrossRef]
Titze, I., 1994, Principles of Voice Production, Prentice-Hall, Englewood Cliffs, NJ.
Pickup, B. A., and Thomson, S. L., 2011, “Identification of Geometric Parameters Influencing the Flow-Induced Vibration of a Two-Layer Self-Oscillating Computational Vocal Fold Model,” J. Acoust. Soc. Am., 129, pp. 2121–2132. [CrossRef] [PubMed]
Smith, S. L., and Thomson, S. L., 2012, “Effect of Inferior Surface Angle on the Self-Oscillation of a Computational Vocal Fold Model,” J. Acoust. Soc. Am., 131, pp. 4062–4075. [CrossRef] [PubMed]
Zhang, K., Siegmund, T., and Chan, R., 2006, “A Constitutive Model of the Human Vocal Fold Cover for Fundamental Frequency Regulation,” J. Acoust. Soc. Am., 119, pp. 1050–1062. [CrossRef] [PubMed]
Chan, R., Siegmund, T., and Zhang, K., 2009, “Biomechanics of Fundamental Frequency Regulation: Constitutive Modeling of the Vocal Fold Lamina Propria,” Logoped. Phoniatr. Vocol., 34(4), pp. 181–189. [CrossRef] [PubMed]
Doyle, J. F., 2001, Nonlinear Analysis of Thin-Walled Structures: Statics, Dynamics, and Stability, Springer-Verlag, New York.
Gunter, H., 2003, “A Mechanical Model of Vocal-Fold Collision With High Spatial and Temporal Resolution,” J. Acoust. Soc. Am., 113, pp. 994–1000. [CrossRef] [PubMed]
Tao, C., Jiang, J., and Zhang, Y., 2006, “Simulation of Vocal Fold Impact Pressures With a Self-Oscillating Finite-Element Model,” J. Acoust. Soc. Am., 119, pp. 3987–3994. [CrossRef] [PubMed]
Goldsmith, W., 2001, Impact: The Theory and Physical Behavior of Colliding Solids, Dover, New York.
Luo, H., Mittal, R., Zheng, X., Bielamowicz, S., Walsh, R., and Hahn, J., 2008, “An Immersed-Boundary Method for Flow–Structure Interaction in Biological Systems With Application to Phonation,” J. Comput. Phys., 227, pp. 9303–9332. [CrossRef] [PubMed]
Zemlin, W., 1997, Speech and Hearing Science: Anatomy & Physiology, 4th ed., Allyn and Bacon, Boston, MA.
Decker, G., and Thomson, S., 2007, “Computational Simulations of Vocal Fold Vibration: Bernoulli Versus Navier-Stokes,” J. Voice, 21(3), pp. 273–284. [CrossRef] [PubMed]
Titze, I. R., 1994, “Mechanical Stress in Phonation,” J. Voice, 8(2), pp. 99–105. [CrossRef] [PubMed]


Grahic Jump Location
Fig. 1

A 2D model for studying the role of geometric nonlinearity in the vocal fold dynamics

Grahic Jump Location
Fig. 2

The lowest four eigenmodes of the current vocal fold model. (a) The first mode at f1 = 86.9 Hz, (b) the second mode at f2 = 192.2 Hz, (c) the third mode at f3 = 217.5 Hz, and (d) the fourth mode at f4 = 372.9 Hz. The dashed line is for the rest shape, and the solid line is for deformed shape.

Grahic Jump Location
Fig. 3

The time history of the glottal gap width for model II with Psub = 0.6 kPa. (a) Initial transition period and (b) the established vibration.

Grahic Jump Location
Fig. 4

The vibration pattern of the vocal folds simulated by model I (a), (b), and (c) and model II (d), (e), and (f). The subglottal pressure is (a) and (d) Psub = 0.6 kPa, (b) and (e) Psub = 0.7 kPa, and (c) and (f) Psub = 0.8 kPa. The solid line and dashed line represent the open and closed phases, respectively. The circle shown in (a) is the marker point that will be used in the PSD analysis.

Grahic Jump Location
Fig. 5

The PSD analysis of the oscillation of a point at the medial surface for Psub =  0.7 and 0.8 kPa. The two dashed lines indicate the eigenfrequencies f2 and f3.

Grahic Jump Location
Fig. 6

(a) The gap width during the open phase versus Psub. (b) The waveform of the gap width at Psub = 0.6 kPa. In both (a) and (b), the dashed line is for model I and the solid line for model II.

Grahic Jump Location
Fig. 7

The impact stress (a) and the total contact force (per unit length in span) (b) during vocal fold collision. The dashed line is for model I and the solid line for model II.

Grahic Jump Location
Fig. 8

(a) The time-averaged flow rate (per unit span) and (b) the waveform of the flow rate at Psub = 0.8 kPa. The dashed line is for model I and the solid line for model II.

Grahic Jump Location
Fig. 10

(a) Illustration of the load tests of the current vocal fold model. The reference point marked by a circle is used to measure the displacement. (b) and (c) p/Psub from the FSI simulation during vocal fold opening for (b) Psub = 0.6 kPa and (c) Psub = 0.8 kPa.

Grahic Jump Location
Fig. 9

Instantaneous vorticity contours during vocal fold opening for (a) model I and (b) model II at Psub = 0.8 kPa

Grahic Jump Location
Fig. 11

Displacement of the reference point in the static load test, the dynamic load test, and the coupled FSI simulation. (a) The total displacement and (b) the y-displacement.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In