Biomechanical modeling of femtosecond laser keyhole endokeratophakia surgery

PURPOSE: To apply a finite element model to endokeratophakia and evaluate anterior and posterior corneal surface changes. METHODS: Spatial elevation data (Pentacam HR; Oculus, Wetzlar, Germany) were obtained for the front and back corneal surfaces of an eye prior to undergoing an endokeratophakia procedure. These were used to warp a spherical template finite element model of the cornea to create a patient-specific finite element mesh and the initial stress distribution was computed with an iterative approach. The finite element model (Optimeyes; Integrated Scientific Services, Biel, Switzerland) included non-linear elastic characteristics of the stroma. The endokeratophakia procedure was recreated in the model: A donor lenticule (-10.50 diopters [D], 5.75-mm zone, 127-μm thick) was inserted into a lamellar pocket (180-μm deep, 6.25-mm diameter) and two 2-mm small incisions were made at 150° and 330°. Anterior and posterior surfaces, computed by the finite element model, were compared to clinical data to assess accuracy and reliability of finite element modeling. RESULTS: The postoperative axial curvature produced by the model closely resembled the patient data; average curvature was 48.01 D clinically and 48.23 D in the simulation, and corneal astigmatism was 3.01 D clinically and 2.88 D in the simulation. The posterior best-fit sphere elevation map also matched the patient data, replicating inward bulging of the posterior surface by approximately 40 μm. Stress distribution modeling predicted a stress increase by 159.94% ± 73% in the cap and a stress decrease by 32.41% ± 21% in the stromal bed. CONCLUSIONS: Finite element modeling of the cornea reproduced the clinically observed anterior and posterior corneal surface changes following an endokeratophakia procedure. This case sets the stage for further study to refine and yield predictive finite element modeling for the evaluation of corneal refractive surgical procedures.