Biomechanics of the keratoconic cornea: Theory, segmentation, pressure distribution, and coupled FE-optimization algorithm

Impacts and Correlations on Corneal Biomechanics, Corneal Optical Density and Intraocular Pressure after Cataract Surgery

The study aimed to investigate the extended effects and interrelations of corneal biomechanics, corneal optical density (COD), corneal thickness (CT), and intraocular pressure (IOP) following cataract surgery. Sixteen eyes were analyzed prospectively. The Corneal Visualization Scheimpflug Technology (Corvis ST) device assessed corneal biomechanics, while the Pentacam AxL® (Pentacam) measured COD and CT. Postoperative data were collected around six months after surgery, with a subgroup analysis of data at nine months. The Pearson correlation was used to examine the relationship between surgical-induced changes in corneal biomechanics and COD. At six months, significant postoperative differences were observed in various biomechanical indices, including uncorrected IOP (IOPuct) and biomechanics-corrected IOP (bIOP). However, many indices lost statistical significance by the nine-month mark, suggesting the reversibility of postoperative corneal changes. Postoperative COD increased at the anterior layer of the 2-6 mm annulus and incision site. The changes in COD correlated with certain biomechanical indices, including maximal (Max) deformative amplitude (DA) and stiffness parameter (SP). In conclusion, despite significant immediate postoperative changes, corneal biomechanics, COD, and IOP experienced a gradual recovery process following cataract surgery. Clinicians should maintain vigilance for any unusual changes during the short-term observation period to detect abnormalities early.

Comparative evaluation of methods for obtaining the stress-free geometry of keratoconus corneas with different levels of severity

A biomechanical model that simulates the physiological pressure load on the cornea without considering the stresses in vivo will result in an overstress or underestimation of the stress field and incorrect deformation of the corneal structure. Therefore, it is essential to propose numerical models that consider the stress-free geometry of the cornea. In this study, the Displacement and Pre-stress methods were compared to obtain the stress-free geometry (S-FG) and the physiological estimated geometry (P-EG), based on the patient-specific geometric behavior and the computational time required to reach each geometry. The same shape and contour conditions were considered in the models obtained from both methods for each of the pathological scenarios analyzed. Both methods behaved differently to obtain the free geometry, and this difference increased with the severity grade of the disease. However, they behaved in a similar way to reach the physiological estimated geometry. The Displacement method required a lower computational cost to reach the free geometry, with both methods presenting a similar computational cost to obtain the physiological geometry. The stress-free geometries obtained by both methods allowed to characterize the existing biomechanical decompensation during the progression of the diseases. In conclusion, the calculation of the stress-free corneal geometry associated to the clinically measured intraocular pressure with the Displacement and Prestress Methods in keratoconus eyes allows the development of accurate and useable models in clinical practice in real time. This displacement method shows some benefits in terms of computational cost.

Study of the Influence of Boundary Conditions on Corneal Deformation Based on the Finite Element Method of a Corneal Biomechanics Model

Implementing in silico corneal biomechanical models for surgery applications can be boosted by developing patient-specific finite element models adapted to clinical requirements and optimized to reduce computational times. This research proposes a novel corneal multizone-based finite element model with octants and circumferential zones of clinical interest for material definition. The proposed model was applied to four patient-specific physiological geometries of keratoconus-affected corneas. Free-stress geometries were calculated by two iterative methods, the displacements and prestress methods, and the influence of two boundary conditions: embedded and pivoting. The results showed that the displacements, stress and strain fields differed for the stress-free geometry but were similar and strongly depended on the boundary conditions for the estimated physiological geometry when considering both iterative methods. The comparison between the embedded and pivoting boundary conditions showed bigger differences in the posterior limbus zone, which remained closer in the central zone. The computational calculation times for the stress-free geometries were evaluated. The results revealed that the computational time was prolonged with disease severity, and the displacements method was faster in all the analyzed cases. Computational times can be reduced with multicore parallel calculation, which offers the possibility of applying patient-specific finite element models in clinical applications.

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.

