Finite element modeling of the complex anisotropic mechanical behavior of the human sclera and pia mater

Comparing continuum and direct fiber models of soft tissues: An ocular biomechanics example reveals that continuum models may artificially disrupt the strains at both the tissue and fiber levels

Collagen fibers are the main load-bearing component of soft tissues but difficult to incorporate into models. Whilst simplified homogenization models suffice for some applications, a thorough mechanistic understanding requires accurate prediction of fiber behavior, including both detailed fiber-level strains and long-distance transmission. Our goal was to compare the performance of a continuum model of the optic nerve head (ONH) built using conventional techniques with a fiber model we recently introduced which explicitly incorporates the complex 3D organization and interaction of collagen fiber bundles [1]. To ensure a fair comparison, we constructed the continuum model with identical geometrical, structural, and boundary specifications as for the fiber model. We found that: 1) although both models accurately matched the intraocular pressure (IOP)-induced globally averaged displacement responses observed in experiments, they diverged significantly in their ability to replicate specific 3D tissue-level strain patterns. Notably, the fiber model faithfully replicated the experimentally observed depth-dependent variability of radial strain, the ring-like pattern of meridional strain, and the radial pattern of circumferential strain, whereas the continuum model failed to do so; 2) the continuum model disrupted the strain transmission along each fiber, a feature captured well by the fiber model. These results demonstrate limitations of the conventional continuum models that rely on homogenization and affine deformation assumptions, which render them incapable of capturing some complex tissue-level and fiber-level deformations. Our results show that the strengths of explicit fiber modeling help capture intricate ONH biomechanics. They potentially also help modeling other fibrous tissues. STATEMENT OF SIGNIFICANCE: Understanding the mechanics of fibrous tissues is crucial for advancing knowledge of various diseases. This study uses the ONH as a test case to compare conventional continuum models with fiber models that explicitly account for the complex fiber structure. We found that the fiber model captured better the biomechanical behaviors at both the tissue level and the fiber level. The insights gained from this study demonstrate the significant potential of fiber models to advance our understanding of not only glaucoma pathophysiology but also other conditions involving fibrous soft tissues. This can contribute to the development of therapeutic strategies across a wide range of applications.

Capturing sclera anisotropy using direct collagen fiber models. Linking microstructure to macroscopic mechanical properties

Because of the crucial role of collagen fibers on soft tissue mechanics, there is great interest in techniques to incorporate them in computational models. Recently we introduced a direct fiber modeling approach for sclera based on representing the long-interwoven fibers. Our method differs from the conventional continuum approach to modeling sclera that homogenizes the fibers and describes them as statistical distributions for each element. At large scale our method captured gross collagen fiber bundle architecture from histology and experimental intraocular pressure-induced deformations. At small scale, a direct fiber model of a sclera sample reproduced equi-biaxial experimental behavior from the literature. In this study our goal was a much more challenging task for the direct fiber modeling: to capture specimen-specific 3D fiber architecture and anisotropic mechanics of four sclera samples tested under equibiaxial and four non-equibiaxial loadings. Samples of sclera from three eyes were isolated and tested in five biaxial loadings following an approach previously reported. Using microstructural architecture from polarized light microscopy we then created specimen-specific direct fiber models. Model fiber orientations agreed well with the histological information (adjusted R2's>0.89). Through an inverse-fitting process we determined model characteristics, including specimen-specific fiber mechanical properties to match equibiaxial loading. Interestingly, the equibiaxial properties also reproduced all the non-equibiaxial behaviors. These results indicate that the direct fiber modeling method naturally accounted for tissue anisotropy within its fiber structure. Direct fiber modeling is therefore a promising approach to understand how macroscopic behavior arises from microstructure.

Comparing continuum and direct fiber models of soft tissues. An ocular biomechanics example reveals that continuum models may artificially disrupt the strains at both the tissue and fiber levels

Collagen fibers are the main load-bearing component of soft tissues but difficult to incorporate into models. Whilst simplified homogenization models suffice for some applications, a thorough mechanistic understanding requires accurate prediction of fiber behavior, including both detailed fiber-level strains and long-distance transmission. Our goal was to compare the performance of a continuum model of the optic nerve head (ONH) built using conventional techniques with a fiber model we recently introduced which explicitly incorporates the complex 3D organization and interaction of collagen fiber bundles [1]. To ensure a fair comparison, we constructed the continuum model with identical geometrical, structural, and boundary specifications as for the fiber model. We found that: 1) although both models accurately matched the intraocular pressure (IOP)-induced globally averaged displacement responses observed in experiments, they diverged significantly in their ability to replicate specific 3D tissue-level strain patterns. Notably, the fiber model faithfully replicated the experimentally observed depth-dependent variability of radial strain, the ring-like pattern of meridional strain, and the radial pattern of circumferential strain, whereas the continuum model failed to do so; 2) the continuum model disrupted the strain transmission along each fiber, a feature captured well by the fiber model. These results demonstrate limitations of the conventional continuum models that rely on homogenization and affine deformation assumptions, which render them incapable of capturing some complex tissue-level and fiber-level deformations. Our results show that the strengths of explicit fiber modeling help capture intricate ONH biomechanics. They potentially also help modeling other fibrous tissues.

