0
Research Papers

# Biomechanical Evaluations of Ocular Injury Risk for Blast LoadingOPEN ACCESS

[+] Author and Article Information
Bahram Notghi

Department of Mechanical Engineering,
Johns Hopkins University,
Baltimore, MD 21218
e-mail: notghi@jhu.edu

Rajneesh Bhardwaj

Department of Mechanical Engineering,
Indian Institute of Technology Bombay,
Powai, Mumbai 400076, India

Shantanu Bailoor, Thao D. Nguyen

Department of Mechanical Engineering,
Johns Hopkins University,
Baltimore, MD 21218

Kimberly A. Thompson

Weapons and Materials Research Directorate,
Army Research Laboratory,
Aberdeen Proving Ground, MD 21005

Ashley A. Weaver, Joel D. Stitzel

VT-WFU Center for Injury Biomechanics,
Wake Forest University School of Medicine,
Winston-Salem, NC 27101

1Corresponding author.

Manuscript received December 25, 2016; final manuscript received May 20, 2017; published online June 28, 2017. Assoc. Editor: Barclay Morrison.This work is in part a work of the U.S. Government. ASME disclaims all interest in the U.S. Government's contributions.

J Biomech Eng 139(8), 081010 (Jun 28, 2017) (9 pages) Paper No: BIO-16-1540; doi: 10.1115/1.4037072 History: Received December 25, 2016; Revised May 20, 2017

## Abstract

Ocular trauma is one of the most common types of combat injuries resulting from the exposure of military personnel with improvised explosive devices. The injury mechanism associated with the primary blast wave is poorly understood. We employed a three-dimensional computational model, which included the main internal ocular structures of the eye, spatially varying thickness of the cornea-scleral shell, and nonlinear tissue properties, to calculate the intraocular pressure and stress state of the eye wall and internal ocular structure caused by the blast. The intraocular pressure and stress magnitudes were applied to estimate the injury risk using existing models for blunt impact and blast loading. The simulation results demonstrated that blast loading can induce significant stresses in the different components of the eyes that correlate with observed primary blast injuries in animal studies. Different injury models produced widely different injury risk predictions, which highlights the need for experimental studies evaluating mechanical and functional damage to the ocular structures caused by the blast loading.

<>

## Introduction

The increasing use of explosive weaponry in military conflicts and terrorist attacks has led to an increase in the incidence of combat-related blast injuries sustained by soldiers and civilians. The ratio of ocular traumatic injuries to all injuries during Operation Desert Storm was nearly six times larger than in World War II [14] and made ocular trauma the fourth most common injury related to military deployment [5]. Blast injuries can be separated into four categories: primary from the blast overpressure, secondary from propelled fragments, tertiary from blunt impact, and quaternary from burns and other effects [6]. While secondary, tertiary, and quaternary injury mechanisms can be identified within the military's casualty care system, mechanisms unique to primary blast injuries are still poorly understood [7]. Computational modeling studies and animal studies show that the risk of ocular blast injury correlates with the high-pressure shock front, referred to as the positive phase, and the subsequent lower subatmospheric pressure, referred to as the negative phase of the blast wave [8], as well as reflections of the blast wave from the orbit [9] and facial features surrounding the eye [10]. There is, however, a dearth of clinical data that could verify these findings and establish the mechanism of the injury. Measuring and assessing the influence of these factors are difficult because survivable primary blast injuries are likely accompanied by injuries from fragments and blunt force trauma and are thus more difficult to distinguish and enumerate. Moreover, the severity of the blast injuries and the distance of the care facility from the injury site mean that often patients are unable to recount the injury event, and witnesses are unavailable. There has been an ongoing effort to identify the mechanism of injuries in recent years. Hines-Beard et al. [11] tested the effects of blast overpressure on mice, by firing short bursts of pressurized air through paintball gun barrels, which inflicted closed-eye injuries with features similar to those seen in patients with ocular blast trauma. Alphonse et al. [12] performed an experiment to investigate the effect of low-pressure blasts from fireworks and gunpowder charges on human cadaver eyes. They reported minor corneal abrasion caused by propelled fragmented and found a low risk of severe ocular damage. Sherwood et al. [13] studied the tissue damage in enucleated porcine eyes caused by a blast wave generated by a shock tube and reported angle recession, internal scleral delamination, cyclodialysis, and peripheral chorioretinal detachments. Bhardwaj et al. [10] demonstrated, through a computational modeling study, the strong influence of facial features on blast pressure loading to the eye and the internal ocular structures on the biomechanical response. Due to the importance of the facial feature and the internal ocular structures on the biomechanical response, and challenges to performing the experimental studies in full scale, computational modeling can improve the understanding of the mechanism of the primary injuries.

Experimental studies of blunt impact have developed statistical models correlating injury risks to the kinetic energy, velocity, mass, and normalized energy of the projectile [2326]. Compared to experimental studies, computational studies can quantify the global responses (i.e., intraocular pressure) of the globe and the local stress and deformation state of the intraocular tissues with more detail [27]. Validated computational models of the eye have been used to analyze the risk of injuries from eye impacts of projectiles with different mass, velocity, and size [15,27,28]. Weaver et al. [27] used a large number of experimental eye impact tests collected from literature to construct a model and analyzed the global and localized responses of the eye to a variety of blunt projectile impacts and evaluated the risk of eye injuries. Alphonse et al. [12] used a correlation between the intraocular pressure (IOP) and normalized energy based on the assumption that the projected area of an unprotected eye is equivalent to a 11.16 mm diameter of a projectile. The authors used the IOP as a metric to calculate the risk of injuries from firework overpressure [10]. Sherwood et al. [13] developed a score-based system using terminology consistent with the Birmingham Eye Trauma Terminology (BETT) [29] based on the specific impulse which was calculated by timewise integration of the positive blast pressure waveform.

The goal of this study is to understand the mechanism of primary ocular injury based on stress analysis and also compare the probability of the injury based on three injury risk models. We extended the computational model of Bhardwaj et al. [10] to include the sclera, cornea, limbus, aqueous, lens, ciliary zonule, ciliary muscle, vitreous, retina, choroid, lamina cribrosa (LC), prelaminar neural tissue (PLNT), and extraocular tissue. We compared the injury risks predicted by Kennedy–Duma [30] and Weaver et al. [20] which used IOP and the maximum principal stress and the injury risk developed by Sherwood et al. [13] which used specific impulse of the positive blast phase.

## Materials and Methods

We developed a three-dimensional computational model composed of a deformable globe with internal ocular structures, rigid skull, and fluid domain. A finite difference compressible flow solver was used to simulate the propagation of the blast wave in the fluid domain [10], while a finite element elastodynamic solver was used to calculate the large deformation and stress response of tissue components of the globe. A sharp-interface immersed boundary method was used to describe the fluid–structure interaction at the interface of the solid and fluid domains [31].

Figure 1 shows a schematic of the computational domain. We considered a 2.0 kg TNT charge mass located Lex = 2.5 m in front of the face and applied the initial boundary condition at Lin = 1.8 m from the charge. The charge mass and location were obtained from the conditions of field blast tests performed by the Army Test Center [32], and these conditions were also used in our previous computational modeling study on the effect of facial features on blast wave reflections around the eye [10]. The rigid skull is meshed using bilinear triangular surface elements, while the deformable eye has trilinear hexahedral elements.

###### Finite Element Models of the Head and Eye.

The model of the rigid skull was described in detail by Bhardwaj et al. [10]. The deformable model for the globe was positioned inside the orbit of the rigid skull model, such that the lateral protrusion (LP = 12.0 mm) of lateral distance (LD = 19.0 mm) of the globe (Fig. 2) was representative of the averaged measurements for a 21-yr-old male [33,34]. The deformable model for the globe was created using CUBIT©, a geometry and mesh generation toolkit developed at Sandia National Laboratories, and meshed using eight-node hexahedral elements. The model included the important intraocular components, including the sclera, limbus, cornea, aqueous humor, lens, ciliary zonule, ciliary muscles, vitreous humor, retina, choroid, LC, PLNT, and the surrounding orbital/fatty tissue (Fig. 2). The surrounding orbital tissue was described as a sphere of a space-filling, homogeneous material, and the attachment location to orbital tissue to the globe was determined from the measurements for a 25-yr-old female, reported by Schutte et al. [35]. The extraocular orbital tissue was assumed to be perfectly bonded to the globe and to the orbit of the rigid skull. The latter was accomplished by fixing the nodes on the outer surface of the extraocular tissue adjacent to the rigid skull model. The finite element model for the individual intraocular components is described below in detail.