The extraction and application of antisymmetric characteristics of the cornea during air-puff perturbations

BACKGROUND

A non-contact tonometer is used to measure intraocular pressure, and studies have primarily relied on apex displacements to assess corneal properties. However, previous studies have overlooked the asymmetric characteristics of lateral corneal perturbations, leading to a gap in understanding of the lateral mechanical properties and its application.

METHOD

To investigate these lateral perturbations, we designed an experiment to sequentially record the corneal profiles when two consecutive air-puffs were applied at the center of the same cornea within a short period. Moreover, we used modal decomposition to decompose anterior surface profiles into symmetric and antisymmetric modes to comprehensively analyze the asymmetric characteristics. To extract mechanical properties, we utilized high-pass frequency analysis (>250 Hz) to filter out noise and errors.

RESULTS

Symmetric modes between the two consecutive air-puffs exhibited major similarities during vibration; however, antisymmetric modes exhibited minor differences in lateral perturbations of asymmetric vibration. The antisymmetric modes might be related to air-puff misalignment and mechanical properties. Through applying frequency analysis, the mechanical properties could be proven at high frequencies and misalignment shown at low frequencies. Furthermore, we compared the corneal vibration profiles of 259 healthy participants and 50 patients with keratoconus. Their properties showed that the antisymmetric modes of the keratoconus group exhibited a completely opposite direction of deformation compared to that in the healthy group.

CONCLUSIONS

Our proposed algorithm not only extracts antisymmetric characteristics but also offers valuable insights into decompose misalignment and mechanical properties of healthy and keratoconus corneas, presenting a new perspective for corneal biomechanics.

Applied use of biomechanical measurements from human tissues for the development of medical skills trainers: a scoping review

OBJECTIVE

The objective of this review was to identify quantitative biomechanical measurements of human tissues, the methods for obtaining these measurements, and the primary motivations for conducting biomechanical research.

INTRODUCTION

Medical skills trainers are a safe and useful tool for clinicians to use when learning or practicing medical procedures. The haptic fidelity of these devices is often poor, which may be because the synthetic materials chosen for these devices do not have the same mechanical properties as human tissues. This review investigates a heterogeneous body of literature to identify which biomechanical properties are available for human tissues, the methods for obtaining these values, and the primary motivations behind conducting biomechanical tests.

INCLUSION CRITERIA

Studies containing quantitative measurements of the biomechanical properties of human tissues were included. Studies that primarily focused on dynamic and fluid mechanical properties were excluded. Additionally, studies only containing animal, in silico , or synthetic materials were excluded from this review.

METHODS

This scoping review followed the JBI methodology for scoping reviews and the Preferred Reporting Items for Systematic Reviews and Meta-Analyses extension for Scoping Reviews (PRISMA-ScR). Sources of evidence were extracted from CINAHL (EBSCO), IEEE Xplore, MEDLINE (PubMed), Scopus, and engineering conference proceedings. The search was limited to the English language. Two independent reviewers screened titles and abstracts as well as full-text reviews. Any conflicts that arose during screening and full-text review were mediated by a third reviewer. Data extraction was conducted by 2 independent reviewers and discrepancies were mediated through discussion. The results are presented in tabular, figure, and narrative formats.

RESULTS

Data were extracted from a total of 186 full-text publications. All of the studies, except for 1, were experimental. Included studies came from 33 countries, with the majority coming from the United States. Ex vivo methods were the predominant approach for extracting human tissue samples, and the most commonly studied tissue type was musculoskeletal. In this study, nearly 200 unique biomechanical values were reported, and the most commonly reported value was Young's (elastic) modulus. The most common type of mechanical test performed was tensile testing, and the most common reason for testing human tissues was to characterize biomechanical properties. Although the number of published studies on biomechanical properties of human tissues has increased over the past 20 years, there are many gaps in the literature. Of the 186 included studies, only 7 used human tissues for the design or validation of medical skills training devices. Furthermore, in studies where biomechanical values for human tissues have been obtained, a lack of standardization in engineering assumptions, methodologies, and tissue preparation may implicate the usefulness of these values.

