Skip to main content
  • Research article
  • Open access
  • Published:

Study on establishment and mechanics application of finite element model of bovine eye



Glaucoma mainly induced by increased intraocular pressure (IOP), it was believed that the pressure that wall of eyeball withstands were determined by material properties of the tissue and stereoscopic geometry of the eyeball. In order to study the pressure changes in different parts of interior eyeball wall, it is necessary to develop a novel eye ball FEM with more accurate geometry and material properties. Use this model to study the stress changes in different parts of eyeball, especially the lamina cribrosa (LC) under normal physiological and pathological IOP, and provide a mathematical model for biomechanical studies of selected retinal ganglion cells (RGCs) death.


(1) Sclera was cut into 3.8-mm wide, 14.5-mm long strips, and cornea was cut into 9.5-mm-wide and 10-mm-long strips; (2) 858 Mini BionixII biomechanical loading instrument was used to stretch sclera and cornea. The stretching rate for sclera was 0.3 mm/s, 3 mm/s, 30 mm/s, 300 mm/s; and for cornea were 0.3 mm/s and 30 mm/s. The deformation-stress curve was recorded; (3) Naso-temporal and longitudinal distance of LC were measured; (4) Micro-CT was used to accurately scan fresh bovine eyes and obtain the geometrical image and data to establish bovine eye model. 3-D reconstruction was performed using these images and data to work out the geometric shape of bovine eye; (5) IOP levels for eyeball FEM was set and the inner wall of eyeball was used taken as load-bearing part. Simulated eyeball FE modeling was run under the IOP level of 10 mmHg, 30 mmHg, 60 mmHg and 100 mmHg, and the force condition of different parts of eyeball was recorded under different IOP levels.


(1) We obtained the material parameters more in line with physiological conditions and established a more realistic eyeball model using reversed engineering of parameters optimization method to calculate the complex nonlinear super-elastic and viscoelastic parameters more accurately; (2) We observed the following phenomenon by simulating increased pressure using FEM: as simulative IOP increased, the stress concentration scope on the posterior half of sclera became narrower; in the meantime, the stress-concentration scope on the anterior half of scleral gradually expanded, and the stress on the central part of LC is highest.


As simulative IOP increased, stress-concentration scope on the posterior half of sclera gradually narrowed; in the meantime, the stress-concentration scope on the anterior half of sclera gradually expanded, and the stress on the LC is mainly concentrated in the central part, suggesting that IOP is mainly concentrated in the anterior part of the eyeball as it increases. This might provide a biomechanical evidence to explain why RGCs in peripheral part die earlier than RGCs in central part under HIOP.

Peer Review reports


Glaucoma is a common blinding eye disease, which is induced by increased intraocular pressure (IOP) leading to problems like optic papilla sag, optic atrophy, visual field defect and loss of retinal ganglion cells (RGCs) [13]. The main risk factors of glaucoma include anatomical abnormalities stemming from inter-individual differences, intraocular pressure, local blood supply, genetic factors and refractive errors, etc. [48]. One study showed that the death of RGCs in glaucoma patients never occurred at the same time. Instead, it is a slow process, in which the subjective visual defects caused by the loss of visual sensitivity in the peripheral part of retina happens earlier than the central part [9]. Animal studies also suggested that the loss of RGCs in the central part occurs later than that of the peripheral part [10, 11]. These phenomena indicated that the damage of different RGCs varies in glaucoma, and the underlying neurobiological and biochemical mechanisms for these differences have been investigated. For example, some reports showed various distribution of excitatory neurotransmitter and neurotrophic factors in different kinds of RGCs [1220], while others showed that the blood supply was different in different parts of retina [21]. However, it remains unclear whether RGCs in different parts of retina bear the same pressure after increased IOP in glaucoma.

It is believed that the pressure on the wall of eyeball is determined by material properties of the tissue and stereoscopic geometry of the eyeball. Sclera, as one of the most important structures of the wall of eyeball, has been the focus of relative studies, which indicated that it was more likely to develop myopia when scleral creep rate was higher [2224]. Philips et al. studied the sclera of chicken with induced monocular myopia, and found that the axial length change in myopia was mainly caused by the change of scleral creep properties [25], suggesting that the geometry and material properties of eyeball may have significant impact on the physiological and pathological changes of the eye. These studies also pointed out that tissues bearing the load of optic papilla, such as sclera around optic papilla, lamina cribrosa (LC), wall of scleral canal, also bear the pressure caused by IOP. In the past two decades, the FE modeling has been widely utilized in IOP studies. For example, Bellezza et al. [26] studied the IOP-induced biomechanical changes of sclera , and Sigal et al. [27] proved that FE model had a good consistency with real eyes in geometric shape, elastic mechanical response and viscous mechanical response. Richard et al. [28] established a “Universal” FEM of human optic papilla which contained LC using the method by Sigal et al. [29]. Sensitivity analysis of the experiment revealed that the thickness of sclera around optic papilla had a strong effect on biomechanics of eyeball, which could influence the responses of LC and altered the displacement of optic papilla under all pressure levels. Due to the anatomical and functional complexity of the eyeball, it is usually difficult to experimentally study the biomechanical responses of eyeball components with detailed optic papilla, sclera and LC under both physiological and HIOP loading. It is therefore necessary to incorporate Finite Element analysis to more accurately study the subtle changes and biomechanical responses of the eyeball.