Cornea and Sclera: The corneoscleral shell is the stiff outer layer of the eye wall that protects the intraocular components from external loading. The sclera and cornea were described as spherical shells with outer radii of 12.0 mm and 7.8 mm, respectively. The center of the sphere representing the cornea was positioned 5.0 mm [15] anterior to the center of the sclera. The thickness of the corneoscleral shell was varied continuously from the cornea to the limbus. The corneal thickness was assumed to increase linearly from 0.52 mm at the apex to 0.66 mm at the limbus [15]. The thickness of the limbus also increased linearly from 0.66 to 0.8 mm [15].

The thickness of the sclera was represented using the micro-MRI measurements of Norman et al. [36] for 11 human eyes. In the experimental study [36], the eyes were sectioned into 15 equal-width slices along the anterior–posterior axis. We used the same method to divide the model of the sclera into 12 sections corresponding to the first 12 slices as in the study of Norman et al. [36]. In order to accommodate the eccentric position of the optic nerve head (ONH), the sectioning planes were rotated such that the anterior-most plane coincided with the plane of the posterior surface of the limbus and the posterior-most plane coincided with the plane of the anterior surface of LC. Intermediate planes were rotated through angles obtained by linear interpolation between the angles of the terminal planes. The nodes of the internal surface of the spherical shell model of the sclera were moved along the normal direction to obtain the average thickness among 11 specimens measured for the given section and quadrant [36]. The nodes on the exterior surface of the scleral shell were fixed, and the internal nodes were moved along the radial direction to achieve the local thickness measured in experiments while maintaining the uniform spacing of the nodes along the radial direction. It should be noted that the model was constructed for a left eye, while the results presented by Norman et al. [36] were for a right eye. Thus, the thickness variation in the model is a mirror image of the measured variation (Fig. 3).

Retina and Choroid: In our model, the retina was extended anteriorly up to the ora serrata, while the choroid was extended up to the ciliary body. We assumed a uniform thickness of 0.24 mm for the retina from the measurements of Wagner-Schuman et al. [37] for Caucasian males and 0.35 mm for the choroid based on measurements of mean subfoveal choroidal thickness for healthy, Japanese individuals by Ikuno et al. [38].

Lens and Zonule: We approximated the geometry of the lens measurements from a young, normal accommodated lens [39]. The lens included two components: the outer cortex and the inner nucleus. Both structures were assumed to be axisymmetric about the optical axis (Fig. 4(a)). A fifth-order polynomial ($y=ax5+bx4+cx3+dx2+f$) with coefficients a = −0.00153, b = 0.01191, c = −0.02032, d = −0.07692, and f = 2.04, for anterior, and a = 0.00375, b = −0.03036, c = 0.06955, d = 0.0943, and f = −2.09, for posterior, was used to define the cortex outline from the central axis to the designated points with 0.39 mm distance from circular end cap, shown in Fig. 4(a). The periphery of the lens was approximated by a circular end cap [39]. The lens nucleus was assumed to have the same shape as the cortex. The thickness of nucleus was scaled to 62% of the thickness of the lens [40]. Finally, the position of the horizontal axis of the nucleus was offset from that of the outer cortex by 0.51 mm [39].

Ciliary Body and Zonule Fibers: In this study, we approximated the zonule fibers as one structure. The anterior attachment location of the zonule fibers to the lens was set to 1.5 mm from the lens periphery [41] (Fig. 4(b)). We also approximated the ciliary body, which includes the ciliary muscle and ciliary processes, as a single structure. The geometry, location, and variation in the thickness of the ciliary body from anterior to posterior were taken from the study of Kao et al. [42]. The ciliary body was attached directly into the choroid, approximately 3 mm anterior to the ora serrata (Fig. 4(b)).

Optic Nerve Head: The dimensions of the ONH were obtained from the literature [43,44] (Fig. 4(c)). The LC and the connective tissue that supports the optic nerve axons as they leave the intraocular space were considered in this study as a disk with an anterior radius of 0.95 mm, a curvature of 0.2 mm [43], and a thickness of 0.4 mm [44]. The shape of the ONH cup was simplified from Ref. [43] to increase mesh size for computational time purposes.

###### Material Properties.

A quasi-incompressible neo-Hookean hyperelastic model was applied to describe the mechanical behavior of all the tissue components of the globe and orbital tissues. The neo-Hookean strain energy function is defined as [45]Display Formula

(1)$ψ(C)=μ2(I1̂−3)+κ2(I3−ln(I3)−1)$

where I1, I2, and I3 are the three invariants of the Cauchy–Green deformation tensor, and $I1̂=I3−13I1$ is the deviatoric part of the first invariant. The parameters μ and κ are the shear and bulk moduli.

Table 1 lists the material properties used for each tissue component. Few studies have measured the mechanical behavior of ocular tissues under high-rate conditions representative of the blast. The majority of the mechanical properties in Table 1 were obtained from quasi-static measurements. The shear modulus for the sclera was obtained from the dynamic inflation experiments of Bisplinghoff et al. [49]. We assumed the bulk modulus of the retina to be the same as that of water. The shear modulus of the retina was obtained from contact angle measurement of Grant et al. [52]. In this study, the properties of the limbus were assumed to be the same as the cornea, and the material parameters for choroid were assumed as those for the retina. The shear modulus of aqueous was considered the same as the vitreous. The density and bulk modulus of the ciliary zonule and ciliary muscle assumed the same values as those reported for the lens. The shear modulus for the LC was estimated from the Young's modulus used by Norman et al. [54], assuming a Poisson's ratio of 0.47 [43]. The density of the PLNT was considered the same as the density of the retina. The Young's modulus and Poisson's ratio of the PLNT are obtained based on the mean value of the ranges prescribed by Sigal [53].

###### Injury Risk.

We applied three different models of an ocular injury risk to analyze the biomechanical outcomes of the blast simulations. Two of the models were constructed from the data for blunt impact tests using either the IOP or the maximum principal stress of the corneal-scleral shell as the dependent variable. The third was developed by Sherwood et al. [13] based on specific impulse of the positive blast phase.

Kennedy–Duma model: Kennedy et al. [30] used a database contains 251 individual eye impact tests reported in the literature to develop injury criteria for hyphema, lens damage, and retinal damage. A survival analysis using maximum likelihood method was performed on the data to generate the following parametric risk functions:Display Formula

(2)$RB=[1−e−xβαβ]100%$

where x is the area normalized energy of the projectile (kJ/m2), β is a dimensionless parameter that depends on the injury type, and the parameter α is expressed in terms of m2/kJ. The values of α and β were obtained from the maximum likelihood parameter estimation method for hyphema, lens damage, and retinal damage [30]. The coefficients of the risk function predicting the corneal abrasions were obtained from logistic regression [30]. To determine the risk of injury, we first determined the maximum IOP from the blast simulation to determine the area-normalized kinetic energy of the projectile using the correlation provided by Duma et al. [56] for an aluminum projectile of 11.16 mm.

VT-WFU CIB Model: Weaver et al. [20] developed an injury risk function for globe rupture using 79 computationally validated models for different types of projectiles [27]. The risk of globe rupture was estimated for maximum stress and maximum IOP for different types of projectiles, and the corresponding normalized energy was calculated based on the globe rupture risk [30]. The normalized energy was then used to estimate the risk of injuries for hyphema, lens dislocation, and retinal damage by using Eq. (2). In the Kennedy–Duma model, the area of the projectile is assumed to be 11 mm to calculate normalized energy. In contrast, the VT-WFU CIB model uses the IOP measured in the simulation to calculate a globe rupture risk and then back calculate the normalized energy and risk of the other injuries.

