Biomechanical model of the human cornea: considering shear stiffness and regional variation of collagen anisotropy and density

A numerical model based on continuum mechanics theory has been developed which represents the 3D anisotropic behaviour of the corneal stroma. Experimental data has been gathered from a number of previous studies to provide the basis and calibration parameters for the numerical modelling. The resulting model introduces numerical representation of collagen fibril density and its related regional variation, interlamellar cohesion and age-related stiffening in an anisotropic model of the human cornea. Further, the model incorporates previous modelling developments including representation of lamellae anisotropy and stiffness of the underlying matrix. Wide angle X-ray scattering has provided measured data which quantifies relative fibril anisotropy in the 2D domain. Accurate numerical description of material response to deformation is essential to providing representative simulations of corneal behaviour. Representing experimentally obtained 2D anisotropy and regional density variation in the 3D domain is an essential component of this accuracy. The constitutive model was incorporated into finite element analysis. Combining with inverse analysis, the model was calibrated to an extensive experimental database of ex vivo corneal inflation tests and ex vivo corneal shear tests. This model represents stiffness of the underlying matrix which is 2–3 orders of magnitude than the mechanical response representing the collagen fibrils in the lamellae. The presented model, along with its age dependent material coefficients, allows finite element modelling for an individual patient with material stiffness approximated based on their age. This has great potential to be used in both daily clinical practice for the planning and optimisation of corrective procedures and in pre-clinical optimisation of diagnostic procedures