Strain by virtual extensometers and video-imaging optical coherence tomography as a repeatable metric for IOP-Induced optic nerve head deformations.

PURPOSE: To determine if in vivo strain response of the Optic Nerve Head (ONH) to IOP elevation visualized using Optical Coherence Tomography (OCT) video imaging and quantified using novel virtual extensometers was able to be provided repeatable measurements of tissue specific deformations. METHODS: The ONHs of 5 eyes from 5 non-human primates (NHPs) were imaged by Spectralis OCT. A vertical and a horizontal B-scan of the ONH were continuously recorded for 60 s at 6 Hz (video imaging mode) during IOP elevation from 10 to 30 mmHg. Imaging was repeated over three imaging sessions. The 2D normal strain was computed by template-matching digital image correlation using virtual extensometers. ANOVA F-test (F) was used to compare inter-eye, inter-session, and inter-tissue variability for the prelaminar, Bruch's membrane opening (BMO), lamina cribrosa (LC) and choroidal regions (against variance the error term). F-test of the ratio between inter-eye to inter-session variability was used to test for strain repeatability across imaging sessions (FIS). RESULTS: Variability of strain across imaging session (F = 0.7263, p = 0.4855) and scan orientation was not significant (F = 1.053, p = 0.3066). Inter session variability of strain was significantly lower than inter-eye variability (FIS = 22.63, p = 0.0428) and inter-tissue variability (FIS = 99.33 p = 0.00998). After IOP elevation, strain was highest in the choroid (-18.11%, p < 0.001), followed by prelaminar tissue (-11.0%, p < 0.001), LC (-3.79%, p < 0.001), and relative change in BMO diameter (-0.57%, p = 0.704). CONCLUSIONS: Virtual extensometers applied to video-OCT were sensitive to the eye-specific and tissue-specific mechanical response of the ONH to IOP and were repeatable across imaging sessions.

Histologic validation of optical coherence tomography-based three-dimensional morphometric measurements of the human optic nerve head: Methodology and preliminary results.

PURPOSE: To compare the three-dimensional (3D) morphology of the deep load-bearing structures of the human optic nerve head (ONH) as revealed in vivo by spectral domain optical coherence tomography (SDOCT) with ex vivo quantitative 3D histology. METHODS: SDOCT imaging of the ONH was performed in six eyes from three brain-dead organ donors on life-support equipment awaiting organ procurement (in vivo conditions). Following organ procurement (ex vivo conditions), the eyes were enucleated and underwent a pars plana vitrectomy followed by pressurization to physiologic IOP and immersion fixation. Ex vivo ONH morphology was obtained from high-fidelity episcopic fluorescent 3D reconstruction. Morphologic parameters of the observed ONH canal geometry and peripapillary choroid, as well as the shape, visibility and depth of the lamina cribrosa were compared between ex vivo and in vivo measurements using custom software to align, scale, and manually delineate the different regions of the ONH. RESULTS: There was significant correspondence between in vivo and ex vivo measurements of the depth and shape of the lamina cribrosa, along with the size and shape of Bruch's membrane opening (BMO) and anterior scleral canal opening (ASCO). Weaker correspondence was observed for choroidal thickness; as expected, a thinner choroid was seen ex vivo due to loss of blood volume upon enucleation (-79.9%, p < 0.001). In addition, the lamina was shallower (-32.3%, p = 0.0019) and BMO was smaller ex vivo (-3.38%, p = 0.026), suggesting post mortem shrinkage of the fixed tissue. On average, while highly variable, only 31% of the anterior laminar surface was visible in vivo with SDOCT (p < 0.001). CONCLUSIONS: Morphologic parameters by SDOCT imaging of the deep ONH showed promising correspondence to histology metrics. Small but significant shrinkage artifact, along with large effects of exsanguination of the choroid, was seen in the ex vivo reconstructions of fixed tissues that may impact the quantification of ex vivo histoarchitecture, and this should be considered when developing models and biomarkers based on ex vivo imaging of fixed tissue. Lack of visibly of most of the lamina surface in SDOCT images is an important limitation to metrics and biomarkers based on in vivo images of the ONH deep tissues.

A Mesh-Free Approach to Incorporate Complex Anisotropic and Heterogeneous Material Properties into Eye-Specific Finite Element Models.

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.

Lamina cribrosa in glaucoma.