Sherwood Model: Sherwood et al. [13] developed an injury risk function using the ordinal logistic regression which predicts the risk of trauma level for a given tissue using the Cumulative Injury Scale Score (CIS) based on the specific impulses (Eq. (4)). The specific impulse of the positive blast phase was calculated by timewise integration of the positive phase of the pressure waveform. The CIS score 1 indicated that the eye has some damage, but should heal fully on its own. The CIS score 2 indicated that the eye has damage that will require surgery to repair, leaving chronic pathology. The CIS 3 indicated that the eye has damage that might be repairable with surgery, with severe visual loss. The CIS 4 indicated that the eye is likely damaged beyond meaningful functional repair [13]Display Formula

(3)

where logit is the inverse of the sigmoidal logistic function, RS is the probability of occurrence of the particular CIS category, j is the individual injury grade, J is the number corresponding to the type of injury, α is the respective intercept, and β is the respective slope.

## Results

The IOP, computed as the average of the three principal stresses, was maximum at a posterior location in the vitreous near the ONH (Fig. 5) throughout the simulation. The time history of the maximum IOP oscillated about an increasing mean. The oscillation was caused by wave reflections, while the overall increase in IOP was caused by the increasing blast overpressure of the positive phase duration. The structures of the posterior eye wall (e.g., sclera, retina, macula, and optic nerve head) experienced the highest maximum principal stress and the highest amplitude oscillations. In contrast, the anterior structures, including the lens, cornea, limbus, ciliary zonule, and ciliary muscle, showed significantly smaller amplitude. The maximum principal stress of the anterior structures increased rapidly to a peak, and then slowly decayed with time, which agrees qualitatively with results reported by Rossi et al. [6]. The maximum principal stress of the sclera initially occurred at a posterior location close to the ONH, while the maximum principal stress of the retina occurred at a location in the fovea. At 1.4 ms, which is 0.5 ms after the initial impact of the blast wave, the location of the maximum principal stress of the sclera shifted to a temporal location where the sclera attached to the extraorbital tissues. The reflection of the blast wave off the facial features produced an asymmetric pressure loading on the eye. This pushed the globe away from the nasal site and caused large stresses to develop at the attachment to the orbital tissue (Fig. 5(a)). This effect can be observed as a rapid increase in the maximum principal stress in the sclera (Fig. 6(b)), cornea, limbus, and ciliary body (Fig. 6(c)) at 1.4 ms. The von Mises stress, which is a distortional stress measure, also experienced a sharp increase in these at this time. The value at 1.7 ms is more than double that at 1.4 ms in the sclera, limbus, choroid, and ciliary body. In contrast, the von Mises stress of the posterior structures, such as the LC, decreased at 1.4 ms. The von Mises stress level in the retina, PLNT, and LC is lower than the other ocular components, and this can be due to the presence of the extraocular tissue which protects the PLNT, LC, and retina from excessive distortion.

###### Implications to the Eye Injury.

We applied the three injury risk models (Sec. 2.2) to estimate the likelihood of the occurrence of different types of injuries. The calculated percentages of injury risks of the VT-WFU CIB model based on maximum IOP value of 0.29 MPa were approximately 1.16%, 78%, 53%, and 48%, for globe rupture, hyphema, lens damage, and retinal damage, respectively. The calculated injuries risk of the Kennedy–Duma model, which was also based on the maximum IOP but assumed a projectile area of 11 mm2, was 8%, 0.01%, 11%, 0%, and 0%, for corneal abrasions, globe rupture, hyphema, lens damage, and retinal damage, respectively. We also calculated the injury risk of the VT-WFU CIB model based on the maximum principal stress of the cornea-scleral shell and found a 0%, 1%, <1%, and <1% risk of injury for globe rupture, hyphema, lens dislocation, and retinal damage, respectively. The summary of the injury risk evaluated using both models is provided in Table 2.

The Sherwood model for the risk of injury was used to calculate the probability of achieving a CIS injury level based on the specific impulse of the positive blast phase. The results for the probability of achieving a CIS level for angle recession, damage to anterior chamber, choroid, optic nerve head, retina, and sclera are summarized in Fig. 7. The probability of occurrence was greater than 55% for all CIS level 1 injury and greater than 15% for the most severe CIS level 4 injuries. The largest decrease in the injury risk with the CIS level was calculated for the anterior chamber, while the injury risk exceeded 98% for all CIS levels for sclera injury. The highest risk of injury was to the sclera, which experienced the highest maximum principal tensile stress and von Mises distortional stress at the location of attachment to the orbital tissues. There was also a high risk of injury to the ONH because of the high stresses that developed in the posterior location in the early stages of blast loading.

###### Sensitivity of the Model to Sclera Thickness.

In this section, we compared the model with nonuniform sclera thickness (baseline model) and model with uniform sclera thickness to show the effect of the variation of the cornea-scleral shell thickness on the risk of injuries and stress distribution. The maximum thickness of the sclera (1 mm) was chosen for the model with uniform thickness. The IOP and stresses of both models showed the same trend and nearly identical magnitudes up to 1.4 ms, when the maximum stress in the sclera shifted from a posterior location near the ONH to a temporal location near the extraocular attachment. After 1.4 ms, the model with nonuniform thickness showed a 6% larger maximum IOP (at 1.63 ms), 34% larger maximum principal stress in the sclera (at 1.75 ms), and 15% larger von Mises stress in the sclera compared to the uniform thickness model (Fig. 8). The scleral thickness at the site of maximum stress was 0.7 mm, which was 80% smaller than for the uniform thickness model.

The von Mises stress in the adjacent tissues of the limbus, choroid, and ciliary muscle was also significantly higher in the nonuniform thickness model (Fig. 9), while there was less than a 5% difference in the von Mises stresses of the two models in the more posterior tissues, such as the PLNT and the LC. The IOP and maximum principal stress of all ocular components for both models showed the same trend and nearly identical magnitudes up to 1.4 ms and after which the model with nonuniform thickness showed higher stress in the sclera and higher IOP in comparison to the model with uniform thickness. These differences between maximum values of the maximum principal stress of the sclera and the IOP lead to the underestimation of risk of injuries based on the VT-WFU CIB risk model and the Kennedy–Duma risk model using the model with uniform thickness (Table 2).

## Discussion

The maximum IOP reached a peak of 0.29 MPa (2175 mm Hg) at 1.63 ms, which is two orders of magnitude larger than the 15 mm Hg physiologic IOP for a healthy eye. Cockerham et al. [57] reported closed-eye injuries such as corneal abrasions, vitreous hemorrhage, retinal detachment, and optic nerve atrophy in 43% out of 46 documented blast-induced traumatic brain injury patients. The blast simulations indicated high tensile (maximum principal) and distortion (von Mises) stresses in the cornea, limbus, sclera, choroid, and ciliary zonule, which were consistent with the observed injuries in the study of Cockerham et al. [57] and ex vivo study of Sherwood et al. [13]. The maximum von Mises stress in the sclera, choroid, and ciliary zonule shows the possibility of internal scleral delamination [13], chorioretinal detachments [6,13], and lens dislocation [15] (Fig. 6(d)). The high maximum principal stress at the extraocular attachment is consistent with the findings of rat studies showing torn extraocular muscles in animals exposed to blast loadings [11]. We evaluated the injury risk to the ocular components using three different injury models. The VT-WFU CIB injury risk model using the maximum IOP generally estimated a higher risk of injury than the Kennedy–Duma model. The injury risks evaluated from VT-WFU CIB model based on the maximum principal stresses were smaller than those evaluated based on the IOP. The deformation and stress state of the ocular tissues under blast loading are different from blunt impact because the pressure rate and spatial distribution are different. The discrepancy in the injury risks evaluated based on different biomechanical outcomes indicated that injury risk models based on blunt impact may not be appropriate for blast conditions. Both the VT-WFU CIB and Kennedy–Duma models showed that the probability of hyphema was higher than the other types of injuries. Moreover, all models based on the maximum IOP and maximum principal stress showed that risk of globe rupture was less than 2%, which was consistent with the experimental observations that exposure to the blast overpressure mainly produced closed globe injuries.

