Automatized Patient-Specific Methodology for Numerical Determination of Biomechanical Corneal Response

This work presents a novel methodology for building a three-dimensional patient-specific eyeball model suitable for performing a fully automatic finite element (FE) analysis of the corneal biomechanics. The reconstruction algorithm fits and smooths the patient’s corneal surfaces obtained in clinic with corneal topographers and creates an FE mesh for the simulation. The patient’s corneal elevation and pachymetry data is kept where available, to account for all corneal geometric features (central corneal thickness–CCT and curvature). Subsequently, an iterative free-stress algorithm including a fiber’s pull-back is applied to incorporate the pre-stress field to the model. A convergence analysis of the mesh and a sensitivity analysis of the parameters involved in the numerical response is also addressed to determine the most influential features of the FE model. As a final step, the methodology is applied on the simulation of a general non-commercial non-contact tonometry diagnostic test over a large set of 130 patients—53 healthy, 63 keratoconic (KTC) and 14 post-LASIK surgery eyes. Results show the influence of the CCT, intraocular pressure (IOP) and fibers (87%) on the numerical corneal displacement (Formula presented.) the good agreement of the (Formula presented.) with clinical results, and the importance of considering the corneal pre-stress in the FE analysis. The potential and flexibility of the methodology can help improve understanding of the eye biomechanics, to help to plan surgeries, or to interpret the results of new diagnosis tools (i.e., non-contact tonometers).