Previous studies on biomechanical model showed that the mechanical stress on optic papilla induced by increased IOP could trigger a series of problems, eventually leading to the dysfunction and death of RGCs [30]. Although the mechanisms have not been completely investigated, researches have found that LC was a vulnerable area to injury [3133], suggesting LC-related damage might be the direct cause of RGCs injury. In addition, previous studies [34, 35] showed that scleral materials, geometry and thickness significantly affected the biomechanics of optic papilla. Finite element studies of primates also showed that the sclera played a significant biomechanical role in the biomechanics of optic papilla [36, 37]. However, these studies all used “spherical sclera shell” models with uniform thickness, which over simplified the physiological and anatomical characteristics of sclera. For further studying the effects of increased IOP on function and death of RGCs, it is necessary to analyze stress changes of interior eye wall, such as different stress distribution under normal and increased IOP, or to what extent the heterogeneous thickness of eye tissues will affect stress distribution. To solve the above problems, it is necessary to establish a new FEM based on heterogeneous thickness of eyeball tissues. Therefore, in order to study the pressure changes in different parts of interior eyeball wall, including LC regions, under different IOP, in the present study we established FEMs of healthy bovine eyes, in which the realistic parameters of eyeball, especially the thickness and geometry of sclera were considered.


Sample collection

Forty bovine eyes were collected from Luxi cattle. Twenty of them, including 8 from two-year-old, 8 from three-year-old and 4 from four-year-old bovine, were used for biomechanical experiment. The other 20 eyes, including 8 from two-year-old, 8 from three-year-oldand 4 from four-year-old bovine, were used for LC measurement. The selected cattle with glossy hair, bright and transparent pupil, and spotless cornea, were healthy and met all physiological criteria. Sample collection time: 10–30 minutes after death. Respective mark above orbital cavity and nasal side were made to locate eyeball before sample collection. All experimental procedures used in the present study were approved by Ethics Committee of Xiangya School of Medicine, in accordance with the NIH guidelines for use and care of laboratory animals. All of our eyeballs were purchased from Luxi cattle slaughter house, and we have gotten the permission from the owner.

Preparation and measurement of sclera and corneal strips

Tissues and muscles around the eyeball were removed (Fig. 1a), then the eyeball was cut open from sclera-corneal junction, corneal separated from sclera and the intraocular contents were disposed. Retina and choroid were cleaned from the interior sclera. Cautions were taken not to hurt the scleral tissues. Sclera and cornea were then cut into strips. Sclera was cut into strips with an average width of 3.8 mm and length of 14.5 mm which radiated from lamina cribiosa. Cornea was cut into strips with an average width of 9.5 mm and length of 10 mm (Fig. 1b and c). It should be noted that retina and choroid are quite weak and fragile. With a thickness of only about 250 μm and Young's modulus of only 20 kPa [37], retina cannot bear load, thus our experiment only analyzed the material characteristics of load bearing tissues like sclera and cornea.

Fig. 1
figure 1

Preparation and measurement of scleral and corneal strips. a: The eyeball that was removed and cleaned up; b: Measurement of scleral width; c: Measurement of corneal width

Measurement of bovine LC

After cleaning up the tissues and muscles around eyeball and removing the retina and choroid, the eyeball was soaked in distilled water for one hour to remove prelaminar tissues and fix the sclera. Vernier caliper was used to measure the naso-temporal width of LC and the length between top and lower LC (Fig. 2).

Fig. 2
figure 2

Picture of LC measurement

Preparation of biomechanical experiment

After cutting the cornea and sclera into regular strips, they were fixed between two stainless steel plates (Fig 1b and c). Sandwiched with sample, steel plates were installed onto 858 Mini BionixII (MTS Co., Minneapolis, MN, USA), a biomechanical loading instrument (Fig. 3). The scleral strips were stretched at different speed of 0.3 mm/s, 3 mm/s, 30 mm/s and 300 mm/s, and the corneal strips were stretched at 0.3 mm/s and 30 mm/s. The deformation-force curve was recorded.

Fig. 3
figure 3

Biomechanical loading picture of eye material

Development of eyeball FE modeling

Micro-CT (ZKKS, Sharp, Osaka, Japan) has been used to accurately scan the fresh bovine eyes, with the slice thickness of 0.07 mm, and segmentation and smoothing of scleral and corneal images from different layers has been finished. Based on these, we did three-dimensional reconstruction using MIMICS (Version 12.0, Materialise Inc., Leuven, Belgium) and finally got the geometry of eyeball with heterogeneous thickness.

We meshed the eyeball geometry using HYPERMESH (Version 10.0, Altair, Troy, MI, USA). The model was divided into three key parts: sclera, cornea and LC. The units had total 14,370 nodes and 19,152 hexahedral elements, and each of elements was 0.7 to 1 mm in size. In order to ensure the accurate operation of the model, we strictly controlled the quality of elements during meshing. The minimum Jacobian Ratio was 0.47 and the minimum angle was 43.53 degrees.

In order to overcome the defects in previous experiments, namely the over simplification of eyes by FEM, such as simplifying the eyeball into a sphere, and scleral thickness into uniform thickness, we segmented 700 transverse sectional images of eyeball (Fig. 4a) and conducted three-dimensional reconstruction (Fig. 4b) based on the pictures from Micro-CT scanning. We accurately measured the thickness of bovine eyeball wall, including sclera and cornea, using MIMICS software and established a three-dimensional FEM of bovine eyes, in which important factors for FE model such as geometry and anatomical structure of eyeballs, scleral and corneal were considered.

Fig. 4
figure 4