Each model for injury risk had its own inherent limitations. The Kennedy–Duma and VT-WFU CIB models were constructed using the data from blunt impact. The Sherwood model was developed from ex vivo experiments on porcine eyes and did not take into account the wave reflections from facial features, which amplify the blast overpressure loading on the eye. Porcine eyes may have higher mechanical strength than human eyes [58], which may lead to an underestimation of the injury risks. Furthermore, we showed the effect of the variation of the corneal-scleral shell thickness on the risk of injuries and stress distribution in ocular components. The model with uniform thickness underestimated the probability of occurrence of a different type of injuries (Table 2). These findings verified the initial hypothesis that the regional thickness variation of the cornea-scleral shell can amplify the stress distribution in the sclera.

The computational model contained a number of limitations. The injury risks were evaluated for a single blast condition with a 2 kg TNT charge mass detonated at a standoff distance of 2.5 m and position in front of the face. Further studies are needed to determine the effect of blast conditions on the injury risks. However, this study determined the distribution of stresses in the intraocular tissues caused by the primary blast wave and provided qualitative information on which ocular components are more susceptible to the damage based on the three different injury models.

The tissues of the cornea and sclera were assumed to be isotropic, spatially homogeneous, and described by a neo-Hookean hyperelastic material model. This neglects the effect of anisotropy and strain stiffening, which can be significant in collagenous tissues of the eye wall, including the cornea, sclera, lamina cribrosa, and choroid. We have extended the constitutive models to include anisotropy and strain stiffening. We are currently evaluating their effects on the tissue stress state and injury risks. We used existing measurements in the literature for the material properties of the ocular components. However, the measurements were performed under quasi-static conditions, and the quasi-static stress response is typically more compliant than the high-rate response elicited under blast conditions. Moreover, the properties vary significantly in the literature for many ocular tissues, such as the cornea, sclera, and choroid. Thus, the material properties are a significant source of uncertainty in the injury risk predictions. We are currently performing a design of experiments to determine the sensitivity of the injury risk predictions to the uncertainty in the material properties. Finally, we did not consider the viscoelastic behavior of ocular tissue, which can influence the stress magnitudes and IOP.

We estimated the risk of ocular injuries using three eye injury models. However, two of those, the Kennedy–Duma model and VT-WFU CIB model, were constructed from the blunt impact test data using either the IOP or the maximum principal stress of the corneal-scleral shell as the dependent variable. The injury risk functions used in Kennedy–Duma model and VT-WFU CIB are high-level correlations of injury risk based upon apparent insult, but are not obtained based on blast experiments. Blunt impact and blast loading have different time scales and mechanics which may result in a different distribution of stresses in the cornea and sclera and IOP. The Sherwood model was constructed based on the shock-tube experiments that subjected enucleated porcine eyes potted in gelatin to a shock wave. The injury model used the specific impulse of the positive phase duration as the dependent variable rather than a biomechanical outcome, such as the IOP or the maximum principal stress of the eye wall. Thus, the Sherwood model does not account for the focusing effects of facial features, which can significantly amplify the blast overpressure [10]. Moreover, postmortem degradation of the enucleated eyes may make the tissues more susceptible to damage from the shock wave. All three models produced widely different injury risk predictions, with the blast-specific Sherwood model providing the most severe predictions. These findings indicate that more realistic blast experiments are needed to evaluate the ocular injury risks. Furthermore, the experiments should be more tightly integrated with computational models to connect the blast conditions (e.g., specific impulse) to the intraocular biomechanical environment (e.g., IOP and tissue-level stress levels) and the injury outcomes.

## Conclusion

We presented a computational model to study the propagation of the blast wave through a deforming human eye. The model included detailed descriptions of a skull with the typical features of a 21-yr-old male, the internal ocular structures, spatially varying thickness, and isotropic, finite deformation behavior of the tissue. The stress analyses in our study correlate with the observed primary blast injuries in animal studies and ex vivo blast experiments on the choroid, retina, lens, sclera, ciliary zonules, and ONH. In addition, our sensitivity analyses showed that the effect of considering variational thickness in the cornea-scleral shell has a noticeable effect on the calculated risk of injuries. The findings of our simulations support the experimental observations of injury to the choroid, retina, lens, sclera, ciliary zonules, and ONH.

## Acknowledgements

This research was supported by US Army Medical Research, Vision Research Program under Grant No. W81XWH10-1-0766. We thank the support from the Department of Defense Science, Mathematics and Research for Transformation (DoD SMART) Scholarship for Service Program. Meshes of the skull were provided by WMRD, US Army Research Laboratory, Aberdeen, MD. We thank Professor R. Mittal for helpful discussions.

• Congressionally Directed Medical Research Programs (W81XWH10-1-0766)

## References

