A novel numerical modelling approach for keratoplasty eye procedure.

Objective of the work is to investigate stress and deformation that conrneal tissue and donor graft undergo during endothelial keratoplasty. In order to attach the donor graft to the cornea, different air bubble pressure profiles acting on the graft are considered. This study is carried out by employing a three-dimensional nonlinear finite element methodology, combined with a contact algorithm. The ocular tissues are treated as isotropic, hyper-elastic and nearly-incompressible materials. The contact algorithm, based on the penalty-based node-to-surface approach, is used to model the donor graft-corneal interface region. First, the proposed computational methodology is tested against benchmark data for bending of the plates over a cylinder. Then, the influence of geometrical and material parameters of the graft on the corneal contact-structural response is investigated. The results are presented in terms of Von Mises stress intensity, displacement and mean contact force. Results clearly indicate that the air bubble pressure plays a key role in the corneal stress and strain, as well as graft stiffness and thickness.