Mechanical impact of epiretinal membranes on the retina utilizing finite element analysis

BACKGROUND AND OBJECTIVE

Epiretinal membrane (ERM) is a transparent membrane that forms on the surface of the neurosensory retina, causing tangential traction on the retinal surface, which may contribute to cell proliferation and contraction. Epiretinal membranes (ERMs) may be asymptomatic in some patients, while in others the membranes can progress, resulting in macular thickening and macular traction, thus distorting and inducing loss of central visual function and metamorphopsia. Currently, treatment options include follow-up or pars plana vitrectomy with an ERM peel, aiming to relieve the macular traction and improve vision and metamorphopsia. No specific criteria exist for predicting which patients might progress and need early surgery to improve and maintain good vision. The decision for surgery is based on the individual's symptoms and the physician's judgment. This study aimed to evaluate the mechanical impact in terms of stress and deformations of the ERM and to qualitatively compare them with the clinical progression of fovea thickening observed through optical coherence tomography (OCT) images.

METHODS

Numerical simulation on a three-dimensional geometrical retina and ERM model was applied to isolate factors that can be used to predict its progression and prognosis. OCT images of 14 patients with ERM were used to derive the fovea thickness progression before and after vitrectomy surgery with ERM peeling.

RESULTS

The results clearly show that the increase in ERM contractility level increases the developed stress at the fovea, which spreads and advances toward its base. The highest stress level (2.1 kPa) was developed at the highest and asymmetric contractility, producing non-uniform distributed deformations that distort the fovea structure.

CONCLUSIONS

These findings imply that high and asymmetric ERM contractility should be evaluated clinically as a factor that might signal the need for early vitrectomy surgery to avoid irreversible visual loss. Moreover, the OCT images revealed that in some cases, the thickness of the fovea indeed remains high, even after ∼12 months postoperatively, which also indicates that the deformation of the fovea in these cases is irreversible.

Fibrous finite element modeling of the optic nerve head region