Stone, W. , 1950, “Ocular Injuries in the Armed Forces,” J. Am. Med. Assoc., 142(3), pp. 151–152. [PubMed]
Weichel, E. D. , Colyer, M. H. , Ludlow, S. E. , Bower, K. S. , and Eiseman, A. S. , 2008, “Combat Ocular Trauma Visual Outcomes During Operations Iraqi and Enduring Freedom,” Ophthalmology, 115(12), pp. 2235–2245. [PubMed]
Heier, J. S. , Enzenauer, R. W. , Wintermeyer, S. F. , Delaney, M. , and LaPiana, F. P. , 1993, “Ocular Injuries and Diseases at a Combat Support Hospital in Support of Operations Desert Shield and Desert Storm,” Arch. Ophthalmol., 111(6), pp. 795–798. [PubMed]
Wong, T. Y. , Smith, G. S. , Lincoln, A. E. , and Tielsch, J. M. , 2000, “Ocular Trauma in the United States Army: Hospitalization Records From 1985 Through 1994,” Am. J. Ophthalmol., 129(5), pp. 645–650. [PubMed]
Weichel, E. D. , Colyer, M. H. , Bautista, C. , Bower, K. S. , and French, L. M. , 2009, “Traumatic Brain Injury Associated With Combat Ocular Trauma,” J. Head Trauma Rehabil., 24(1), pp. 41–50. [PubMed]
Rossi, T. , Boccassini, B. , Esposito, L. , Clemente, C. , Iossa, M. , Placentino, L. , and Bonora, N. , 2012, “Primary Blast Injury to the Eye and Orbit: Finite Element Modeling,” Invest. Ophthalmol. Visual Sci., 53(13), pp. 8057–8066.
Thach, A. B. , Johnson, A. J. , Carroll, R. B. , Huchun, A. , Ainbinder, D. J. , Stutzman, R. D. , Blaydon, S. M. , DeMartelaere, S. L. , Mader, T. H. , Slade, C. S. , George, R. K., Ritchey, J. P., Barnes, S. D., and Fannin L. A., 2008, “Severe Eye Injuries in the War in Iraq, 2003–2005,” Ophthalmology, 115(2), pp. 377–382. [PubMed]
Duke-Elder, S. , 1954, Concussion Injuries. Text-Book of Ophthalmology, Vol. VI: Injuries, Henry Kimptonm, London.
Scott, R. , 2011, “The Injured Eye,” Philos. Trans. R. Soc., B, 366(1562), pp. 251–260.
Bhardwaj, R. , Ziegler, K. , Seo, J. H. , Ramesh, K. , and Nguyen, T. D. , 2014, “A Computational Model of Blast Loading on the Human Eye,” Biomech. Model. Mechanobiol., 13(1), pp. 123–140. [PubMed]
Hines-Beard, J. , Marchetta, J. , Gordon, S. , Chaum, E. , Geisert, E. E. , and Rex, T. S. , 2012, “A Mouse Model of Ocular Blast Injury That Induces Closed Globe Anterior and Posterior Pole Damage,” Exp. Eye Res., 99, pp. 63–70. [PubMed]
Alphonse, V. D. , Kemper, A. R. , Strom, B. T. , Beeman, S. M. , and Duma, S. M. , 2012, “Mechanisms of Eye Injuries From Fireworks,” JAMA, 308(1), pp. 33–34. [PubMed]
Sherwood, D. , Sponsel, W. E. , Lund, B. J. , Gray, W. , Watson, R. , Groth, S. L. , Thoe, K. , Glickman, R. D. , and Reilly, M. A. , 2014, “Anatomical Manifestations of Primary Blast Ocular Trauma Observed in a Postmortem Porcine Modelprimary Blast Ocular Trauma,” Invest. Ophthalmol. Visual Sci., 55(2), pp. 1124–1132.
Uchio, E. , Ohno, S. , Kudoh, J. , Aoki, K. , and Kisielewicz, L. T. , 1999, “Simulation Model of an Eyeball Based on Finite Element Analysis on a Supercomputer,” Br. J. Ophthalmol., 83(10), pp. 1106–1111. [PubMed]
Stitzel, J. D. , Duma, S. M. , Cormier, J. M. , and Herring, I. P. , 2002, “A Nonlinear Finite Element Model of the Eye With Experimental Validation for the Prediction of Globe Rupture,” Stapp Car Crash J., 46, pp. 81–102.
Al-Sukhun, J. , Kontio, R. , and Lindqvist, C. , 2006, “Orbital Stress Analysis—Part I: Simulation of Orbital Deformation Following Blunt Injury by Finite Element Analysis Method,” J. Oral Maxillofac. Surg., 64(3), pp. 434–442. [PubMed]
Al-Sukhun, J. , Lindqvist, C. , and Kontio, R. , 2006, “Modelling of Orbital Deformation Using Finite-Element Analysis,” J. R. Soc., Interface, 3(7), pp. 255–262.
Cirovic, S. , Bhola, R. M. , Hose, D. R. , Howard, I. C. , Lawford, P. , Marr, J. E. , and Parsons, M. A. , 2006, “Computer Modelling Study of the Mechanism of Optic Nerve Injury in Blunt Trauma,” Br. J. Ophthalmol., 90(6), pp. 778–783. [PubMed]
Stitzel, J. D. , and Weaver, A. A. , 2012, “Computational Simulations of Ocular Blast Loading and Prediction of Eye Injury Risk,” ASME Paper No. SBC2012-80792.
Weaver, A. A. , Stitzel, S. M. , and Stitzel, J. D. , 2017, “Injury Risk Prediction From Computational Simulations of Ocular Blast Loading,” Biomech. Model. Mechanobiol., 16(2), pp. 463–477. [PubMed]
Esposito, L. , Clemente, C. , Bonora, N. , and Rossi, T. , 2015, “Modelling Human Eye Under Blast Loading,” Comput. Methods Biomech. Biomed. Eng., 18(2), pp. 107–115.
Bailoor, S. , Bhardwaj, R. , and Nguyen, T. , 2015, “Effectiveness of Eye Armor During Blast Loading,” Biomech. Model. Mechanobiol., 14(6), pp. 1227–1237. [PubMed]
Duma, S. M. , and Crandall, J. R. , 2000, “Eye Injuries From Airbags With Seamless Module Covers,” J. Trauma Acute Care Surg., 48(4), pp. 786–789.
Kennedy, E. A. , Ng, T. P. , McNally, C. , Stitzel, J. D. , and Duma, S. M. , 2006, “Risk Functions for Human and Porcine Eye Rupture Based on Projectile Characteristics of Blunt Objects,” Stapp Car Crash J., 50, pp. 651–671. [PubMed]
Scott, W. R. M. , Lloyd, W. C. , Benedict, J. V. , and Meredith, R. , 2000, “Ocular Injuries Due to Projectile Impacts,” Annual Conference for the Association for the Advancement of Automotive Medicine, Chicago, IL, Oct. 2–4, Vol. 44, pp. 205–217.
Duma, S. M. , Ng, T. P. , Kennedy, E. A. , Stitzel, J. D. , Herring, I. P. , and Kuhn, F. , 2005, “Determination of Significant Parameters for Eye Injury Risk From Projectiles,” J. Trauma Acute Care Surg., 59(4), pp. 960–964.
Weaver, A. A. , Kennedy, E. A. , Duma, S. M. , and Stitzel, J. D. , 2011, “Evaluation of Different Projectiles in Matched Experimental Eye Impact Simulations,” ASME J. Biomech. Eng., 133(3), p. 031002.
Gayzik, F. , Moreno, D. , Geer, C. , Wuertzer, S. , Martin, R. , and Stitzel, J. , 2011, “Development of a Full Body Cad Dataset for Computational Modeling: A Multi-Modality Approach,” Ann. Biomed. Eng., 39(10), pp. 2568–2583. [PubMed]
Kuhn, F. , Morris, R. , Witherspoon, C. , and Mester, V. , 2004, “The Birmingham Eye Trauma Terminology System (BETT),” J. Fr. d'Ophtalmologie, 27(2), pp. 206–210.
Kennedy, E. , Duma, S. , and Street, S. , 2011, “Final Report: Eye Injury Risk Functions for Human and Focus Eyes: Hyphema, Lens Dislocation, and Retinal Damage,” U.S. Army Medical Research and Materiel Command, Fort Detrick, MD, Award No. W81XWH-05-2-0055.
Mittal, R. , Dong, H. , Bozkurttas, M. , Najjar, F. , Vargas, A. , and von Loebbecke, A. , 2008, “A Versatile Sharp Interface Immersed Boundary Method for Incompressible Flows With Complex Boundaries,” J. Comput. Phys., 227(10), pp. 4825–4852. [PubMed]
Bentz, V. , and Grimm, G. , 2013, “Joint Live Fire (JLF) Final Report for Assessment of Ocular Pressure as a Result of Blast for Protected and Unprotected Eyes,” U.S. Army Aberdeen Test Center, Aberdeen Proving Ground, MD, Report No. JLF-TR-13-01.
Weaver, A. A. , Loftis, K. L. , Tan, J. C. , Duma, S. M. , and Stitzel, J. D. , 2010, “CT Based Three-Dimensional Measurement of Orbit and Eye Anthropometry,” Invest. Ophthalmol. Visual Sci., 51(10), pp. 4892–4897.
Weaver, A. A. , Loftis, K. L. , Duma, S. M. , and Stitzel, J. D. , 2011, “Biomechanical Modeling of Eye Trauma for Different Orbit Anthropometries,” J. Biomech., 44(7), pp. 1296–1303. [PubMed]
Schutte, S. , van den Bedem, S. P. , van Keulen, F. , van der Helm, F. C. , and Simonsz, H. J. , 2006, “A Finite-Element Analysis Model of Orbital Biomechanics,” Vision Res., 46(11), pp. 1724–1731. [PubMed]
Norman, R. E. , Flanagan, J. G. , Rausch, S. M. , Sigal, I. A. , Tertinegg, I. , Eilaghi, A. , Portnoy, S. , Sled, J. G. , and Ethier, C. R. , 2010, “Dimensions of the Human Sclera: Thickness Measurement and Regional Changes With Axial Length,” Exp. Eye Res., 90(2), pp. 277–284. [PubMed]
Wagner-Schuman, M. , Dubis, A. M. , Nordgren, R. N. , Lei, Y. , Odell, D. , Chiao, H. , Weh, E. , Fischer, W. , Sulai, Y. , Dubra, A. , and Carroll, J. , 2011, “Race- and Sex-Related Differences in Retinal Thickness and Foveal Pit Morphology,” Invest. Ophthalmol. Visual Sci., 52(1), pp. 625–634.
Ikuno, Y. , Kawaguchi, K. , Nouchi, T. , and Yasuno, Y. , 2010, “Choroidal Thickness in Healthy Japanese Subjects,” Invest. Ophthalmol. Visual Sci., 51(4), pp. 2173–2176.
Burd, H. J. , 2009, “A Structural Constitutive Model for the Human Lens Capsule,” Biomech. Model. Mechanobiol., 8(3), pp. 217–231. [PubMed]
Schachar, R. A. , Abolmaali, A. , and Le, T. , 2006, “Insights Into the Age-Related Decline in the Amplitude of Accommodation of the Human Lens Using a Non-Linear Finite-Element Model,” Br. J. Ophthalmol., 90(10), pp. 1304–1309. [PubMed]
Lanchares, E. , Navarro, R. , and Calvo, B. , 2012, “Hyperelastic Modelling of the Crystalline Lens: Accommodation and Presbyopia,” J. Optom., 5(3), pp. 110–120.
Kao, C.-Y. , Richdale, K. , Sinnott, L. T. , Ernst, L. E. , and Bailey, M. D. , 2011, “Semi-Automatic Extraction Algorithm for Images of the Ciliary Muscle,” Optom. Vision Sci., 88(2), pp. 275–289.
Sigal, I. A. , Flanagan, J. G. , and Ethier, C. R. , 2005, “Factors Influencing Optic Nerve Head Biomechanics,” Invest. Ophthalmol. Visual Sci., 46(11), pp. 4189–4199.
Jonas, J. B. , Berenshtein, E. , and Holbach, L. , 2004, “Lamina Cribrosa Thickness and Spatial Relationships Between Intraocular Space and Cerebrospinal Fluid Space in Highly Myopic Eyes,” Invest. Ophthalmol. Visual Sci., 45(8), pp. 2660–2665.
Holzapfel, G. A. , 2000, Nonlinear Solid Mechanics, Vol. 24, Wiley, Chichester, UK.
Duck, F. A. , 2013, Physical Properties of Tissues: A Comprehensive Reference Book, Academic Press, London.
Bereiter-Hahn, J. , 1995, “Probing Biological Cells and Tissues With Acoustic Microscopy,” Advances in Acoustic Microscopy, Springer, New York, pp. 79–115.
Steinert, R. F. , 2010, Cataract Surgery, Elsevier Health Sciences, Irvine, CA.
Bisplinghoff, J. A. , McNally, C. , Manoogian, S. J. , and Duma, S. M. , 2009, “Dynamic Material Properties of the Human Sclera,” J. Biomech., 42(10), pp. 1493–1497. [PubMed]
Nickerson, C. S. , Park, J. , Kornfield, J. A. , and Karageozian, H. , 2008, “Rheological Properties of the Vitreous and the Role of Hyaluronic Acid,” J. Biomech., 41(9), pp. 1840–1846. [PubMed]
Mikielewicz, M. , Michael, R. , Montenegro, G. , Pinilla Cortés, L. , and Barraquer, R. , 2013, “Elastic Properties of Human Lens Zonules as a Function of Age in Presbyopes,” Acta Ophthalmol., 53(10), pp. 6109–6114.
Grant, C. A. , Twigg, P. C. , Savage, M. D. , Woon, W. H. , Wilson, M. , and Greig, D. , 2013, “Estimating the Mechanical Properties of Retinal Tissue Using Contact Angle Measurements of a Spreading Droplet,” Langmuir, 29(16), pp. 5080–5084. [PubMed]
Sigal, I. A. , 2009, “Interactions Between Geometry and Mechanical Properties on the Optic Nerve Head,” Invest. Ophthalmol. Visual Sci., 50(6), pp. 2785–2795.
Norman, R. E. , Flanagan, J. G. , Sigal, I. A. , Rausch, S. M. , Tertinegg, I. , and Ethier, C. R. , 2011, “Finite Element Modeling of the Human Sclera: Influence on Optic Nerve Head Biomechanics and Connections With Glaucoma,” Exp. Eye Res., 93(1), pp. 4–12. [PubMed]
Schoemaker, I. , Hoefnagel, P. , Mastenbroek, T. , Kolff, F. , Picken, S. , van der Helm, F. , and Simonsz, H. , 2004, “Elasticity, Viscosity and Deformation of Retrobulbar Fat in Eye Rotation,” Invest. Ophthalmol. Visual Sci., 45(13), pp. 5020–5020.
Duma, S. M. , Bisplinghoff, J. A. , Senge, D. M. , McNally, C. , and Alphonse, V. D. , 2012, “Evaluating the Risk of Eye Injuries: Intraocular Pressure During High Speed Projectile Impacts,” Curr. Eye Res., 37(1), pp. 43–49. [PubMed]
Cockerham, G. C. , Rice, T. A. , Hewes, E. H. , Cockerham, K. P. , Lemke, S. , Wang, G. , Lin, R. C. , Glynn-Milley, C. , and Zumhagen, L. , 2011, “Closed-Eye Ocular Injuries in the Iraq and Afghanistan Wars,” N. Engl. J. Med., 364(22), pp. 2172–2173. [PubMed]
Kennedy, E. A. , McNally, C. , and Duma, S. M. , 2007, “Experimental Techniques for Measuring the Biomechanical Response of the Eye During Impact,” Biomed. Sci. Instrum., 43, pp. 7–12. [PubMed]
View article in PDF format.