CONCLUSIONS

This review is the first of its kind to give a broad overview of the biomechanics of human tissues in the published literature. With respect to high-fidelity haptics, there is a large gap in the published literature. Even in instances where biomechanical values are available, comparing or using these values is difficult. This is likely due to the lack of standardization in engineering assumptions, testing methodology, and reporting of the results. It is recommended that journals and experts in engineering fields conduct further research to investigate the feasibility of implementing reporting standards.

REVIEW REGISTRATION

Open Science Framework https://osf.io/fgb34.

Biomechanics of the JCT and SC Inner Wall Endothelial Cells with Their Basement Membrane Using 3D Serial Block-Face Scanning Electron Microscopy

BACKGROUND

More than ~70% of the aqueous humor exits the eye through the conventional aqueous outflow pathway that is comprised of the trabecular meshwork (TM), juxtacanalicular tissue (JCT), the inner wall endothelium of Schlemm's canal (SC). The flow resistance in the JCT and SC inner wall basement membrane is thought to play an important role in the regulation of the intraocular pressure (IOP) in the eye, but current imaging techniques do not provide enough information about the mechanics of these tissues or the aqueous humor in this area.

METHODS

A normal human eye was perfusion-fixed and a radial wedge of the TM tissue from a high-flow region was dissected. The tissues were then sliced and imaged using serial block-face scanning electron microscopy. Slices from these images were selected and segmented to create a 3D finite element model of the JCT and SC cells with an inner wall basement membrane. The aqueous humor was used to replace the intertrabecular spaces, pores, and giant vacuoles, and fluid-structure interaction was employed to couple the motion of the tissues with the aqueous humor.

RESULTS

Higher tensile stresses (0.8-kPa) and strains (25%) were observed in the basement membrane beneath giant vacuoles with open pores. The volumetric average wall shear stress was higher in SC than in JCT/SC. As the aqueous humor approached the inner wall basement membrane of SC, the velocity of the flow decreased, resulting in the formation of small eddies immediately after the flow left the inner wall.

CONCLUSIONS

Improved modeling of SC and JCT can enhance our understanding of outflow resistance and funneling. Serial block-face scanning electron microscopy with fluid-structure interaction can achieve this, and the observed micro-segmental flow patterns in ex vivo perfused human eyes suggest a hypothetical mechanism.

Developing an experimental-computational workflow to study the biomechanics of the human conventional aqueous outflow pathway