The optic nerve head (ONH) region at the posterior pole of the eye is supported by a fibrous structure of collagen fiber bundles. Discerning how the fibrous structure determines the region biomechanics is crucial to understand normal physiology, and the roles of biomechanics on vision loss. The fiber bundles within the ONH structure exhibit complex three-dimensional (3D) organization and continuity across the various tissue components. Computational models of the ONH, however, usually represent collagen fibers in a homogenized fashion without accounting for their continuity across tissues, fibers interacting with each other and other fiber-specific effects in a fibrous structure. We present a fibrous finite element (FFE) model of the ONH that incorporates discrete collagen fiber bundles and their histology-based 3D organization to study ONH biomechanics as a fibrous structure. The FFE model was constructed using polarized light microscopy data of porcine ONH cryosections, representing individual fiber bundles in the sclera, dura and pia maters with beam elements and canal tissues as continuum structures. The FFE model mimics the histological in-plane orientation and width distributions of collagen bundles as well as their continuity across different tissues. Modeling the fiber bundles as linear materials, the FFE model predicts the nonlinear ONH response observed in an inflation experiment from the literature. The model also captures important microstructural mechanisms including fiber interactions and long-range strain transmission among bundles that have not been considered before. The FFE model presented here advances our understanding of the role of fibrous collagen structure in the ONH biomechanics. STATEMENT OF SIGNIFICANCE: The microstructure and mechanics of the optic nerve head (ONH) are central to ocular physiology. Histologically, the ONH region exhibits a complex continuous fibrous structure of collagen bundles. Understanding the role of the fibrous collagen structure on ONH biomechanics requires high-fidelity computational models previously unavailable. We present a computational model of the ONH that incorporates histology-based fibrous collagen structure derived from polarized light microscopy images. The model predictions agree with experiments in the literature, and provide insight into important microstructural mechanisms of fibrous tissue biomechanics, such as long-range strain transmission along fiber bundles. Our model can be used to study the microstructural basis of biomechanical damage and the effects of collagen remodeling in glaucoma.

Segmental biomechanics of the normal and glaucomatous human aqueous outflow pathway

The conventional aqueous outflow pathway, encompassing the trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and inner wall endothelium of Schlemm's canal (SC), governs intraocular pressure (IOP) regulation. This study targets the biomechanics of low-flow (LF) and high-flow (HF) regions within the aqueous humor outflow pathway in normal and glaucomatous human donor eyes, using a combined experimental and computational approach. LF and HF TM/JCT/SC complex tissues from normal and glaucomatous eyes underwent uniaxial tensile testing. Dynamic motion of the TM/JCT/SC complex was recorded using customized green-light optical coherence tomography during SC pressurization in cannulated anterior segment wedges. A hyperviscoelastic model quantified TM/JCT/SC complex properties. A fluid-structure interaction model simulated tissue-aqueous humor interaction. FluoSpheres were introduced into the pathway via negative pressure in the SC, with their motion tracked using two-photon excitation microscopy. Tensile test results revealed that the elastic moduli of the LF and HF regions in glaucomatous eyes are 3.5- and 1.5-fold stiffer than the normal eyes, respectively. The FE results also showed significantly larger shear moduli in the TM, JCT, and SC of the glaucomatous eyes compared to the normal subjects. The LF regions in normal eyes demonstrated larger elastic moduli compared to the HF regions in glaucomatous eyes. The resultant strain in the outflow tissues and velocity of the aqueous humor in the FSI models were in good agreement with the digital volume correlation and 3D particle image velocimetry data, respectively. This study uncovers stiffer biomechanical responses in glaucomatous eyes, with LF regions stiffer than HF regions in both normal and glaucomatous eyes. STATEMENT OF SIGNIFICANCE: This study delves into the biomechanics of the conventional aqueous outflow pathway, a crucial regulator of intraocular pressure and ocular health. By analyzing mechanical differences in low-flow and high-flow regions of normal and glaucomatous eyes, this research unveils the stiffer response in glaucomatous eyes. The distinction between regions' properties offers insights into aqueous humor outflow regulation, while the integration of experimental and computational methods enhances credibility. These findings have potential implications for disease management and present a vital step toward innovative ophthalmic interventions. This study advances our understanding of glaucoma's biomechanical basis and its broader impact on ocular health.

Typical localised element-specific finite element anterior eye model

Purpose

The study presents an averaged anterior eye geometry model combined with a localised material model that is straightforward, appropriate and amenable for implementation in finite element (FE) modelling.

Methods