## References

Stone, W. , 1950, “Ocular Injuries in the Armed Forces,” J. Am. Med. Assoc., 142(3), pp. 151–152. [PubMed]
Weichel, E. D. , Colyer, M. H. , Ludlow, S. E. , Bower, K. S. , and Eiseman, A. S. , 2008, “Combat Ocular Trauma Visual Outcomes During Operations Iraqi and Enduring Freedom,” Ophthalmology, 115(12), pp. 2235–2245. [PubMed]
Heier, J. S. , Enzenauer, R. W. , Wintermeyer, S. F. , Delaney, M. , and LaPiana, F. P. , 1993, “Ocular Injuries and Diseases at a Combat Support Hospital in Support of Operations Desert Shield and Desert Storm,” Arch. Ophthalmol., 111(6), pp. 795–798. [PubMed]
Wong, T. Y. , Smith, G. S. , Lincoln, A. E. , and Tielsch, J. M. , 2000, “Ocular Trauma in the United States Army: Hospitalization Records From 1985 Through 1994,” Am. J. Ophthalmol., 129(5), pp. 645–650. [PubMed]
Weichel, E. D. , Colyer, M. H. , Bautista, C. , Bower, K. S. , and French, L. M. , 2009, “Traumatic Brain Injury Associated With Combat Ocular Trauma,” J. Head Trauma Rehabil., 24(1), pp. 41–50. [PubMed]
Rossi, T. , Boccassini, B. , Esposito, L. , Clemente, C. , Iossa, M. , Placentino, L. , and Bonora, N. , 2012, “Primary Blast Injury to the Eye and Orbit: Finite Element Modeling,” Invest. Ophthalmol. Visual Sci., 53(13), pp. 8057–8066.
Thach, A. B. , Johnson, A. J. , Carroll, R. B. , Huchun, A. , Ainbinder, D. J. , Stutzman, R. D. , Blaydon, S. M. , DeMartelaere, S. L. , Mader, T. H. , Slade, C. S. , George, R. K., Ritchey, J. P., Barnes, S. D., and Fannin L. A., 2008, “Severe Eye Injuries in the War in Iraq, 2003–2005,” Ophthalmology, 115(2), pp. 377–382. [PubMed]
Duke-Elder, S. , 1954, Concussion Injuries. Text-Book of Ophthalmology, Vol. VI: Injuries, Henry Kimptonm, London.
Scott, R. , 2011, “The Injured Eye,” Philos. Trans. R. Soc., B, 366(1562), pp. 251–260.
Bhardwaj, R. , Ziegler, K. , Seo, J. H. , Ramesh, K. , and Nguyen, T. D. , 2014, “A Computational Model of Blast Loading on the Human Eye,” Biomech. Model. Mechanobiol., 13(1), pp. 123–140. [PubMed]
Hines-Beard, J. , Marchetta, J. , Gordon, S. , Chaum, E. , Geisert, E. E. , and Rex, T. S. , 2012, “A Mouse Model of Ocular Blast Injury That Induces Closed Globe Anterior and Posterior Pole Damage,” Exp. Eye Res., 99, pp. 63–70. [PubMed]
Alphonse, V. D. , Kemper, A. R. , Strom, B. T. , Beeman, S. M. , and Duma, S. M. , 2012, “Mechanisms of Eye Injuries From Fireworks,” JAMA, 308(1), pp. 33–34. [PubMed]
Sherwood, D. , Sponsel, W. E. , Lund, B. J. , Gray, W. , Watson, R. , Groth, S. L. , Thoe, K. , Glickman, R. D. , and Reilly, M. A. , 2014, “Anatomical Manifestations of Primary Blast Ocular Trauma Observed in a Postmortem Porcine Modelprimary Blast Ocular Trauma,” Invest. Ophthalmol. Visual Sci., 55(2), pp. 1124–1132.
Uchio, E. , Ohno, S. , Kudoh, J. , Aoki, K. , and Kisielewicz, L. T. , 1999, “Simulation Model of an Eyeball Based on Finite Element Analysis on a Supercomputer,” Br. J. Ophthalmol., 83(10), pp. 1106–1111. [PubMed]
Stitzel, J. D. , Duma, S. M. , Cormier, J. M. , and Herring, I. P. , 2002, “A Nonlinear Finite Element Model of the Eye With Experimental Validation for the Prediction of Globe Rupture,” Stapp Car Crash J., 46, pp. 81–102.
Al-Sukhun, J. , Kontio, R. , and Lindqvist, C. , 2006, “Orbital Stress Analysis—Part I: Simulation of Orbital Deformation Following Blunt Injury by Finite Element Analysis Method,” J. Oral Maxillofac. Surg., 64(3), pp. 434–442. [PubMed]
Al-Sukhun, J. , Lindqvist, C. , and Kontio, R. , 2006, “Modelling of Orbital Deformation Using Finite-Element Analysis,” J. R. Soc., Interface, 3(7), pp. 255–262.
Cirovic, S. , Bhola, R. M. , Hose, D. R. , Howard, I. C. , Lawford, P. , Marr, J. E. , and Parsons, M. A. , 2006, “Computer Modelling Study of the Mechanism of Optic Nerve Injury in Blunt Trauma,” Br. J. Ophthalmol., 90(6), pp. 778–783. [PubMed]
Stitzel, J. D. , and Weaver, A. A. , 2012, “Computational Simulations of Ocular Blast Loading and Prediction of Eye Injury Risk,” ASME Paper No. SBC2012-80792.
Weaver, A. A. , Stitzel, S. M. , and Stitzel, J. D. , 2017, “Injury Risk Prediction From Computational Simulations of Ocular Blast Loading,” Biomech. Model. Mechanobiol., 16(2), pp. 463–477. [PubMed]
Esposito, L. , Clemente, C. , Bonora, N. , and Rossi, T. , 2015, “Modelling Human Eye Under Blast Loading,” Comput. Methods Biomech. Biomed. Eng., 18(2), pp. 107–115.
Bailoor, S. , Bhardwaj, R. , and Nguyen, T. , 2015, “Effectiveness of Eye Armor During Blast Loading,” Biomech. Model. Mechanobiol., 14(6), pp. 1227–1237. [PubMed]
Duma, S. M. , and Crandall, J. R. , 2000, “Eye Injuries From Airbags With Seamless Module Covers,” J. Trauma Acute Care Surg., 48(4), pp. 786–789.
Kennedy, E. A. , Ng, T. P. , McNally, C. , Stitzel, J. D. , and Duma, S. M. , 2006, “Risk Functions for Human and Porcine Eye Rupture Based on Projectile Characteristics of Blunt Objects,” Stapp Car Crash J., 50, pp. 651–671. [PubMed]
Scott, W. R. M. , Lloyd, W. C. , Benedict, J. V. , and Meredith, R. , 2000, “Ocular Injuries Due to Projectile Impacts,” Annual Conference for the Association for the Advancement of Automotive Medicine, Chicago, IL, Oct. 2–4, Vol. 44, pp. 205–217.
Duma, S. M. , Ng, T. P. , Kennedy, E. A. , Stitzel, J. D. , Herring, I. P. , and Kuhn, F. , 2005, “Determination of Significant Parameters for Eye Injury Risk From Projectiles,” J. Trauma Acute Care Surg., 59(4), pp. 960–964.
Weaver, A. A. , Kennedy, E. A. , Duma, S. M. , and Stitzel, J. D. , 2011, “Evaluation of Different Projectiles in Matched Experimental Eye Impact Simulations,” ASME J. Biomech. Eng., 133(3), p. 031002.
Gayzik, F. , Moreno, D. , Geer, C. , Wuertzer, S. , Martin, R. , and Stitzel, J. , 2011, “Development of a Full Body Cad Dataset for Computational Modeling: A Multi-Modality Approach,” Ann. Biomed. Eng., 39(10), pp. 2568–2583. [PubMed]
Kuhn, F. , Morris, R. , Witherspoon, C. , and Mester, V. , 2004, “The Birmingham Eye Trauma Terminology System (BETT),” J. Fr. d'Ophtalmologie, 27(2), pp. 206–210.
Kennedy, E. , Duma, S. , and Street, S. , 2011, “Final Report: Eye Injury Risk Functions for Human and Focus Eyes: Hyphema, Lens Dislocation, and Retinal Damage,” U.S. Army Medical Research and Materiel Command, Fort Detrick, MD, Award No. W81XWH-05-2-0055.
Mittal, R. , Dong, H. , Bozkurttas, M. , Najjar, F. , Vargas, A. , and von Loebbecke, A. , 2008, “A Versatile Sharp Interface Immersed Boundary Method for Incompressible Flows With Complex Boundaries,” J. Comput. Phys., 227(10), pp. 4825–4852. [PubMed]
Bentz, V. , and Grimm, G. , 2013, “Joint Live Fire (JLF) Final Report for Assessment of Ocular Pressure as a Result of Blast for Protected and Unprotected Eyes,” U.S. Army Aberdeen Test Center, Aberdeen Proving Ground, MD, Report No. JLF-TR-13-01.
Weaver, A. A. , Loftis, K. L. , Tan, J. C. , Duma, S. M. , and Stitzel, J. D. , 2010, “CT Based Three-Dimensional Measurement of Orbit and Eye Anthropometry,” Invest. Ophthalmol. Visual Sci., 51(10), pp. 4892–4897.
Weaver, A. A. , Loftis, K. L. , Duma, S. M. , and Stitzel, J. D. , 2011, “Biomechanical Modeling of Eye Trauma for Different Orbit Anthropometries,” J. Biomech., 44(7), pp. 1296–1303. [PubMed]
Schutte, S. , van den Bedem, S. P. , van Keulen, F. , van der Helm, F. C. , and Simonsz, H. J. , 2006, “A Finite-Element Analysis Model of Orbital Biomechanics,” Vision Res., 46(11), pp. 1724–1731. [PubMed]
Norman, R. E. , Flanagan, J. G. , Rausch, S. M. , Sigal, I. A. , Tertinegg, I. , Eilaghi, A. , Portnoy, S. , Sled, J. G. , and Ethier, C. R. , 2010, “Dimensions of the Human Sclera: Thickness Measurement and Regional Changes With Axial Length,” Exp. Eye Res., 90(2), pp. 277–284. [PubMed]
Wagner-Schuman, M. , Dubis, A. M. , Nordgren, R. N. , Lei, Y. , Odell, D. , Chiao, H. , Weh, E. , Fischer, W. , Sulai, Y. , Dubra, A. , and Carroll, J. , 2011, “Race- and Sex-Related Differences in Retinal Thickness and Foveal Pit Morphology,” Invest. Ophthalmol. Visual Sci., 52(1), pp. 625–634.
Ikuno, Y. , Kawaguchi, K. , Nouchi, T. , and Yasuno, Y. , 2010, “Choroidal Thickness in Healthy Japanese Subjects,” Invest. Ophthalmol. Visual Sci., 51(4), pp. 2173–2176.
Burd, H. J. , 2009, “A Structural Constitutive Model for the Human Lens Capsule,” Biomech. Model. Mechanobiol., 8(3), pp. 217–231. [PubMed]
Schachar, R. A. , Abolmaali, A. , and Le, T. , 2006, “Insights Into the Age-Related Decline in the Amplitude of Accommodation of the Human Lens Using a Non-Linear Finite-Element Model,” Br. J. Ophthalmol., 90(10), pp. 1304–1309. [PubMed]
Lanchares, E. , Navarro, R. , and Calvo, B. , 2012, “Hyperelastic Modelling of the Crystalline Lens: Accommodation and Presbyopia,” J. Optom., 5(3), pp. 110–120.
Kao, C.-Y. , Richdale, K. , Sinnott, L. T. , Ernst, L. E. , and Bailey, M. D. , 2011, “Semi-Automatic Extraction Algorithm for Images of the Ciliary Muscle,” Optom. Vision Sci., 88(2), pp. 275–289.
Sigal, I. A. , Flanagan, J. G. , and Ethier, C. R. , 2005, “Factors Influencing Optic Nerve Head Biomechanics,” Invest. Ophthalmol. Visual Sci., 46(11), pp. 4189–4199.
Jonas, J. B. , Berenshtein, E. , and Holbach, L. , 2004, “Lamina Cribrosa Thickness and Spatial Relationships Between Intraocular Space and Cerebrospinal Fluid Space in Highly Myopic Eyes,” Invest. Ophthalmol. Visual Sci., 45(8), pp. 2660–2665.
Holzapfel, G. A. , 2000, Nonlinear Solid Mechanics, Vol. 24, Wiley, Chichester, UK.
Duck, F. A. , 2013, Physical Properties of Tissues: A Comprehensive Reference Book, Academic Press, London.
Bereiter-Hahn, J. , 1995, “Probing Biological Cells and Tissues With Acoustic Microscopy,” Advances in Acoustic Microscopy, Springer, New York, pp. 79–115.
Steinert, R. F. , 2010, Cataract Surgery, Elsevier Health Sciences, Irvine, CA.
Bisplinghoff, J. A. , McNally, C. , Manoogian, S. J. , and Duma, S. M. , 2009, “Dynamic Material Properties of the Human Sclera,” J. Biomech., 42(10), pp. 1493–1497. [PubMed]
Nickerson, C. S. , Park, J. , Kornfield, J. A. , and Karageozian, H. , 2008, “Rheological Properties of the Vitreous and the Role of Hyaluronic Acid,” J. Biomech., 41(9), pp. 1840–1846. [PubMed]
Mikielewicz, M. , Michael, R. , Montenegro, G. , Pinilla Cortés, L. , and Barraquer, R. , 2013, “Elastic Properties of Human Lens Zonules as a Function of Age in Presbyopes,” Acta Ophthalmol., 53(10), pp. 6109–6114.
Grant, C. A. , Twigg, P. C. , Savage, M. D. , Woon, W. H. , Wilson, M. , and Greig, D. , 2013, “Estimating the Mechanical Properties of Retinal Tissue Using Contact Angle Measurements of a Spreading Droplet,” Langmuir, 29(16), pp. 5080–5084. [PubMed]
Sigal, I. A. , 2009, “Interactions Between Geometry and Mechanical Properties on the Optic Nerve Head,” Invest. Ophthalmol. Visual Sci., 50(6), pp. 2785–2795.
Norman, R. E. , Flanagan, J. G. , Sigal, I. A. , Rausch, S. M. , Tertinegg, I. , and Ethier, C. R. , 2011, “Finite Element Modeling of the Human Sclera: Influence on Optic Nerve Head Biomechanics and Connections With Glaucoma,” Exp. Eye Res., 93(1), pp. 4–12. [PubMed]
Schoemaker, I. , Hoefnagel, P. , Mastenbroek, T. , Kolff, F. , Picken, S. , van der Helm, F. , and Simonsz, H. , 2004, “Elasticity, Viscosity and Deformation of Retrobulbar Fat in Eye Rotation,” Invest. Ophthalmol. Visual Sci., 45(13), pp. 5020–5020.
Duma, S. M. , Bisplinghoff, J. A. , Senge, D. M. , McNally, C. , and Alphonse, V. D. , 2012, “Evaluating the Risk of Eye Injuries: Intraocular Pressure During High Speed Projectile Impacts,” Curr. Eye Res., 37(1), pp. 43–49. [PubMed]
Cockerham, G. C. , Rice, T. A. , Hewes, E. H. , Cockerham, K. P. , Lemke, S. , Wang, G. , Lin, R. C. , Glynn-Milley, C. , and Zumhagen, L. , 2011, “Closed-Eye Ocular Injuries in the Iraq and Afghanistan Wars,” N. Engl. J. Med., 364(22), pp. 2172–2173. [PubMed]
Kennedy, E. A. , McNally, C. , and Duma, S. M. , 2007, “Experimental Techniques for Measuring the Biomechanical Response of the Eye During Impact,” Biomed. Sci. Instrum., 43, pp. 7–12. [PubMed]