The aqueous humor actively interacts with the trabecular meshwork (TM), juxtacanalicular tissue (JCT), and Schlemm's canal (SC) through a dynamic fluid-structure interaction (FSI) coupling. Despite the fact that intraocular pressure (IOP) undergoes significant fluctuations, our understanding of the hyperviscoelastic biomechanical properties of the aqueous outflow tissues is limited. In this study, a quadrant of the anterior segment from a normal human donor eye was dynamically pressurized in the SC lumen, and imaged using a customized optical coherence tomography (OCT). The TM/JCT/SC complex finite element (FE) with embedded collagen fibrils was reconstructed based on the segmented boundary nodes in the OCT images. The hyperviscoelastic mechanical properties of the outflow tissues' extracellular matrix with embedded viscoelastic collagen fibrils were calculated using an inverse FE-optimization method. Thereafter, the 3D microstructural FE model of the TM, with adjacent JCT and SC inner wall, from the same donor eye was constructed using optical coherence microscopy and subjected to a flow load-boundary from the SC lumen. The resultant deformation/strain in the outflow tissues was calculated using the FSI method, and compared to the digital volume correlation (DVC) data. TM showed larger shear modulus (0.92 MPa) compared to the JCT (0.47 MPa) and SC inner wall (0.85 MPa). Shear modulus (viscoelastic) was larger in the SC inner wall (97.65 MPa) compared to the TM (84.38 MPa) and JCT (56.30 MPa). The conventional aqueous outflow pathway is subjected to a rate-dependent IOP load-boundary with large fluctuations. This necessitates addressing the biomechanics of the outflow tissues using hyperviscoelastic material-model. STATEMENT OF SIGNIFICANCE: While the human conventional aqueous outflow pathway is subjected to a large-deformation and time-dependent IOP load-boundary, we are not aware of any studies that have calculated the hyperviscoelastic mechanical properties of the outflow tissues with embedded viscoelastic collagen fibrils. A quadrant of the anterior segment of a normal humor donor eye was dynamically pressurized from the SC lumen with relatively large fluctuations. The TM/JCT/SC complex were OCT imaged and the mechanical properties of the tissues with embedded collagen fibrils were calculated using the inverse FE-optimization algorithm. The resultant displacement/strain in the FSI outflow model was validated versus the DVC data. The proposed experimental-computational workflow may significantly contribute to understanding of the effects of different drugs on the biomechanics of the conventional aqueous outflow pathway.

A position- and time-dependent pressure profile to model viscoelastic mechanical behavior of the brain tissue due to tumor growth

This study proposed a computational framework to calculate the resultant position- and time-dependent pressure profile on the brain tissue due to tumor growth. A finite element (FE) patch of the brain tissue was constructed and an inverse dynamic FE-optimization algorithm was used to calculate its viscoelastic mechanical properties under compressive uniaxial loading. Two patient-specific post-tumor resection FE models were input to the FE-optimization algorithm to calculate the optimized 3rd-order position-dependent and normal distribution time-dependent pressure profile parameters. The optimized viscoelastic material properties, the most suitable simulation time, and the optimized 3rd-order position- and -time-dependent pressure profiles were calculated.

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.

Biomechanical analysis of ocular diseases and its in vitro study methods

Ocular diseases are closely related to the physiological changes in the eye sphere and its contents. Using biomechanical methods to explore the relationship between the structure and function of ocular tissue is beneficial to reveal the pathological processes. Studying the pathogenesis of various ocular diseases will be helpful for the diagnosis and treatment of ocular diseases. We provide a critical review of recent biomechanical analysis of ocular diseases including glaucoma, high myopia, and diabetes. And try to summarize the research about the biomechanical changes in ocular tissues (e.g., optic nerve head, sclera, cornea, etc.) associated with those diseases. The methods of ocular biomechanics research in vitro in recent years are also reviewed, including the measurement of biomechanics by ophthalmic equipment, finite element modeling, and biomechanical analysis methods. And the preparation and application of microfluidic eye chips that emerged in recent years were summarized. It provides new inspiration and opportunity for the pathogenesis of eye diseases and personalized and precise treatment.

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.

A detailed methodology to model the Non Contact Tonometry: a Fluid Structure Interaction study