Both right and left eye profile data of 118 subjects (63 females and 55 males) aged 22-67 years (38.5 ± 7.6) were used to build an averaged geometry model. Parametric representation of the averaged geometry model was achieved through two polynomials dividing the eye into three smoothly connected volumes. This study utilised the collagen microstructure x-ray data of 6 ex-vivo healthy human eyes, 3 right eyes and 3 left eyes in pairs from 3 donors, 1 male and 2 females aged between 60 and 80 years, to build a localised element-specific material model for the eye.

Results

Fitting the cornea and the posterior sclera sections to a 5th-order Zernike polynomial resulted in 21 coefficients. The averaged anterior eye geometry model recorded a limbus tangent angle of 37° at a radius of 6.6 mm from the corneal apex. In terms of material models, the difference between the stresses generated in the inflation simulation up to 15 mmHg in the ring-segmented material model and localised element-specific material model were significantly different (p < 0.001) with the ring-segmented material model recording average Von-Mises stress 0.0168 ± 0.0046 MPa and the localised element-specific material model recording average Von-Mises stress 0.0144 ± 0.0025 MPa.

Conclusions

The study illustrates an averaged geometry model of the anterior human eye that is easy to generate through two parametric equations. This model is combined with a localised material model that can be used either parametrically through a Zernike fitted polynomial or non-parametrically as a function of the azimuth angle and the elevation angle of the eye globe. Both averaged geometry and localised material models were built in a way that makes them easy to implement in FE analysis without additional computation cost compared to the limbal discontinuity so-called idealised eye geometry model or ring-segmented material model.

Modeling the Endothelial Glycocalyx Layer in the Human Conventional Aqueous Outflow Pathway

A layer of proteoglycans and glycoproteins known as glycocalyx covers the surface of the trabecular meshwork (TM), juxtacanalicular tissue (JCT), and Schlemm's canal (SC) inner wall of the conventional aqueous outflow pathway in the eye. This has been shown to play a role in the mechanotransduction of fluid shear stress and in the regulation of the outflow resistance. The outflow resistance in the conventional outflow pathway is the main determinant of the intraocular pressure (IOP) through an active, two-way, fluid-structure interaction coupling between the outflow tissues and aqueous humor. A 3D microstructural finite element (FE) model of a healthy human eye TM/JCT/SC complex with interspersed aqueous humor was constructed. A very thin charged double layer that represents the endothelial glycocalyx layer covered the surface of the elastic outflow tissues. The aqueous humor was modeled as electroosmotic flow that is charged when it is in contact with the outflow tissues. The electrical-fluid-structure interaction (EFSI) method was used to couple the charged double layer (glycocalyx), fluid (aqueous humor), and solid (outflow tissues). When the IOP was elevated to 15 mmHg, the maximum aqueous humor velocity in the EFSI model was decreased by 2.35 mm/s (9%) compared to the fluid-structure interaction (FSI) model. The charge or electricity in the living human conventional outflow pathway generated by the charged endothelial glycocalyx layer plays a minor biomechanical role in the resultant stresses and strains as well as the hydrodynamics of the aqueous humor.

The Effect of Intraocular Pressure Load Boundary on the Biomechanics of the Human Conventional Aqueous Outflow Pathway

BACKGROUND

Aqueous humor outflow resistance in the trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and Schlemm's canal (SC) endothelium of the conventional outflow pathway actively contribute to intraocular pressure (IOP) regulation. Outflow resistance is actively affected by the dynamic outflow pressure gradient across the TM, JCT, and SC inner wall tissues. The resistance effect implies the presence of a fluid-structure interaction (FSI) coupling between the outflow tissues and the aqueous humor. However, the biomechanical interactions between viscoelastic outflow tissues and aqueous humor dynamics are largely unknown.

METHODS

A 3D microstructural finite element (FE) model of a healthy human eye TM/JCT/SC complex was constructed with elastic and viscoelastic material properties for the bulk extracellular matrix and embedded elastic cable elements. The FE models were subjected to both idealized and a physiologic IOP load boundary using the FSI method.

