Abstract:Purpose: To improve the effectivity of patient-specific finite element analysis (FEA) to predict refractive power change (RPC) in rigid Ortho-K contact lens fitting. Novel eyelid boundary detection is introduced to the FEA model to better model the effects of the lid on lens performance, and stress and strain outcomes are investigated to identify the most effective FEA components to use in modelling. Methods: The current study utilises fully anonymised records of 249 eyes, 132 right eyes, and 117 left eyes from subjects aged 14.1 ± 4.0 years on average (range 9 to 38 years), which were selected for secondary analysis processing. A set of custom-built MATLAB codes was built to automate the process from reading Medmont E300 height and distance files to processing and displaying FEA stress and strain outcomes. Measurements from before and after contact lens wear were handled to obtain the corneal surface change in shape and power. Tangential refractive power maps were constructed from which changes in refractive power pre- and post-Ortho-K wear were determined as the refractive power change (RPC). A total of 249 patient-specific FEA with innovative eyelid boundary detection and 3D construction analyses were automatically built and run for every anterior eye and lens combination while the lens was located in its clinically detected position. Maps of four stress components: contact pressure, Mises stress, pressure, and maximum principal stress were created in addition to maximum principal logarithmic strain maps. Stress and strain components were compared to the clinical RPC maps using the two-dimensional (2D) normalised cross-correlation and structural similarity (SSIM) index measure. Results: On the one hand, the maximum principal logarithmic strain recorded the highest moderate 2D cross-correlation area of 8.6 ± 10.3%, and contact pressure recorded the lowest area of 6.6 ± 9%. Mises stress recorded the second highest moderate 2D cross-correlation area with 8.3 ± 10.4%. On the other hand, when the SSIM index was used to compare the areas that were most similar to the clinical RPC, maximum principal stress was the most similar, with an average strong similarity percentage area of 26.5 ± 3.3%, and contact pressure was the least strong similarity area of 10.3 ± 7.3%. Regarding the moderate similarity areas, all components were recorded at around 34.4% similarity area except the contact pressure, which was down to 32.7 ± 5.8%. Conclusions: FEA is an increasingly effective tool in being able to predict the refractive outcome of Ortho-K treatment. Its accuracy depends on identifying which clinical and modelling metrics contribute to the most accurate prediction of RPC with minimal ocular complications. In terms of clinical metrics, age, Intra-ocular pressure (IOP), central corneal thickness (CCT), surface topography, lens decentration and the 3D eyelid effect are all important for effective modelling. In terms of FEA components, maximum principal stress was found to be the best FEA barometer that can be used to predict the performance of Ortho-K lenses. In contrast, contact pressure provided the worst stress performance. In terms of strain, the maximum principal logarithmic strain was an effective strain barometer.

Finite Element Analysis of LASIK

The finite element analysis of the pathology and treatment of high myopia esotropia

Analytic Modeling and Simulating of the Cornea with Finite Element Method

Effects of Laser In Situ Keratomileusis and Small-Incision Lenticule Extraction on Corneal Biomechanical Behavior: A Finite Element Analysis

Biomechanical Simulation of Stress Concentration and Intraocular Pressure in Corneas Subjected to Myopic Refractive Surgical Procedures

Small Incision Lenticule Extraction and Femtosecond-Assisted Laser in Situ Keratomileusis in Patients with Deep Corneal Opacity: Case Series.

Biomechanical model of cornea deformation and simulation

Theoretical Analysis of Wave-Front Aberrations Induced from Conventional Laser Refractive Surgery in a Biomechanical Finite Element Model

Comparison of corneaI biomechanicaI changes at earIy stage after femtosecond LASIK and microkeratome LASIK

Comparison on biomechanical stability between femtosecond laser assisted LASIK and LASEK with ocular response analyzer

Characteristics and Correlative Factors of Posterior Corneal Surface Changes after Laser in Situ Keratomileusis

FEA-Based Stress-Strain Barometers as Forecasters for Corneal Refractive Power Change in Orthokeratology

Research of Corneal Biomechanical Properties after LASIK Surgery

Simulated biomechanical effect of aspheric transition zone ablation profiles after conventional hyperopia refractive surgery

[Advances in Biomechanical Effects of Laser Corneal Refractive Surgery].

Three-dimensional biomechanical modeling and simulation of trephine cutting cornea for keratoplasty.

A review of human cornea finite element modeling: geometry modeling, constitutive modeling, and outlooks

Customized eye modeling for optical quality assessment in myopic femto-LASIK surgery

Biomechanical and refractive behaviors of keratoconic cornea based on three-dimensional anisotropic hyperelastic models.

The Biomechanical Response of the Cornea in Orthokeratology

Quantitative assessment of corneal elasticity distribution after FS‐LASIK using optical coherence elastography