Micro-CT scanning image of bovine eyes and the established FEM model. a: Image of Micro-CT scanning of bovine eyes and scleral thickness measurement; b: Three dimensional FEM based on bovine eyes

Methods optimization

Most previous studies used linear elastic material when selecting the material parameter of the eyeball model [3841], however [42], the human body is actually made up of both hyperelastic and nonlinear viscous materials [42]. It performs as hyperelastic material under low speed conditions and its stress–strain curve changes along with the alteration of loading conditions, with high velocity loadings. In order to reflect these two characteristics, our research combines the nonlinear elastic material theory of Ogden [43] with viscoelastic mechanics theory of Christensen [44]. Ogden’s theory takes α and μ as material parameters and can adequately reflect the non-linear stress–strain curve under Quasi static state; Christensen’s theory takes β and G as material parameters and can reflect the mechanical properties that stress and strain change with the alteration of loading velocity, and the mathematic derivation process is as follows:

Based on Ogden’ theory, the energy density equation for incompressible materials can be expressed as follows:

$$ W={\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=1}^n\frac{\mu_j}{\alpha_j}\left({\lambda}_i^{\alpha_j}-1\right)+\frac{1}{2}K{\left(J-1\right)}^2={\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=1}^n\frac{\mu_j}{\alpha_j}\left({\lambda}_i^{\alpha_j}-1\right)}}}} $$

Where I set as the direction of deformation, J the order for parameter calculation and n can be any integer value from 1 to 8. We only represent three perpendicular directions in this equation.

$$ W={\displaystyle \sum_{i=1}^3{\displaystyle \sum_{j=1}^n\frac{\mu_j}{\alpha_j}\left({\lambda}_i^{\alpha_j}-1\right)={\displaystyle \sum_{j=1}^n\frac{\mu_j}{\alpha_j}\left({\lambda}_x^{\alpha_j}+{\lambda}_y^{\alpha_j}+{\lambda}_z^{\alpha_j}-3\right)}}} $$

The following can be derived from equation (1) and (2):

$$ W={\displaystyle \sum_{j=1}^n\frac{\mu_j}{\alpha_j}\left({\lambda}_x^{\alpha_j}+2{\lambda}_x^{\frac{-1}{2}{\alpha}_j}-3\right)} $$

Therefore the hyperelastic stress can be expressed as:

$$ {\sigma}_{hyper}={\mu}_j\left({\lambda}_x^{\alpha_j}-{\lambda}_x^{\frac{-{\alpha}_j}{2}}\right) $$

Based on the viscoelastic theory of Christensen, the effect of strain rate on the viscoelastic stress can be expressed by a nonlinear viscoelastic all-integral equation:

$$ {\sigma}_{ij}={\displaystyle {\int}_0^t{g}_{ijkl}}\left(t-\tau \right)\frac{\partial {\varepsilon}_{kj}}{\partial \tau }d\tau $$

The total stress on the tested material is combined from both hyperelastic part and viscoelastic part, as shown in equation (6):

$$ {\sigma}_{total}={\mu}^{*}{\lambda}_x^{\alpha -1}+{\displaystyle \underset{0}{\overset{t}{\int }}{G}^{*}{e}^{-\beta \left(t-\tau \right)}}\frac{\partial \varepsilon }{\partial \tau }d\tau $$

Perform the reverse engineering of material mechanics on the experimental samples. Firstly, establish the mechanical model of experimental samples (Fig. 5):

Fig. 5
figure 5

FEM of sclera and cornea. a: sclera; b: cornea

Simulate corneal stress based on results of tension experiments and perform the curve-fitting of displacement versus force curve. Obtain hyperelastic-viscoelastic mechanical parameters of materials by using optimal calculation of reverse engineering methods showed as follows. Finish the optimal reverse engineering of scleral and corneal and obtain the materials parameters using Non-dominated Sorting Genetic Algorithm-II (NSGA-II) [45], set α, μ, β, Gas input variables, μ and G stands for hyperelastic modulus and viscoelastic modulus, respectively, measured in units of MPa. β stands for hyperelastic parameter and viscoelastic parameter firstly.

Secondly, finish the design of experiment (DOE) and define the scope of each input variable (Fig. 6):

Fig. 6
figure 6

NSGA-II optimiazation Algorithm process design

Use Mode FRONTIER software (Version 4.0, ESTECO, Trieste, Italy) to optimize four Ogden material models in FEM: parameter α, μ, β and G. Steps are as follows: The upper and lower limits of parameter α, μ, β and G are set initial values according to values recorded in pertinent literature, use Mode FRONTIER to randomly choose parameter values within its ranges as input value, calculate with formula used in theoretical analysis of continuum medium mechanics, take the calculation results with minimum error value as the optimal results. When designing the experiment, use Sobol’s method to randomly select a series of samples as estimated input variables and establish the sample space. Samples were uniformly distributed in the generation region. Maximum Number of Iteration is 100. Repeatedly calls the simulative running of LS-DYNA (LS-DYNA3D 971, LSTC, Livermore, CA, USA) FEM, the number of call can be adjusted according to error range. After finishing the optimization, enter the best estimating parameters as the material parameter into the FEM and perform analog computation using LS-DYNA.

Mechanical condition settings of IOP

In order to fully simulate the impact of IOP on the inner wall of eyeball, we chose nodes on the inner wall of eyeball as the setting area that bear stress; the direction of IOP was perpendicular to each node. Perform analog computation after setting pressure value at 10 mmHg, 30 mmHg, 60 mmHg and 100 mmHg, where 10 mmHg is the IOP value for normal human eyes; 30 mmHg is the IOP value for low or middle IOP glaucoma [46]; 60 mmHg is the IOP value for HIOP glaucoma [47]; 100 mmHg is the commonly used IOP value for animal glaucoma experiments [48].