RESULTS

The elastic material model for both the idealized and physiologic IOP load boundary at equal IOPs showed similar stresses and strains in the outflow tissues as well as pressure in the aqueous humor. However, outflow tissues with viscoelastic material properties were sensitive to the IOP load rate, resulting in different mechanical and hydrodynamic responses in the tissues and aqueous humor.

CONCLUSIONS

Transient IOP fluctuations may cause a relatively large IOP difference of ~20 mmHg in a very short time frame of ~0.1 s, resulting in a rate stiffening in the outflow tissues. Rate stiffening reduces strains and causes a rate-dependent pressure gradient across the outflow tissues. Thus, the results suggest it is necessary to use a viscoelastic material model in outflow tissues that includes the important role of IOP load rate.

Viscoelastic Biomechanical Properties of the Conventional Aqueous Outflow Pathway Tissues in Healthy and Glaucoma Human Eyes

BACKGROUND

Although the tissues comprising the ocular conventional outflow pathway have shown strong viscoelastic mechanical response to aqueous humor pressure dynamics, the viscoelastic mechanical properties of the trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and Schlemm's canal (SC) inner wall are largely unknown.

METHODS

A quadrant of the anterior segment from two human donor eyes at low- and high-flow (LF and HF) outflow regions was pressurized and imaged using optical coherence tomography (OCT). A finite element (FE) model of the TM, the adjacent JCT, and the SC inner wall was constructed and viscoelastic beam elements were distributed in the extracellular matrix (ECM) of the TM and JCT to represent anisotropic collagen. An inverse FE-optimization algorithm was used to calculate the viscoelastic properties of the ECM/beam elements such that the TM/JCT/SC model and OCT imaging data best matched over time.

RESULTS

The ECM of the glaucoma tissues showed significantly larger time-dependent shear moduli compared to the heathy tissues. Significantly larger shear moduli were also observed in the LF regions of both the healthy and glaucoma eyes compared to the HF regions.

CONCLUSIONS

The outflow tissues in both glaucoma eyes and HF regions are stiffer and less able to respond to dynamic IOP.

Relative Contributions of Intraocular and Cerebrospinal Fluid Pressures to the Biomechanics of the Lamina Cribrosa and Laminar Neural Tissues

Purpose

The laminar region of the optic nerve head (ONH), thought to be the site of damage to the retinal ganglion cell axons in glaucoma, is continuously loaded on its anterior and posterior surfaces by dynamic intraocular pressure (IOP) and orbital cerebrospinal fluid pressure (CSFP), respectively. Thus, translaminar pressure (TLP; TLP = IOP-CSFP) has been proposed as a glaucoma risk factor.

Methods

Three eye-specific finite element models of the posterior human eye were constructed, including full 3D microstructures of the load-bearing lamina cribrosa (LC) with interspersed laminar neural tissues (NTs), and heterogeneous, anisotropic, hyperelastic material formulations for the surrounding peripapillary sclera and adjacent pia. ONH biomechanical responses were simulated using three combinations of IOP and CSFP loadings consistent with posture change from sitting to supine.

Results

Results show that tensile, compressive, and shear stresses and strains in the ONH were higher in the supine position compared to the sitting position (P < 0.05). In addition, LC beams bear three to five times more TLP-driven stress than interspersed laminar NT, whereas laminar NT exhibit three to five times greater strain than supporting LC (P < 0.05). Compared with CSFP, IOP drove approximately four times greater stress and strain in the LC, NT, and peripapillary sclera, normalized per mm Hg pressure change. In addition, IOP drove approximately three-fold greater scleral canal expansion and anterior-posterior laminar deformation than CSFP per mm Hg (P < 0.05).

Conclusions

Whereas TLP has been hypothesized to play a prominent role in ONH biomechanics, the IOP and CSFP effects are not equivalent, as IOP-driven stress, strain, and deformation play a more dominant role than CSFP effects.