Understanding the corneal mechanical properties has great importance in the study of corneal pathologies and the prediction of refractive surgery outcomes. Non-Contact Tonometry (NCT) is a non-invasive diagnostic tool intended to characterize the corneal tissue response in vivo by applying a defined air-pulse. The biomarkers inferred from this test can only be considered as indicators of the global biomechanical behaviour rather than the intrinsic biomechanical properties of the corneal tissue. A possibility to isolate the mechanical response of the corneal tissue is the use of an inverse finite element method, which is based on accurate and reliable modelling. Since a detailed methodology is still missing in the literature, this paper aims to construct a high-fidelity finite-element model of an idealized 3D eye for in silico NCT. A fluid-structure interaction (FSI) simulation is developed to virtually apply a defined air-pulse to a 3D idealized eye model comprising cornea, limbus, sclera, lens and humors. Then, a sensitivity analysis is performed to examine the influence of the intraocular pressure (IOP) and the structural material parameters on three biomarkers associated with corneal deformation. The analysis reveals the requirements for the in silico study linked to the correct reproduction of three main aspects: the air pressure over the cornea, the biomechanical properties of the tissues, and the IOP. The adoption of an FSI simulation is crucial to capture the correct air pressure profile over the cornea as a consequence of the air-jet. Regarding the parts of the eye, an anisotropic material should be used for the cornea. An important component is the sclera: the stiffer the sclera, the lower the corneal deformation due to the air-puff. Finally, the fluid-like behavior of the humors should be considered in order to account for the correct variation of the IOP during the test which will, otherwise, remain constant. The development of a strong FSI tool amenable to model coupled structures and fluids provides the basis to find the biomechanical properties of the corneal tissue in vivo.

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.

Computational Modeling of Ophthalmic Procedures: Computational Modeling of Ophthalmic Procedures

PURPOSE

To explore how finite-element calculations can continue to contribute to diverse problems in ophthalmology and vision science, we describe our recent work on modeling the force on the peripheral retina in intravitreal injections and how that force increases with shorter, smaller gauge needles. We also present a calculation that determines the location and stress on a retinal pigment epithelial detachment during an intravitreal injection, the possibility that stress induced by the injection can lead to a tear of the retinal pigment epithelium.

BACKGROUND

Advanced computational models can provide a critical insight into the underlying physics in many surgical procedures, which may not be intuitive.

METHODS

The simulations were implemented using COMSOL Multiphysics. We compared the monkey retinal adhesive force of 18 Pa with the results of this study to quantify the maximum retinal stress that occurs during intravitreal injections.

CONCLUSIONS

Currently used 30-gauge needles produce stress on the retina during intravitreal injections that is only slightly below the limit that can create retinal tears. As retina specialists attempt to use smaller needles, the risk of complications may increase. In addition, we find that during an intravitreal injection, the stress on the retina in a pigment epithelial detachment occurs at the edge of the detachment (found clinically), and the stress is sufficient to tear the retina. These findings may guide physicians in future clinical research. NOTE: Publication of this article is sponsored by the American Ophthalmological Society.

Biomechanics of human trabecular meshwork in healthy and glaucoma eyes via dynamic Schlemm's canal pressurization

BACKGROUND AND OBJECTIVE

The trabecular meshwork (TM) consists of extracellular matrix (ECM) with embedded collagen and elastin fibers providing its mechanical support. TM stiffness is considerably higher in glaucoma eyes. Emerging data indicates that the TM moves dynamically with transient intraocular pressure (IOP) fluctuations, implying the viscoelastic mechanical behavior of the TM. However, little is known about TM viscoelastic behavior. We calculated the viscoelastic mechanical properties of the TM in n = 2 healthy and n = 2 glaucoma eyes.

METHODS

A quadrant of the anterior segment was submerged in a saline bath, and a cannula connected to an adjustable saline reservoir was inserted into Schlemm's canal (SC). A spectral domain-OCT (SD-OCT) provided continuous cross-sectional B-scans of the TM/JCT/SC complex during pressure oscillation from 0 to 30 mmHg at two locations. The TM/JCT/SC complex boundaries were delineated to construct a 20-µm-thick volume finite element (FE) mesh. Pre-tensioned collagen and elastin fibrils were embedded in the model using a mesh-free penalty-based cable-in-solid algorithm. SC pressure was represented by a position- and time-dependent pressure boundary; floating boundary conditions were applied to the other cut edges of the model. An FE-optimization algorithm was used to adjust the ECM/fiber mechanical properties such that the TM/JCT/SC model and SD-OCT imaging data best matched over time.

RESULTS

