May 2008
Volume 49, Issue 13
ARVO Annual Meeting Abstract  |   May 2008
Finite Element Model of the Human Cornea for Simulation of Intrastromal Laser Refractive Surgery
Author Affiliations & Notes
  • J. F. Bille
    Physics, University of Heidelberg, Mannheim, Germany
  • H. Zhang
    Physics, University of Heidelberg, Mannheim, Germany
  • Footnotes
    Commercial Relationships  J.F. Bille, None; H. Zhang, None.
  • Footnotes
    Support  None.
Investigative Ophthalmology & Visual Science May 2008, Vol.49, 668. doi:
  • Views
  • Share
  • Tools
    • Alerts
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      J. F. Bille, H. Zhang; Finite Element Model of the Human Cornea for Simulation of Intrastromal Laser Refractive Surgery. Invest. Ophthalmol. Vis. Sci. 2008;49(13):668.

      Download citation file:

      © ARVO (1962-2015); The Authors (2016-present)

  • Supplements

Purpose: : To develop a three-dimensional, non-linear, non-homogeneous finite element model of the human cornea to estimate three-dimensional stresses and finite strains in the cornea.

Methods: : Finite element modeling method is used to study the biomechanical properties of the cornea. Three-dimensional spherical polar coordinates are used to avoid membrane or shell theories and minimize the number of degrees-of-freedom needed to obtain converged solutions. 144 three-dimensional elements are used in the current model and each element is associated with 8 nodes. The cornea is assumed to be a slightly compressible, non-linear, non-homogenous material. The whole cornea is divided into 9 layers (3 layers for Bowman’s layer and 6 layers for stroma) and Bowman’s layer is assumed to be stiffer than stroma. Intraocular pressure is applied homogeneously to the posterior layer of the cornea during the simulation. A new boundary condition is proposed and applied in the current model that only Bowman’s layer is fully fixed instead of the whole cornea at the limbus.

Results: : Intraocular pressure (IOP) is simulated to be increasingly applied to the cornea from 0 kPa (stress-free state) to 4 kPa and the corresponding apical rise is calculated by finite element analysis. The "Apical-rise vs. IOP" curve shows a good agreement with the published experimental data. Three-dimensional stresses and strains at the center of each element are calculated under normal IOP (2 kPa) and abnormal IOP (4 kPa). Three locations (apex, mid-section and near-limbus) in the cornea are chosen for further quantitative evaluation. 3-D stresses and strains within Bowman’s layer and stroma are also evaluated. The boundary conditions at the limbus are varied to study how such clamping conditions affect the stress and strain estimates.

Conclusions: : A new three-dimensional, non-linear, heterogeneous finite element model of the human cornea has been established to study the biomechanics of the cornea, which is important for choosing proper surgical parameters and predicting post-operative outcomes.

Keywords: computational modeling • cornea: clinical science • cornea: stroma and keratocytes 

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.