LC size and corneal and scleral thickness

It is necessary to measure the LC size of bovine eyes since the Micro-CT is not able to show its size and position. Twenty processed fresh bovine eyes have been selected; the LC size and scleral thickness around optic papilla of each sample were measured using a vernier caliper (Table 1). The results are as follows: the naso-temporal length range of LC was: 4.37 to 4.71 mm, the longitudinal length range of LC was: 3.49 to 3.84 mm, range of scleral thickness around optic papilla: 1.55 to 1.86 mm.

Table 1 The measurement of size of LC

Deformation-force curve

Figures 7 and 8 showed the measured displacement versus force curves under different tensile velocity. The deformation-force curve with lowest speed was at quasi-static state in this loading, and the hyper-elastic material parameters α and μ could be obtained by optimal NSGA-II reverse engineering. In order to obtain a more accurate viscoelastic parameter, we performed tension tests at multiple strain rates and the obtained data were used for calculation of viscoelastic material parameter β and G (Table 2). The results showed that the stress was smallest under low speed (0.03 mm/s) condition; the maximum force was 16 N from the obtained experimental curve; when the stretching speed is 300 mm/s, the maximum force reached 30 N. This phenomenon defined the time-dependent viscoelastic properties of scleral material.

Fig. 7
figure 7

Deformation-stress curve of sclera under tensile rate of 0.3 mm/s-300 mm/s. a: 0.3 mm/s; b: 3 mm/s; c: 30 mm/s; d: 300 mm/s

Fig. 8
figure 8

Deformation-stress curve of cornea under tensile rate of 0.3 mm/s-300 mm/s. a: 0.3 mm/s; b: 30 mm/s

Table 2 Hyperelastic parameters and viscoelastic parameters after optimization

Hyperelastic parameters and viscoelastic parameters after optimization

The reverse optimization method of FEM that we used was based on optimization method simulated by small FEM. Compared with the traditional curve fitting methods, this method performed specimen-specific optimization therefore can more accurately obtain complicated nonlinear hyperelastic parameters and viscoelastic parameters at the same time. NSGA-II optimization method was widely adopted by researchers in the field of FEM optimization [49] and has made significant contributions in improving the calculation accuracy in biomechanical field [50]. By applying this method to the establishment of bovine eye FEM, we not only obtained more accurate material parameters, as shown in Table 2, but also established a more realistic eyeball model and provided a reliable basis to truly simulate the biomechanical responses when IOP was elevated.

FEM-simulated IOP change

According to the principle of continuum mechanics, we simulated the morphologic strain and mechanical stress of cornea and sclera under different IOP, and obtained the following equivalent strain and stress contours of bovine eye. Figures 9, 10, 11, 12 were the equivalent strain nephogram, while Figs. 13, 14, 15, 16 were the equivalent stress contours. Strain or stress value was increased as the color turned redder, and got closer to zero as color became bluer. By analyzing equivalent stress contours, we found that as simulated IOP was increased, the blue area in the posterior half of sclera expanded while the green area narrowed, indicating that as IOP increased, the stress change in the anterior half of sclera was not obvious while the stress in the posterior half was lowered. By analyzing the corneal equivalent stress nephogram (Fig. 15) we found that no matter how we changed the IOP, the stress of the cornea was always most concentrated in the central part. Moreover, by analyzing equivalent stress contours of LC (Fig. 16), we found that LC stress was mainly concentrated in the central part as IOP was increased.

Fig. 9
figure 9

Equivalent strain nephogram of eyeball. a: 10 mmHg; b: 30 mmHg: c: 60 mmHg; d: 100 mmHg

Fig. 10
figure 10

Equivalent strain nephogram of sclera. a: 10 mmHg; b: 30 mmHg; c: 60 mmHg; d: 100 mmHg

Fig. 11
figure 11

Equivalent strain nephogram of cornea. a: 10 mmHg; b: 30 mmHg; c: 60 mmHg; d: 100 mmHg

Fig. 12
figure 12

Equivalent strain nephogram of LC. a: 10 mmHg; b: 30 mmHg; c: 60 mmHg; d: 100 mmHg

Fig. 13
figure 13

Equivalent stress nephogram of eyeball. a: 10 mmHg; b: 30 mmHg; c: 60 mmHg; d: 100 mmHg

Fig. 14
figure 14

Equivalent stress nephogram of sclera. a: 10 mmHg; b: 30 mmHg; c: 60 mmHg; d: 100 mmHg

Fig. 15
figure 15

Equivalent stress nephogram of cornea. a: 10 mmHg; b: 30 mmHg; c: 60 mmHg; d: 100 mmHg

Fig. 16
figure 16

Equivalent stress nephogram of LC. a: 10 mmHg; b: 30 mmHg; c: 60 mmHg; d: 100 mmHg


Acute HIOP glaucoma has some important pathological features such as the damage and atrophy of optic nerve, progressive loss of the RGCs, and so on. It is characterized by the most important clinical symptom-tubular- like visual field. Patients with tubular visual field normally lose peripheral visual field first, which expands to the central field as the disease deteriorates, eventually causing tubular-like visual field. In other words, the visual sensitivity loss and subjective visual sensory defect happens later in the central part of retina than the peripheral part in the progression of glaucoma. Preliminary clinical study showed that RGCs death did not occur at the same time in glaucoma, instead, peripheral RGCs died earlier than the central part [51]. This phenomenon is well known as the selected RGCs death in glaucoma patients [52, 53].