PURPOSE OF REVIEW: This article presents, summarizes, and interprets the most recent advances in the study and understanding of the lamina cribrosa in glaucoma, in the context of previous work. RECENT FINDINGS: The lamina is an active living structure that responds to strain by changing morphology at the micro-scale and macro-scale in glaucoma. Changes in lamina cribrosa morphology in glaucoma include posteriorization of the laminar insertion into the sclera, increased cupping or depth of the lamina cribrosa, and the development of focal lamina cribrosa defects. These lamina cribrosa changes are associated with disk hemorrhages and visual field damage, and are detectable with clinical imaging techniques such as optical coherence tomography. Glaucomatous changes in the lamina cribrosa are thought to be driven by cellular processes mediated by focal cyclical mechanical strain. Strain is eye specific and mediated by intraocular pressure, cerebrospinal fluid pressure, scleral and lamina cribrosa morphology, and structural stiffness; deleterious lamina cribrosa strains can occur at all levels of mean intraocular pressure. SUMMARY: Laminar morphology is ever changing in health and disease, and recent studies have identified several promising morphological changes that are indicative of glaucoma susceptibility, onset, and progression.

IOP telemetry in the nonhuman primate.

This review is focused on continuous IOP monitoring using telemetry systems in the nonhuman primate (NHP), presented in the context that IOP fluctuations at various timescales may be involved in glaucoma pathogenesis and progression. We use glaucoma as the primary framework to discuss how the dynamic nature of IOP might change with age, racial heritage, and disease in the context of glaucoma susceptibility and progression. We focus on the limited work that has been published in IOP telemetry in NHPs, as well as the emerging data and approaches. We review the ongoing efforts to measure continuous IOP, and the strengths, weaknesses and general pitfalls of the various approaches.

Optic nerve head biomechanics in aging and disease.

This nontechnical review is focused upon educating the reader on optic nerve head biomechanics in both aging and disease along two main themes: what is known about how mechanical forces and the resulting deformations are distributed in the posterior pole and ONH (biomechanics) and what is known about how the living system responds to those deformations (mechanobiology). We focus on how ONH responds to IOP elevations as a structural system, insofar as the acute mechanical response of the lamina cribrosa is confounded with the responses of the peripapillary sclera, prelaminar neural tissues, and retrolaminar optic nerve. We discuss the biomechanical basis for IOP-driven changes in connective tissues, blood flow, and cellular responses. We use glaucoma as the primary framework to present the important aspects of ONH biomechanics in aging and disease, as ONH biomechanics, aging, and the posterior pole extracellular matrix (ECM) are thought to be centrally involved in glaucoma susceptibility, onset and progression.

Intraocular pressure magnitude and variability as predictors of rates of structural change in non-human primate experimental glaucoma.

The purpose of this study is to determine the effects of intraocular pressure (IOP) mean, maximum and variability on the rate of structural change in experimental glaucoma. Data were taken retrospectively from 59 non-human primates involved in ongoing studies of experimental glaucoma. IOP was measured by tonometry every 1-3 weeks, and these readings split into non-overlapping fixed-length windows. First, different characterizations of IOP variability were tested to find the one that was least correlated with the mean IOP within the same window. Next, the rates of change of the Mean Position of the Disc (MPD) from confocal scanning laser tomography, and Retinal Nerve Fiber Layer Thickness (RNFLT) from spectral domain ocular coherence tomography, were calculated over each window. Mixed effects models were formed to predict these rates based on the characterizations of IOP. Normalized root mean squared residual (RMSR) from the trend of IOP during windows of five IOP measurements provided a characterization of variability showing lowest correlation with mean IOP (r < 0.001). In univariate analyses, rate of change of MPD and RNFLT were predicted by mean IOP (p < 0.001 for both) and maximum IOP (p < 0.001 for both). IOP variability did not significantly predict change in MPD (p = 0.129) or RNFLT (p = 0.438). In bivariate models, maximum IOP was the most significant predictor of change. We conclude that normalized RMSR allows the effects of IOP variability to be assessed independently of mean IOP. Maximum IOP provided the best predictability of structural change, either causally or because it captures the contributions of both mean and variability.

Lamina Cribrosa Thickening in Early Glaucoma Predicted by a Microstructure Motivated Growth and Remodeling Approach.

