Study of the Influence of Boundary Conditions on Corneal Deformation Based on the Finite Element Method of a Corneal Biomechanics Model

Abstract

Implementing in silico corneal biomechanical models for surgery applications can be boosted by developing patient-specific finite element models adapted to clinical requirements and optimized to reduce computational times. This research proposes a novel corneal multizone-based finite element model with octants and circumferential zones of clinical interest for material definition. The proposed model was applied to four patient-specific physiological geometries of keratoconus-affected corneas. Free-stress geometries were calculated by two iterative methods, the displacements and prestress methods, and the influence of two boundary conditions: embedded and pivoting. The results showed that the displacements, stress and strain fields differed for the stress-free geometry but were similar and strongly depended on the boundary conditions for the estimated physiological geometry when considering both iterative methods. The comparison between the embedded and pivoting boundary conditions showed bigger differences in the posterior limbus zone, which remained closer in the central zone. The computational calculation times for the stress-free geometries were evaluated. The results revealed that the computational time was prolonged with disease severity, and the displacements method was faster in all the analyzed cases. Computational times can be reduced with multicore parallel calculation, which offers the possibility of applying patient-specific finite element models in clinical applications.