Significantly larger short- and long-time ECM shear moduli (p = 0.0032), and collagen (1.82x) and elastin (2.72x) fibril elastic moduli (p = 0.0001), were found in the TM of glaucoma eyes compared to healthy controls.

CONCLUSIONS

These findings provide additional clarity on the mechanical property differences in healthy and glaucomatous outflow pathway under dynamic loading. Understanding the viscoelastic properties of the TM may serve as a new biomarker in early diagnosis of glaucoma.

Modeling the biomechanics of the conventional aqueous outflow pathway microstructure in the human eye

BACKGROUND AND OBJECTIVE

Intraocular pressure (IOP) is determined by aqueous humor outflow resistance, which is a function of the combined resistance of Schlemm's canal (SC) endothelium and the trabecular meshwork (TM) and their interactions in the juxtacanalicular connective tissue (JCT) region. Aqueous outflow in the conventional outflow pathway results in pressure gradient across the TM, JCT, and SC inner wall, and induces mechanical stresses and strains that influence the geometry and homeostasis of the outflow system. The outflow resistance is affected by alteration in tissues' geometry, so there is potential for active, two-way, fluid-structure interaction (FSI) coupling between the aqueous humor (fluid) and the TM, JCT, and SC inner wall (structure). However, our understanding of the biomechanical interactions of the aqueous humor with the outflow connective tissues and its contribution to the outflow resistance regulation is incomplete.

METHODS

In this study, a microstructural finite element (FE) model of a human eye TM, JCT, and SC inner wall was constructed from a segmented, high-resolution histologic 3D reconstruction of the human outflow system. Three different elastic moduli (0.004, 0.128, and 51.5 MPa based on prior reports) were assigned to the TM/JCT complex while the elastic modulus of the SC inner wall was kept constant at 0.00748 MPa. The hydraulic conductivity was programmed separately for the TM, JCT, and SC inner wall using a custom subroutine. Cable elements were embedded into the TM and JCT extracellular matrix to represent the directional stiffness imparted by anisotropic collagen fibril orientation. The resultant stresses and strains in the outflow system were calculated using fluid-structure interaction method.

RESULTS

The higher TM/JCT stiffness resulted in larger stresses, but smaller strains in the outflow connective tissues, and resulted in a 4- and 5-fold larger pressure drop across the SC inner wall, respectively, compared to the most compliant model. Funneling through µm-sized SC endothelial pores was evident in the models at lower tissue stiffness, but aqueous flow was more turbulent in models with higher TM/JCT stiffness.

CONCLUSIONS

The mechanical properties of the outflow tissues play a crucial role in the hydrodynamics of the aqueous humor in the conventional outflow system.

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

BACKGROUND AND OBJECTIVE

Accurate finite element (FE) simulation of the optic nerve head (ONH) depends on accurate mechanical properties of the load-bearing tissues. The peripapillary sclera in the ONH exhibits a depth-dependent, anisotropic, heterogeneous collagen fiber distribution. This study proposes a novel cable-in-solid modeling approach that mimics heterogeneous anisotropic collagen fiber distribution, validates the approach against published experimental biaxial tensile tests of scleral patches, and demonstrates its effectiveness in a complex model of the posterior human eye and ONH.

METHODS

A computational pipeline was developed that defines control points in the sclera and pia mater, distributes the depth-dependent circumferential, radial, and isotropic cable elements in the sclera and pia in a pattern that mimics collagen fiber orientation, and couples the cable elements and solid matrix using a mesh-free penalty-based cable-in-solid algorithm. A parameter study was performed on a model of a human scleral patch subjected to biaxial deformation, and computational results were matched to published experimental data. The new approach was incorporated into a previously published eye-specific model to test the method; results were then interpreted in relation to the collagen fibers' (cable elements) role in the resultant ONH deformations, stresses, and strains.

RESULTS