Glaucoma is among the leading causes of blindness worldwide. The ocular disease is characterized by irreversible damage of the retinal ganglion cell axons at the level of the lamina cribrosa (LC). The LC is a porous, connective tissue structure whose function is believed to provide mechanical support to the axons as they exit the eye on their path from the retina to the brain. Early experimental glaucoma studies have shown that the LC remodels into a thicker, more posterior structure which incorporates more connective tissue after intraocular pressure (IOP) elevation. The process by which this occurs is unknown. Here we present a microstructure motivated growth and remodeling (G&R) formulation to explore a potential mechanism of these structural changes. We hypothesize that the mechanical strain experienced by the collagen fibrils in the LC stimulates the G&R response at the micro-scale. The proposed G&R algorithm controls collagen fibril synthesis/degradation and adapts the residual strains between collagen fibrils and the surrounding tissue to achieve biomechanical homeostasis. The G&R algorithm was applied to a generic finite element model of the human eye subjected to normal and elevated IOP. The G&R simulation underscores the biomechanical need for a LC at normal IOP. The numerical results suggest that IOP elevation leads to LC thickening due to an increase in collagen fibril mass, which is in good agreement with experimental observations in early glaucoma monkey eyes. This is the first study to demonstrate that a biomechanically-driven G&R mechanism can lead to the LC thickening observed in early experimental glaucoma.

Perspectives on biomechanical growth and remodeling mechanisms in glaucoma().

Glaucoma is a blinding diseases in which damage to the axons results in loss of retinal ganglion cells. Experimental evidence indicates that chronic intraocular pressure elevation initiates axonal insult at the level of the lamina cribrosa. The lamina cribrosa is a porous collagen structure through which the axons pass on their path from the retina to the brain. Recent experimental studies revealed the extensive structural changes of the lamina cribrosa and its surrounding tissues during the development and progression of glaucoma. In this perspective paper we review the experimental evidence for growth and remodeling mechanisms in glaucoma including adaptation of tissue anisotropy, tissue thickening/thinning, tissue elongation/shortening and tissue migration. We discuss the existing predictive computational approaches that try to elucidate the potential biomechanical basis of theses growth and remodeling mechanisms and highlight open questions, challenges, and avenues for further development.

Ocular biomechanics during improvised explosive device blast: A computational study using eye-specific models.

BACKGROUND: Eye injuries comprise 10-13% of civilian improvised explosive device (IED) injuries. The bomb blast wave induces a normal and shear forces on the tissues, causing a large acute IOP elevation. This study calculated the biomechanical stresses and strains in the eye due to IED explosion via eye-specific fluid-structure interaction (FSI) models. METHODS: Blast occurred at 2, 3, and 4 m from the front and side of the victim and the weights of the IED were 1 and 2 kg. The ground was covered with the deformable soil to mimic the realistic IED explosion condition and reflect the blast wave. RESULTS: The IOP elevation of ∼6,000-48,000 mmHg was observed in the eyes while the highest IOP was occurred with the IED weight and distance of 2 kg and 2 m (front) and the lowest was occurred with the IED weight and distance of 1 kg and 4 m (side). Our findings suggest the importance of the victim location and orientation concerning the blast wave when it comes to ocular injury assessment. IOP elevation of ∼2900 and ∼2700 mmHg were observed in ∼1.6 ms after the blast for the IEDS weight of 2 kg and a victim distance of 2 m in front and side blasts, respectively, in consistence with the literature. Nonetheless, IOPs were considerably higher after ∼1.6 ms due to the merging of the bomb blast wave and its reflection off the ground. CONCLUSIONS: The stresses and strains were highest for the frontal blast. Both side and frontal blasts caused higher stresses and strains at the rectus muscle insertions where the sclera is thinnest and prone to rupture. Blast angle has no considerable role in the resultant IOP. Front blast with a heavier IED resulted a higher stresses and deformations in the eye connective tissues compared to the side blast.

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.

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.

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.

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.

Ocular Pulse Amplitude Correlates With Ocular Rigidity at Native IOP Despite the Variability in Intraocular Pulse Volume With Each Heartbeat.

Purpose: The purpose of this study was to assess ocular coat mechanical behavior using controlled ocular microvolumetric injections (MVI) of 15 µL of balanced salt solution (BSS) infused over 1 second into the anterior chamber (AC) via a syringe pump. Methods: Intraocular pressure (IOP) was continuously recorded at 200 Hz with a validated implantable IOP telemetry system in 7 eyes of 7 male rhesus macaques (nonhuman primates [NHPs]) during 5 MVIs in a series at native (3 trials), 15 and 20 mm Hg baseline IOPs, repeated in 2 to 5 sessions at least 2 weeks apart. Ocular rigidity coefficients (K) and ocular pulse volume (PV) were calculated for each trial. Data were averaged across sessions within eyes; PV was analyzed with a three-level nested ANOVA, and parameter relationships were analyzed with Pearson Correlation and linear regression. Results: After MVI at native baseline IOP of 10.4 ± 1.6 mm Hg, IOP increased by 9.1 ± 2.8 mm Hg (∆IOP) at a 9.6 ± 2.7 mm Hg/s slope, ocular pulse amplitude (OPA) was 0.70 ± 0.13 mm Hg on average; the average K was 0.042 ± 0.010 µL-1 and average PV was 1.16 ± 0.43 µL. PV varied significantly between trials, days, and eyes (P ≤ 0.05). OPA was significantly correlated with K at native IOP: Pearson coefficients ranged from 0.71 to 0.83 (P ≤ 0.05) and R2 ranged from 0.50 to 0.69 (P ≤ 0.05) during the first trial. Conclusions: The MVI-driven ∆IOP and slope can be used to assess ocular coat mechanical behavior and measure ocular rigidity. Translational Relevance: Importantly, OPA at native IOP is correlated with ocular rigidity despite the significant variability in PV between heartbeats.