The underlying mechanisms of selected RGCs death in glaucoma have been extensively studied. Some evidences implied that the death of RGC was affected by its size. For example, Quigley et al. found that large RGCs were relatively easy to die in glaucoma [10, 54, 55]. However, this argument is controversial as other researchers found that there was no significant difference between the number of lost small and large RGCs examined by retro-tracing RGCs in experimental glaucoma eyes using horseradish peroxidase [53]. Another potential mechanism of selected RGCs death might be related to the expression of functional tyrosine kinase in retinal ganglion cells [1220], however this hypothesis is contradictory to optic nerve transection model and chronic HIOP model [56, 57]. Other experiments using subunit of glutamate receptor to analyze the selected RGCs death failed to provide promising explanation [58, 59]. Recent studies suggested that neural degeneration in glaucoma was also related to immune system [6062], and they pointed out that the loss of complement activity in both acquired and innate immune process was involved in RGCs damage [63, 64]. However the roles of the complement system in RGCs damage remain to be further studied. Tong et al. studied the selected RGCs death from the perspective of retinal blood supply [21]. Their results showed that the number of lost RGCs in peripheral part was much more than the central part in rats with acute HIOP and its loss rate was obviously related to local blood supply. This result could be explained by the fact that compared with blood vessels in central part, peripheral blood vessel possess a thinner wall and is less elastic, thus prone to developing ischemia by oppression. Although there exist various hypotheses for selected RGCs death, its mechanisms are yet to be further investigated.

Our present experiment established the mathematical FEM of bovine eyeball and studied it from a different perspective. Equivalent stress nephogram of sclera obtained by simulating different IOP (Fig. 14) showed that as simulated IOP increased, the blue area in posterior portion of scleral FEM gradually expanded and in the meantime, the scope of green area gradually narrowed. This result showed that with the gradual increase of IOP, the stress change of anterior half of sclera, which was in front of middle axis of eyeball and approximately represented the areas around retina, was not obvious, while stress on the posterior half of sclera, which approximately represented the middle and central area of retina, becomes smaller. That means, with the increase of simulated IOP, the stress is mainly concentrated in the peripheral part rather than the central or middle part, suggesting that important pressure-bearing tissues such as sclera will deform as IOP increases, leading to corresponding changes of local stress. Such stress can indirectly reflect the pressure from inner side of wall of eyeball. These above results are consistent with Tong's study by establishing mathematical models and enrich his theory by further discussing reasons for selected RGCs death from the perspective of local stress rather than being limited to vascular factors. Another interesting finding from the simulation of increasing IOP by FEM is that, as IOP increased, most stress is concentrated in central part of LC (Fig. 16). This result is similar to Sigal’s [65], showing that optic nerve passing through central part of LC bears the maximum pressure, which is much higher than those taken by peripheral optic nerve. We believe that the scientific significance behind this interesting phenomenon is worth further study.

The results of our FEM simulation of increased IOP bring up another problem worthy of attention. Regardless of how we change the simulative IOP, corneal stress is always most concentrated in the central part (Fig. 15), suggesting that IOP in the central area of cornea is the highest. Clinical researches found that after cutting cornea by laser, as is used in laser-assisted in situ Keratomileusis operation for myopia, patients are prone to develop secondary keratoconus [66]. It is generally believed that keratoconus is associated with the thinning of central portion of cornea, which can be verified by our FEM since corneal stress is most concentrated in the central part.

It should be noted that there are some limitations in our study. Firstly, we only studied the in vitro response of the eyeball materials, more in vivo biomechanical experiments of eyeball will be studied in our future work. Secondly, the prediction of HIOP in FE modeling simulation need to be validated in in vivo experiment, so the clinical value of this FE modeling will be greatly improved.