Results show that the cable-in-solid approach can mimic the full range of scleral mechanical behavior measured experimentally. Disregarding the collagen fibers/cable elements in the posterior eye model resulted in ∼20-60% greater tensile and shear stresses and strains, and ∼30% larger posterior deformations in the lamina cribrosa and peripapillary sclera.

CONCLUSIONS

The cable-in-solid approach can easily be implemented into commercial FE packages to simulate the heterogeneous and anisotropic mechanical properties of collagenous biological tissues.

Exploring the Biomechanical Properties of the Human Cornea In Vivo Based on Corvis ST

Purpose: The aim of this study was to provide a method to determine corneal nonlinear viscoelastic properties based on the output data of corneal visualization Scheimpflug technology (Corvis ST).

Methods: The Corvis ST data from 18 eyes of 12 healthy humans were collected. Based on the air-puff pressure and the corneal displacement from the Corvis ST test of normal human eyes, the work done by the air-puff attaining the whole corneal displacement was obtained. By applying a visco-hyperelastic strain energy density function of the cornea, in which the first-order Prony relaxation function and the first-order Ogden strain energy were employed, the corneal strain energy during the Corvis ST test was calculated. Then the work done by the air-puff attaining the whole corneal displacement was completely regarded as the strain energy of the cornea. The identification of the nonlinear viscoelastic parameters was carried out by optimizing the sum of difference squares of the work and the strain energy using the genetic algorithm.

Results: The visco-hyperelastic model gave a good fit to the data of corneal strain energy with time during the Corvis ST test (R² > 0.95). The determined Ogden model parameter μ ranged from 0.42 to 0.74 MPa, and α ranged from 32.76 to 55.63. The parameters A and τ in the first-order Prony function were 0.09–0.36 and 1.21–1.95 ms, respectively.

Conclusion: It is feasible to determine the corneal nonlinear viscoelastic properties based on the corneal contour information and air-puff pressure of the Corvis ST test.

Inverse dynamic finite element-optimization modeling of the brain tumor mass-effect using a variable pressure boundary

BACKGROUND AND OBJECTIVE

Statistical atlases of brain structure can potentially contribute in the surgical and radiotherapeutic treatment planning for the brain tumor patients. However, the current brain image-registration methods lack of accuracy when it comes to the mass-effect caused by tumor growth. Numerical simulations, such as finite element method (FEM), allow us to calculate the resultant pressure and deformation in the brain tissue due to tumor growth, and to predict the mass-effect. To date, however, the pressure boundary in the brain tissue due to tumor growth has been simply presented as a constant profile throughout the entire tumor outer surface that resulted in discrepancy between the patient imaging data and brain atlases.

METHODS

In this study, we employed a fully-coupled inverse dynamic FE-optimization method to estimate the resultant variable pressure boundary due to tumor resection surgery. To do that, magnetic resonance imaging data of two patients' pre- and post-tumor resection surgery were registered, segmented, volume-meshed, and prepared for fully-coupled inverse dynamic FE-optimization simulations. Two different pressure boundaries were defined on the brain cavity after tumor resection including: a) a constant pressure boundary and b) a variable pressure boundary. The inverse FE-optimization algorithm was used to find the optimum constant and variable pressure boundaries that result in the least distance between the surface-nodes of the post-surgery brain cavity and pre-surgery tumor.

RESULTS

The results revealed that a variable pressure boundary causes a considerably lower mean percentage error compared to a constant pressure one; hence, it can more effectively address the realistic boundary in tumor resection surgery and predict the mass-effect.

CONCLUSIONS

The proposed variable pressure boundary can be a robust tool that allows batch processing to register the brains with tumors to statistical atlases of normal brains and construction of brain tumor atlases. This approach is also computationally inexpensive and can be coupled to any FE software to run. The findings of this study have implications for not only predicting the accurate pressure boundary and mass-effect before tumor resection surgery, but also for predicting some clinical symptoms of brain cancers and presenting useful tools for APPLICATIONs in image-guided neurosurgery.