A mesh-free approach to incorporate complex anisotropic and heterogeneous material properties into eye-specific finite element models

• Published in: Computer methods in applied mechanics and engineering Vol. 358; p. 112654
• Main Authors: Grytz, Rafael, Krishnan, Kapil, Whitley, Ryan, Libertiaux, Vincent, Sigal, Ian A., Girkin, Christopher A., Downs, J. Crawford
• Format: Journal Article
• Language: English
• Published: Netherlands Elsevier B.V 01.01.2020; Elsevier BV
• Subjects: Anisotropy; Biological properties; Boundary value problems; Collagen; Eye (anatomy); Finite element method; Material properties; Mathematical models; Mesh-free; Meshless methods; Ocular biomechanics; Parameter estimation; Patient-specific modeling; Tissues; Ocular biomechanics; Mesh-free; Patient-specific modeling; patient-specific modeling; mesh-free; ocular biomechanics
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/j.cma.2019.112654

Commercial finite element modeling packages do not have the tools necessary to effectively incorporate the complex anisotropic and heterogeneous material properties typical of the biological tissues of the eye. We propose a mesh-free approach to incorporate realistic material properties into finite element models of individual human eyes. The method is based on the idea that material parameters can be estimated or measured at so called control points, which are arbitrary and independent of the finite element mesh. The mesh-free approach approximates the heterogeneous material parameters at the Gauss points of each finite element while the boundary value problem is solved using the standard finite element method. The proposed method was applied to an eye-specific model a human posterior pole and optic nerve head. We demonstrate that the method can be used to effectively incorporate experimental measurements of the lamina cribrosa micro-structure into the eye-specific model. It was convenient to define characteristic material orientations at the anterior and posterior scleral surface based on the eye-specific geometry of each sclera. The mesh-free approach was effective in approximating these characteristic material directions with smooth transitions across the sclera. For the first time, the method enabled the incorporation of the complex collagen architecture of the peripapillary sclera into an eye-specific model including the recently discovered meridional fibers at the anterior surface and the depth dependent width of circumferential fibers around the scleral canal. The model results suggest that disregarding the meridional fiber region may lead to an underestimation of local strain concentrations in the retina. The proposed approach should simplify future studies that aim to investigate collagen remodeling in the sclera and optic nerve head or in other biological tissues with similar challenges. •Combining advantages of mesh-free and mesh-based computational methods.•Effective approximation of complex anisotropic and heterogeneous material properties.•Eye-specific finite element modeling of the human eye.