As simulative IOP is increased, the stress-concentration scope of the posterior half of sclera gradually narrows, at the same time, the stress-concentration scope on the anterior half of sclera gradually expands, and stress on the LC is mainly concentrated in the central part, suggesting that IOP becomes mainly concentrated in the anterior part of the eyeball as it increases. These findings provide a biomechanical evidence to explain why RGCs in peripheral part die earlier than RGCs in central part under HIOP.


  1. Huang JF, Shang L, Zhang MQ, Wang H, Chen D, Tong JB, et al. Differential neuronal expression of receptor interacting protein 3 in rat retina: involvement in ischemic stress response. BMC Neurosci. 2013;14:16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Jiang SH, Shang L, Xue LX, Ding W, Chen S, Ma RF, et al. The effect and underlying mechanism of Timosaponin B-II on RGC-5 necroptosis induced by hydrogen peroxide. BMC Complementary Altern Med. 2014;14:459.

    Article  Google Scholar 

  3. Shang L, Huang JF, Ding W, Chen S, Xue LX, Ma RF, et al. Calpain: a molecule to induce AIF-mediated necroptosis in RGC-5 following elevated hydrostatic pressure. BMC Neurosci. 2014;15:63.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Quigley HA, Vitale S. Models of open-angle glaucoma prevalence and incidence in the United States. Invest Ophthalmol Vis Sci. 1997;38(1):83–91.

    CAS  PubMed  Google Scholar 

  5. Quigley HA, Broman AT. The number of people with glaucoma worldwide in 2010 and 2020. Br J Ophthalmol. 2006;90(3):262–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Gupta N, Yucel YH. Glaucoma as a neurodegenerative disease. Curr Opin Ophthalmol. 2007;18(2):110–4.

    Article  PubMed  Google Scholar 

  7. Ying.Tian: Reviews on early diagnosis therapy of primary glaucoma. Chinese Journal of Practical Ophthalmology 2002(01):3–7.

  8. Li FM. Chinese Ophthalmology. Beijing: People’s Medical Publishing House; 2004.

    Google Scholar 

  9. Quigley HA, Broman AT. The number of people with glaucoma worldwide in 2010 and 2020. Brit J Ophthalmol. 2006;90(3):262–7.

    Article  CAS  Google Scholar 

  10. Quigley HA. Neuronal death in glaucoma. Prog Retin Eye Res. 1999;18(1):39–57.

    Article  CAS  PubMed  Google Scholar 

  11. Quigley HA. Open-angle glaucoma. N Engl J Med. 1993;328(15):1097–106.

    Article  CAS  PubMed  Google Scholar 

  12. Kuehn MH, Fingert JH, Kwon YH. Retinal ganglion cell death in glaucoma: mechanisms and neuroprotective strategies. Ophthalmol Clin North Am. 2005;18(3):383–95. vi.

    Article  PubMed  Google Scholar 

  13. Quigley HA, McKinnon SJ, Zack DJ, Pease ME, Kerrigan-Baumrind LA, Kerrigan DF, et al. Retrograde axonal transport of BDNF in retinal ganglion cells is blocked by acute IOP elevation in rats. Invest Ophthalmol Vis Sci. 2000;41(11):3460–6.

    CAS  PubMed  Google Scholar 

  14. Pease ME, McKinnon SJ, Quigley HA, Kerrigan-Baumrind LA, Zack DJ. Obstructed axonal transport of BDNF and its receptor TrkB in experimental glaucoma. Invest Ophthalmol Vis Sci. 2000;41(3):764–74.

    CAS  PubMed  Google Scholar 

  15. Ko ML, Hu DN, Ritch R, Sharma SC, Chen CF. Patterns of retinal ganglion cell survival after brain-derived neurotrophic factor administration in hypertensive eyes of rats. Neurosci Lett. 2001;305(2):139–42.

    Article  CAS  PubMed  Google Scholar 

  16. Martin KRG, Quigley HA, Zack DJ, Levkovitch-Verbin H, Kielczewski J, Valenta D, et al. Gene therapy with brain-derived neurotrophic factor as a protection: Retinal ganglion cells in a rat glaucoma model. Invest Ophthalmol Vis Sci. 2003;44(10):4357–65.

    Article  PubMed  Google Scholar 

  17. Ji JZ, Elyaman W, Yip HK, Lee VWH, Yick LW, Hugon J, et al. CNTF promotes survival of retinal ganglion cells after induction of ocular hypertension in rats: the possible involvement of STAT3 pathway. Eur J Neurosci. 2004;19(2):265–72.

    Article  PubMed  Google Scholar 

  18. Van Adel BA, Kostic C, Deglon N, Ball AK, Arsenijevic Y. Delivery of ciliary neurotrophic factor via lentiviral-mediated transfer protects axotomized retinal ganglion cells for an extended period of time. Hum Gene Ther. 2003;14(2):103–15.

    Article  PubMed  Google Scholar 

  19. Wang H, Liu QL, Luo XG, Wen JY. Protection of the retinal ganglion cells from experimental glaucoma by brain derived neurotrophic factor. Anatomy Research. 2002;24(2):119–22. +165.

    CAS  Google Scholar 

  20. Huang JF, Jiang LZ, Tong JB, Chen D, Xiong K, Zeng LP: Expression of TrkB in rat retina following acute high intraocular pressure induced by brain-derived neurotrophic factor. Anatomy Research 2006(06):734-737 + 740.

  21. Tong JB, Chen D, Zeng LP, Mo XY, Wang H, Huang JF, et al. Differential Changes of Local Blood Supply in Rat Retinae Are Involved in the Selective Loss of Retinal Ganglion Cells Following the Acute High Intraocular Pressure. Curr Eye Res. 2010;35(5):425–34.

    Article  PubMed  Google Scholar 

  22. St Helen R, Mc EW. Rheology of the human sclera. 1. Anelastic behavior. Am J Ophthalmol. 1961;52:539–48.

    Article  CAS  PubMed  Google Scholar 

  23. Arciniegas A, Amaya LE. Mechanical behavior of the sclera. Ophthalmologica. 1986;193(1–2):45–55.

    Article  CAS  PubMed  Google Scholar 

  24. Eberhardt AW, Grissom JD, Siegwart Jr JT, Norton TT: Analytical and finite element models of scleral creep in tree shrews with deprivation-induced myopia. In: The 1996 3 rd Biennial Joint Conference on Engineering Systems Design and Analysis, ESDA Part 5(of 9): 1996; 1996: 189–193.

  25. Phillips JR, Khalaj M, McBrien NA. Induced myopia associated with increased scleral creep in chick and tree shrew eyes. Invest Ophthalmol Vis Sci. 2000;41(8):2028–34.

    CAS  PubMed  Google Scholar 

  26. Bellezza AJ, Hart RT, Burgoyne CF. The optic nerve head as a biomechanical structure: initial finite element modeling. Investig Ophthalmol Vis Sci. 2000;41(10):2991–3000.

    CAS  Google Scholar 

  27. Sigal IA, Flanagan JG, Tertinegg I, Ethier CR. Reconstruction of human optic nerve heads for finite element modeling. Technology and health. 2005;13(4):313–29.

    Google Scholar 

  28. Norman RE, Flanagan JG, Sigal IA, Rausch SM, Tertinegg I, Ethier CR. Finite element modeling of the human sclera: influence on optic nerve head biomechanics and connections with glaucoma. Exp Eye Res. 2011;93(1):4–12.

    Article  CAS  PubMed  Google Scholar 

  29. Sigal IA, Flanagan JG, Tertinegg I, Ethier CR. Finite element modeling of optic nerve head biomechanics. Invest Ophthalmol Vis Sci. 2004;45(12):4378–87.

    Article  PubMed  Google Scholar 

  30. Burgoyne CF, Downs JC, Bellezza AJ, Suh JK, Hart RT. The optic nerve head as a biomechanical structure: a new paradigm for understanding the role of IOP-related stress and strain in the pathophysiology of glaucomatous optic nerve head damage. Prog Retin Eye Res. 2005;24(1):39–73.

    Article  PubMed  Google Scholar 

  31. Quigley HA, Flower RW, Addicks EM, McLeod DS. The mechanism of optic nerve damage in experimental acute intraocular pressure elevation. Invest Ophthalmol Vis Sci. 1980;19(5):505–17.

    CAS  PubMed  Google Scholar 

  32. Quigley HA, Addicks EM, Green WR, Maumenee AE. Optic nerve damage in human glaucoma. II. The site of injury and susceptibility to damage. Arch Ophthalmol. 1981;99(4):635–49.

    Article  CAS  PubMed  Google Scholar 

  33. Quigley HA, Hohman RM, Addicks EM, Massof RW, Green WR. Morphologic changes in the lamina cribrosa correlated with neural loss in open-angle glaucoma. Am J Ophthalmol. 1983;95(5):673–91.

    Article  CAS  PubMed  Google Scholar 

  34. Sigal IA, Flanagan JG, Ethier CR. Factors influencing optic nerve head biomechanics. Invest Ophthalmol Vis Sci. 2005;46(11):4189–99.

    Article  PubMed  Google Scholar 

  35. Cahane M, Bartov E. Axial length and scleral thickness effect on susceptibility to glaucomatous damage: a theoretical model implementing Laplace's law. Ophthalmic Res. 1992;24(5):280–4.

    Article  CAS  PubMed  Google Scholar 

  36. Girard MJ, Downs JC, Bottlang M, Burgoyne CF, Suh JK. Peripapillary and posterior scleral mechanics--part II: experimental and inverse finite element characterization. J Biomech Eng. 2009;131(5):051012.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Girard MJ, Downs JC, Burgoyne CF, Suh JK. Peripapillary and posterior scleral mechanics--part I: development of an anisotropic hyperelastic constitutive model. J Biomech Eng. 2009;131(5):051011.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Uchio E, Ohno S, Kudoh J, Aoki K, Kisielewicz LT. Simulation model of an eyeball based on finite element analysis on a supercomputer. Br J Ophthalmol. 1999;83(10):1106–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Schutte S, van den Bedem SP, van Keulen F, van der Helm FC, Simonsz HJ. A finite-element analysis model of orbital biomechanics. Vision Res. 2006;46(11):1724–31.

    Article  PubMed  Google Scholar 

  40. Rangarajan N, Kamalakkannan SB, Hasija V, Shams T, Jenny C, Serbanescu I, et al. Finite element model of ocular injury in abusive head trauma. J AAPOS. 2009;13(4):364–9.

    Article  PubMed  Google Scholar 

  41. Asejczyk-Widlicka M, Śródka D, Kasprzak H, Iskander D. Influence of intraocular pressure on geometrical properties of a linear model of the eyeball: Effect of optical self-adjustment. Optik-International Journal for Light and Electron Optics. 2004;115(11):517–24.

    Article  Google Scholar 

  42. Feng YZ. biomechanics. Beijing: Science Press; 1983.

    Google Scholar 

  43. Ogden R. Large deformation isotropic elasticity-on the correlation of theory and experiment for incompressible rubberlike solids. Proceedings of the Royal Society of London A Mathematical and Physical Sciences. 1972;326(1567):565–84.

    Article  CAS  Google Scholar 

  44. Christensen RM. Theory of viscoelasticity. 2nd ed. Mineola, N.Y.: Dover Publications; 2003.

    Google Scholar 

  45. Deb K, Pratap A, Agarwal S, Meyarivan T. A fast and elitist multiobjective genetic algorithm: NSGA-II. Evolutionary Computation, IEEE Transactions on. 2002;6(2):182–97.

    Article  Google Scholar 

  46. Agar A, Li S, Agarwal N, Coroneo MT, Hill MA. Retinal ganglion cell line apoptosis induced by hydrostatic pressure. Brain Res. 2006;1086:191–200.

    Article  CAS  PubMed  Google Scholar 

  47. Tan SZ, Sampat K, Rasool S, Nolan D: Unilateral acute angle closure glaucoma. BMJ case reports 2013, 2013:bcr2012007836

  48. Huang JF, Zhou LH, Wang H, Luo J, Zeng LP, Xiong K, et al. Distribution of thrombospondins and their neuronal receptor alpha 2 beta 1 in the rat retina. Exp Eye Res. 2013;111:36–49.

    Article  CAS  PubMed  Google Scholar 

  49. Guan FJ, Han X, Mao HJ, Wagner C, Yeni YN, Yang KH. Application of Optimization Methodology and Specimen-Specific Finite Element Models for Investigating Material Properties of Rat Skull. Ann Biomed Eng. 2011;39(1):85–95.

    Article  PubMed  Google Scholar 

  50. Hu JW, Klinich KD, Miller CS, Nazmi G, Pearlman MD, Schneider LW, et al. Quantifying dynamic mechanical properties of human placenta tissue using optimization techniques with specimen-specific finite-element models. J Biomech. 2009;42(15):2528–34.

    Article  PubMed  Google Scholar 

  51. Stamper RL. The effect of glaucoma on central visual function. Trans Am Ophthalmol Soc. 1984;82:792–826.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. Laquis S, Chaudhary P, Sharma SC. The patterns of retinal ganglion cell death in hypertensive eyes. Brain Res. 1998;784(1–2):100–4.

    Article  CAS  PubMed  Google Scholar 

  53. Urcola JH, Hernandez M, Vecino E. Three experimental glaucoma models in rats: comparison of the effects of intraocular pressure elevation on retinal ganglion cell size and death. Exp Eye Res. 2006;83(2):429–37.

    Article  CAS  PubMed  Google Scholar 

  54. Quigley HA. Ganglion cell death in glaucoma: pathology recapitulates ontogeny. Aust N Z J Ophthalmol. 1995;23(2):85–91.

    Article  CAS  PubMed  Google Scholar 

  55. Glovinsky Y, Quigley HA, Pease ME. Foveal ganglion cell loss is size dependent in experimental glaucoma. Invest Ophthalmol Vis Sci. 1993;34(2):395–400.

    CAS  PubMed  Google Scholar 

  56. Cui Q, Tang LS, Hu B, So KF, Yip HK. Expression of trkA, trkB, and trkC in injured and regenerating retinal ganglion cells of adult rats. Invest Ophthalmol Vis Sci. 2002;43(6):1954–64.

    PubMed  Google Scholar 

  57. Rudzinski M, Wong TP, Saragovi HU. Changes in retinal expression of neurotrophins and neurotrophin receptors induced by ocular hypertension. J Neurobiol. 2004;58(3):341–54.

    Article  CAS  PubMed  Google Scholar 

  58. Hof PR, Lee PY, Yeung G, Wang RF, Podos SM, Morrison JH. Glutamate receptor subunit GluR2 and NMDAR1 immunoreactivity in the retina of macaque monkeys with experimental glaucoma does not identify vulnerable neurons. Exp Neurol. 1998;153(2):234–41.

    Article  CAS  PubMed  Google Scholar 

  59. Dijk F, Kraal-Muller E, Kamphuis W. Ischemia-induced changes of AMPA-type glutamate receptor subunit expression pattern in the rat retina: a real-time quantitative PCR study. Invest Ophthalmol Vis Sci. 2004;45(1):330–41.

    Article  PubMed  Google Scholar 

  60. Tezel G, Wax MB. Glaucoma. Chem Immunol Allergy. 2007;92:221–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Wax MB, Tezel G. Immunoregulation of retinal ganglion cell fate in glaucoma. Exp Eye Res. 2009;88(4):825–30.

    Article  CAS  PubMed  Google Scholar 

  62. Tezel G. The role of glia, mitochondria, and the immune system in glaucoma. Invest Ophthalmol Vis Sci. 2009;50(3):1001–12.

    Article  PubMed  Google Scholar 

  63. Wax MB, Tezel G, Yang J, Peng G, Patil RV, Agarwal N, et al. Induced autoimmunity to heat shock proteins elicits glaucomatous loss of retinal ganglion cell neurons via activated T-cell-derived fas-ligand. J Neurosci. 2008;28(46):12085–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Tezel G, Hernandez R, Wax MB. Immunostaining of heat shock proteins in the retina and optic nerve head of normal and glaucomatous eyes. Arch Ophthalmol. 2000;118(4):511–8.

    Article  CAS  PubMed  Google Scholar 

  65. Sigal IA, Flanagan JG, Tertinegg I, Ethier CR. Modeling individual-specific human optic nerve head biomechanics. Part I: IOP-induced deformations and influence of geometry. Biomech Model Mechanobiol. 2009;8(2):85–98.

    Article  PubMed  Google Scholar 

  66. Gefen A, Shalom R, Elad D, Mandel Y. Biomechanical analysis of the keratoconic cornea. J Mech Behav Biomed Mater. 2009;2(3):224–36.

    Article  PubMed  Google Scholar 

Download references


This study was supported by National Natural Science Foundation of China (81371011); National Key Technologies Research and Development Program of China (2012BAK14B03); Wu Jie-Ping Medical Foundation of the Minister of Health of China (320.6750.14118); Natural Science Foundation of Hunan Province (2015JJ2187) and Teacher Research Foundation of Central South University (2014JSJJ026).

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Na Li or Kun Xiong.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

KX, NL and JFH designed the experiment, YHC, WW and NL performed the experiment, NL and SL analyzed the data, YHC and SYC drafted the manuscript, KX and NL participated in paper modification and revised the manuscript for English writing, all authors participated in critical revision of the manuscript and approved the final manuscript. All authors read and approved the final manuscript.

Yan-Hui Cui and Ju-Fang Huang contributed equally to this work.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cui, YH., Huang, JF., Cheng, SY. et al. Study on establishment and mechanics application of finite element model of bovine eye. BMC Ophthalmol 15, 101 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: