Multiscale mechanical model based on patient-specific geometry: application to early keratoconus development

Abstract: Keratoconus is a pathology of the cornea associated with a tissue thinning and a weakening of its mechanical properties. However, it remains elusive which aspect is the leading cause of the disease. To investigate this question, we combined a multiscale model with a patient-geometry in order to simulate the mechanical response of healthy and pathological corneas under intraocular pressure. The constitutive behavior of the cornea is described through an energy function which takes into account the isotropic matrix of the cornea, the geometric structure of collagen lamellae and the quasi-incompressibility of the tissue. A micro-sphere description is implemented to take into account the typical features of the collagen lamellae as obtained experimentally, namely their orientation, their stiffness and their dispersion, as well as the their unfolding stretch, at which they start to provide a significant force. A set of reference parameters is obtained to fit experimental inflation data of the literature. We show that the most sensitive parameter is the unfolding stretch, as a small variation of this parameter induces a major change in the corneal apex displacement. The keratoconus case is then studied by separating the impact of the geometry and the one of the mechanics. We computed the evolution of the SimK (a clinical indicator of cornea curvature) and elevation maps: we were able to reproduce the reported changes of SimK with pressure only by a mechanical weakening, and not by a change in geomtry. More specifically, the weakening has to target the lamellae and not the matrix. The mechanical weakening leads to elevations close to early stage keratoconus, but our model lacks the remodeling component to couple the change in mechanics with changes in geometry.