## Figures

Fig. 1

The computational domain including the fluid domain, rigid skull, and the deformable eye. The fluid domain is discretized by a structured grid. The rigid skull is discretized using bilinear triangular surface elements, and the deformable eye is discretized using trilinear hexahedral elements. Initial boundary conditions representing the charge located in front of the face with Lex = 2.5 m standoff distance (Lin) were applied Lin = 1.8 m from charge.

Fig. 2

Ocular components in exploded view including sclera, limbus, cornea, aqueous humor, lens, ciliary zonule, ciliary muscles, vitreous humor, retina, choroid, lamina cribrosa (LC), prelaminar neural tissue (PLNT), and the surrounding orbital/fatty tissue. Inset: Assembled view of the eye depicting the cornea, limbus, sclera, intraocular, and extraocular tissues.

Fig. 3

(a) Reconstructed sclera pseudosurface with contours of thickness in (left) temporal and (right) posterior views. Sectioning planes for the current eye model plotted with the sclera model in (b) nasal view. The units are in millimeter.

Fig. 4

The geometric features of (a) lens, (b) zonules, ciliary body, and (c) the optic nerve head. (a) A fifth-order polynomial was used to define the lens outline from the central axis. (I) The geometric dimensions of the axisymmetric section of the lens including cortex and nucleus. (II) Top view of the lens. (b) (I) A cross section of the revolved finite element mesh of the lens, zonules, and ciliary body. (II) A section of our model showing the connection of the lens, zonules, and ciliary body. The units are in millimeter. (c) The optic nerve head (ONH) consists of the prelaminar tissue, the lamina cribrosa, and the neural tissue. (I) The geometric dimensions of the tissues that make up the ONH are given in detail in Ref. [43]. (II) Cross section of finite element mesh of the ONH.