Ocular biomechanics due to ground blast reinforcement.

BACKGROUND AND OBJECTIVE: Bomb blast injuries exerts a shearing force on the air-tissue interfaces, causing devastating ocular injury from the blast wave. Improvised explosive devices (IEDs) are usually placed at different heights from the ground to induce more severe injury through ground blast reinforcement (GBR). However, there is still a lack of knowledge of the role of GBR and IED height from the ground on ocular biomechanics, and how they can affect the intraocular pressure (IOP) in the eye. This study aimed to estimate the IOP due to frontal IED explosion at different heights from the ground using a fluid-structure interaction model with and without GBR effects. METHODS: A 2 kg IED was placed within 5 m of the victim at 5 different heights from the ground, including 0, 0.42, 0.85, 1.27, and 1.70 m. Two different blast formulations were used to simulate the IED explosion: (a) spherical airburst, with no amplification of the initial shock wave due to interaction with the ground-surface, and (b) hemispherical surface-burst, where the initial blast wave is immediately reflected and reinforced by the ground (GBR). RESULTS: Results revealed that the blast wave due to GBR reaches to the skull prior to the IED blast itself. The GBR also reached to the skull ∼ 0.6 ms earlier when the IED was on the ground compared to the height of 1.70 m. The highest and lowest IOPs of ∼ 17,000 and ∼ 15,000 mmHg were observed at the IED heights of 1.70 and 0 m from the ground considering GBR. However, when the role of the GBR is ignored, IOP of ∼ 9,000 mmHg was observed regardless of the IED height from the ground. The deformation in the apex of the cornea was higher when considering the GBR (∼ 0.75 cm) versus no GBR (∼ 0.65 cm). Considering GBR led to higher stresses and strains in the sclera. CONCLUSIONS: When the role of GBR was ignored, the results showed similar patterns and magnitudes of stresses and deformations in the skull and eye regardless of the height of the IED from the ground, which was not the case when GBR was considered. The findings of this study suggest the critical role of GBR in ocular blast simulations.

Analysis of the effects of finite element type within a 3D biomechanical model of a human optic nerve head and posterior pole.

BACKGROUND AND OBJECTIVE: Biomechanical stresses and strains can be simulated in the optic nerve head (ONH) using the finite element (FE) method, and various element types have been used. This study aims to investigate the effects of element type on the resulting ONH stresses and strains. METHODS: A single eye-specific model was constructed using 3D delineations of anatomic surfaces in a high-resolution, fluorescent, 3D reconstruction of a human posterior eye, then meshed using our simple meshing algorithm at various densities using 4- and 10-noded tetrahedral elements, as well as 8- and 20-noded hexahedral elements. A mesh-free approach was used to assign heterogeneous, anisotropic, hyperelastic material properties to the lamina cribrosa, sclera and pia. The models were subjected to elevated IOP of 45 mmHg after pre-stressing from 0 to 10 mmHg, and solved in the open-source FE package Calculix; results were then interpreted in relation to computational time and simulation accuracy, using the quadratic hexahedral model as the reference standard. RESULTS: The 10-noded tetrahedral and 20R-noded hexahedral elements exhibited similar scleral canal and laminar deformations, as well as laminar and scleral stress and strain distributions; the quadratic tetrahedral models ran significantly faster than the quadratic hexahedral models. The linear tetrahedral and hexahedral elements were stiffer compared to the quadratic element types, yielding much lower stresses and strains in the lamina cribrosa. CONCLUSIONS: Prior studies have shown that 20-noded hexahedral elements yield the most accurate results in complex models. Results show that 10-noded tetrahedral elements yield very similar results to 20-noded hexahedral elements and so they can be used interchangeably, with significantly lower computational time. Linear element types did not yield acceptable results.