This article has Open Peer Review reports available.
Study on establishment and mechanics application of finite element model of bovine eye
- Yan-Hui Cui†1,
- Ju-Fang Huang†1,
- Si-Ying Cheng2,
- Wei Wei3,
- Lei Shang1,
- Na Li4Email author and
- Kun Xiong1Email author
© Cui et al. 2015
Received: 14 November 2014
Accepted: 14 July 2015
Published: 13 August 2015
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.
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) [1–3]. 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. [4–8]. 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 . 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 [12–20], while others showed that the blood supply was different in different parts of retina . 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 [22–24]. 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 , 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.  studied the IOP-induced biomechanical changes of sclera , and Sigal et al.  proved that FE model had a good consistency with real eyes in geometric shape, elastic mechanical response and viscous mechanical response. Richard et al.  established a “Universal” FEM of human optic papilla which contained LC using the method by Sigal et al. . 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 . Although the mechanisms have not been completely investigated, researches have found that LC was a vulnerable area to injury [31–33], 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.
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
Measurement of bovine LC
Preparation of biomechanical experiment
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.
Most previous studies used linear elastic material when selecting the material parameter of the eyeball model [38–41], however , the human body is actually made up of both hyperelastic and nonlinear viscous materials . 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  with viscoelastic mechanics theory of Christensen . 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:
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) , 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.
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 ; 60 mmHg is the IOP value for HIOP glaucoma ; 100 mmHg is the commonly used IOP value for animal glaucoma experiments .
LC size and corneal and scleral thickness
The measurement of size of LC
Naso-temporal width of LC (mm)
Longitudinal length of LC (mm)
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  and has made significant contributions in improving the calculation accuracy in biomechanical field . 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
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 . 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 . Another potential mechanism of selected RGCs death might be related to the expression of functional tyrosine kinase in retinal ganglion cells [12–20], 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 [60–62], 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 . 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 , 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 . 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.
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).
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.PubMedGoogle Scholar
- Quigley HA, Broman AT. The number of people with glaucoma worldwide in 2010 and 2020. Br J Ophthalmol. 2006;90(3):262–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Gupta N, Yucel YH. Glaucoma as a neurodegenerative disease. Curr Opin Ophthalmol. 2007;18(2):110–4.View ArticlePubMedGoogle Scholar
- Ying.Tian: Reviews on early diagnosis therapy of primary glaucoma. Chinese Journal of Practical Ophthalmology 2002(01):3–7.Google Scholar
- Li FM. Chinese Ophthalmology. Beijing: People’s Medical Publishing House; 2004.Google Scholar
- Quigley HA, Broman AT. The number of people with glaucoma worldwide in 2010 and 2020. Brit J Ophthalmol. 2006;90(3):262–7.View ArticleGoogle Scholar
- Quigley HA. Neuronal death in glaucoma. Prog Retin Eye Res. 1999;18(1):39–57.View ArticlePubMedGoogle Scholar
- Quigley HA. Open-angle glaucoma. N Engl J Med. 1993;328(15):1097–106.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- St Helen R, Mc EW. Rheology of the human sclera. 1. Anelastic behavior. Am J Ophthalmol. 1961;52:539–48.View ArticlePubMedGoogle Scholar
- Arciniegas A, Amaya LE. Mechanical behavior of the sclera. Ophthalmologica. 1986;193(1–2):45–55.View ArticlePubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMedGoogle Scholar
- 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.Google Scholar
- 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
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Sigal IA, Flanagan JG, Ethier CR. Factors influencing optic nerve head biomechanics. Invest Ophthalmol Vis Sci. 2005;46(11):4189–99.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Feng YZ. biomechanics. Beijing: Science Press; 1983.Google Scholar
- 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.View ArticleGoogle Scholar
- Christensen RM. Theory of viscoelasticity. 2nd ed. Mineola, N.Y.: Dover Publications; 2003.Google Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Tan SZ, Sampat K, Rasool S, Nolan D: Unilateral acute angle closure glaucoma. BMJ case reports 2013, 2013:bcr2012007836Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Stamper RL. The effect of glaucoma on central visual function. Trans Am Ophthalmol Soc. 1984;82:792–826.PubMedPubMed CentralGoogle Scholar
- Laquis S, Chaudhary P, Sharma SC. The patterns of retinal ganglion cell death in hypertensive eyes. Brain Res. 1998;784(1–2):100–4.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Quigley HA. Ganglion cell death in glaucoma: pathology recapitulates ontogeny. Aust N Z J Ophthalmol. 1995;23(2):85–91.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Tezel G, Wax MB. Glaucoma. Chem Immunol Allergy. 2007;92:221–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Wax MB, Tezel G. Immunoregulation of retinal ganglion cell fate in glaucoma. Exp Eye Res. 2009;88(4):825–30.View ArticlePubMedGoogle Scholar
- Tezel G. The role of glia, mitochondria, and the immune system in glaucoma. Invest Ophthalmol Vis Sci. 2009;50(3):1001–12.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Gefen A, Shalom R, Elad D, Mandel Y. Biomechanical analysis of the keratoconic cornea. J Mech Behav Biomed Mater. 2009;2(3):224–36.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.