Fig. 5

(a) Maximum principal stress in the transverse plane of the orbit, (b) maximum intraocular pressure in the transverse plane of the orbit, and (c) maximum von Misses stress in the transverse plane of the orbit

Fig. 6

(a) Time-varying maximum intraocular pressure (IOP), (b) high amplitude time-varying maximum principal stress, (c) low amplitude time-varying maximum principal stress, and (d) time-varying maximum von Mises in different ocular structures

Fig. 7

The probability of occurrence of injuries based on calculated reflected specific impulse for different levels of CIS score level conditions

Fig. 8

Comparison of the baseline model and the model with uniform thickness for time-varying maximum (left) intraocular pressure, IOP, and (right) sclera principal stress, s1

Fig. 9

Comparison of time-varying maximum von Mises stress for the baseline model with the model with uniform thickness in sclera, limbus, choroid, and ciliary muscle

## Tables

Table 1 Material properties used for the various ocular components (w—assumed equal to that of water, r—assumed equal to that of retina, s—assumed equal to that of sclera, and l—assumed equal to that of lens)
Table 2 Injury risks calculated for model with uniform thickness (uniform) and model with nonuniform thickness (nonuniform) based on maximum IOP (MPa)

## Errata

Some tools below are only available to our subscribers or users with an online account.

### Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related Proceedings Articles
Related eBook Content
Topic Collections