The mechanical response of the cornea subjected to a non-contact air-jet tonometry diagnostic test represents an interplay between its geometry, the corneal material behavior and the loading. The objective is to study this interplay to better understand and interpret the results obtained with a non-contact tonometry test. A patient-specific finite element model of a healthy eye, accounting for the load free configuration, was used. The corneal tissue was modeled as an anisotropic hyperelastic material with two preferential directions. Three different sets of parameters within the human experimental range obtained from inflation tests were considered. The influence of the IOP was studied by considering four pressure levels (10–28 mmHg) whereas the influence of corneal thickness was studied by inducing a uniform variation (300–600 microns). A Computer Fluid Dynamics (CFD) air-jet simulation determined pressure loading exerted on the anterior corneal surface. The maximum apex displacement showed a linear variation with IOP for all materials examined. On the contrary, the maximum apex displacement followed a cubic relation with corneal thickness. In addition, a significant sensitivity of the apical displacement to the corneal stiffness was also obtained. Explanation to this behavior was found in the fact that the cornea experiences bending when subjected to an air-puff loading, causing the anterior surface to work in compression whereas the posterior surface works in tension. Hence, collagen fibers located at the anterior surface do not contribute to load bearing. Non-contact tonometry devices give useful information that could be misleading since the corneal deformation is the result of the interaction between the mechanical properties, IOP, and geometry. Therefore, a non-contact tonometry test is not sufficient to evaluate their individual contribution and a complete in-vivo characterization would require more than one test to independently determine the membrane and bending corneal behavior.
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 (U(Num)) the good agreement of the U(Num) 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).
Keratoconus is an idiopathic, non-inflammatory and degenerative corneal disease characterised by a loss of the organisation in the corneal collagen fibrils. As a result, keratoconic corneas present a localised thinning and conical protrusion with irregular astigmatism and high myopia that worsen visual acuity. Intracorneal ring segments (ICRSs) are used in clinic to regularise the corneal surface and to prevent the disease from progressing. Unfortunately, the post-surgical effect of the ICRS is not explicitly accounted beforehand. Traditional treatments rely on population-based nomograms and the experience of the surgeon. In this vein, in silico models could be a clinical aid tool for clinicians to plan the intervention, or to test the post-surgical impact of different clinical scenarios. A semi-automatic computational methodology is presented in order to simulate the ICRS surgical operation and to predict the post-surgical optical outcomes. For the sake of simplicity, circular cross section rings, average corneas and an isotropic hyperelastic material are used. To determine whether the model behaves physiologically and to carry out a sensitivity analysis, a [Formula: see text] full-factorial analysis is carried out. In particular, how the stromal depth insertion, horizontal distance of ring insertion (hDRI) and diameter of the ring's cross section ([Formula: see text]) are impacting in the spherical and cylindrical power of the cornea is analysed. Afterwards, the kinematics, mechanics and optics of keratoconic corneas after the ICRS insertion are analysed. Based on the parametric study, we can conclude that our model follows clinical trends previously reported. In particular and although there is an improvement in defocus, all corneas presented a change in their optical aberrations. The stromal depth insertion is the parameter that affects the corneal optics the most, whereas hDRI and [Formula: see text] are less important. Not only that, but it is almost impossible to achieve an optimal trade-off between spherical and cylindrical correction. Regarding the mechanical behaviour, inserting the rings at 65% depth or above will cause the cornea to slightly bend. This abnormal stress distribution greatly distorts the corneal optics and, more importantly, could be the cause of clinical problems such as corneal extrusion. Not only that, but our model also supports that rings are acting as restraint elements which relax the stresses of the corneal stroma in the cone of the disease. However, depending on the exact spatial location of the keratoconus, the insertion of rings could promote its evolution instead of preventing it. ICRS inserted deeper will prevent keratoconus in the posterior stroma from growing (relaxation of posterior surface), but will promote its growing if they are located in the anterior surface (increment of stress). In conclusion, the methodology proposed is suitable for simulating long-term mechanical and optical effects of ICRS insertion.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.