Ex vivo testing of intact eye globes under inflation conditions to determine regional variation of mechanical stiffness
© The Author(s). 2016
Received: 27 September 2015
Accepted: 23 July 2016
Published: 10 August 2016
The eye globe exhibits significant regional variation of mechanical behaviour. The aim of this present study is to develop a new experimental technique for testing intact eye globes in a form that is representative of in vivo conditions, and therefore suitable for determining the material properties of the complete outer ocular tunic.
A test rig has been developed to provide closed-loop control of either applied intra-ocular pressure or resulting apical displacement; measurement of displacements across the external surface of the eye globe using high-resolution digital cameras and digital image correlation software; prevention of rigid-body motion and protection of the ocular surface from environmental drying. The method has been demonstrated on one human and one porcine eye globe, which were cyclically loaded. Finite element models based on specimen specific tomography, free from rotational symmetry, were used along with experimental pressure-displacement data in an inverse analysis process to derive the mechanical properties of tissue in different regions of the eye’s outer tunic.
The test method enabled monitoring of mechanical response to intraocular pressure variation across the surface of the eye globe. For the two eyes tested, the method showed a gradual change in the sclera’s stiffness from a maximum at the limbus to a minimum at the posterior pole, while in the cornea the stiffness was highest at the centre and lowest in the peripheral zone. Further, for both the sclera and cornea, the load–displacement behaviour did not vary significantly between loading cycles.
The first methodology capable of mechanically testing intact eye globes, with applied loads and boundary conditions that closely represent in vivo conditions is introduced. The method enables determination of the regional variation in mechanical behaviour across the ocular surface.
KeywordsOcular biomechanics Experimental testing Digital image correlation
Simulating the response of the ocular vessel to external and internal forces is of particular interest as it can lead to developments in diagnostic and corrective procedures for various degenerative visual conditions such as glaucoma, keratoconus and myopia [1–6]. Accurate simulations require knowledge of the tissue’s material properties, geometry, loading and boundary conditions, in addition to the numerical tools, with which to describe these details. Simulations can be carried out using numerical tools in the form of finite element models, and the material properties described within these models take the form of constitutive numerical relationships, which can encompass multiple facets of material behaviour.
The multi-faceted behaviour of corneoscleral tissue can be described as anisotropic, hyperelastic [7, 8] and viscoelastic [9, 10]. The stroma is the primary contributor to this behaviour and in the human cornea constitutes 90 % of its thickness. It consists of stacked lamellae of collagen fibrils embedded in a hydrated matrix of proteoglycans, glycoproteins, and keratocytes [11–13]. Similarly, the sclera consists predominantly of collagen fibrils [8, 13] combined as either interlaced bundles or defined lamellae. The fibrils are heterotypic structures as they combine high proportions of collagen types I and III with lower proportions of types V and VI [14, 15]. The straightening of collagen fibrils under stress [16, 17] is a main cause of hyperelastic behaviour where gradual stiffening of tissue under stress continues until most fibres have become taut and able to dominate the behaviour with a linear response to loading. On the other hand, viscoelastic behaviour results from the large fluid component of corneoscleral tissue, which introduces a long-term response to loading and strain-rate dependency. These two observations illustrate the dependence of behaviour on the ocular microstructure, which means that reconfiguration of the microstructure (or re-alignment of the collagen fibres), such as that due to surgical manipulation or changes in the anisotropic strain distribution [9, 18, 19], should be avoided in experimental testing, otherwise the mechanical behaviour would not be relevant to the dominant in vivo conditions.
Several biomechanical studies have been conducted on ex vivo ocular tissues employing mainly uniaxial tension tests [20–26] and inflation tests on separated corneas [27–30] and part scleras [31–33]. Comparative studies of the mechanical stiffness (as measured by the tangent modulus) obtained from uniaxial and inflation tests observed a significant increase in values obtained from the former tests [34, 35]. The higher stiffness estimation was attributed in part to the non-physiologic loading conditions and the higher strain rates commonly employed in uniaxial tension tests . Further, inflation tests of separated corneas or part scleras have a number of drawbacks, namely that (1) the behaviour trends obtained for the corneas and scleras cannot be correlated since the specimens rarely come from the same donor, (2) the anterior sclera is not characterised since it is usually the clamp site for both cornea and sclera specimens, and (3) providing the cornea and sclera specimens with rigid edges during the test creates edge conditions and local stress distributions that are not physiological.
In this study, a case is made for applying the inflation concept on intact eye globes to avoid the shortfalls of both the uniaxial testing of tissue strips and the inflation testing of corneal buttons or scleral cups. While globe inflation testing has been attempted before, the studies were limited to observations of parts of the eye globe; such as posterior sclera , anterior segment [37, 38] or around a single horizon . Further, earlier studies have employed boundary conditions that were dissimilar to in vivo conditions and may therefore have affected to the globe’s response to internal pressure loading.
The present study extends behaviour observations to the whole ocular surface, attempts to use more physiologic boundary conditions and illustrates the new method through an initial experimental test programme involving cyclic inflation of intact eye globes. Attention has been given, as much as possible, to ensure consistency with the in vivo conditions of the eye in order to avoid the main sources of error affecting the reliability of the target material properties of the ocular tissue.
In addition to the experimental testing method, the procedure includes a numerical inverse analysis method that uses the experimental pressure-deformation data to determine the tissue’s hyperelastic stress–strain properties across the ocular surface.
A 2 × 50 mm hypodermic needle was inserted into the ocular cavity through the posterior pole, passed through the back wall of the glass box and its support system and connected through a pipe network to a motor-driven syringe pump, which provided changes in the applied pressure. Fixed at the same elevation as the centre of the eye was a pressure transducer (FDW 060-K262-01, RDP Electronics, USA) that measured the applied pressure with a 0.1 mmHg resolution.
The deformation of the eye that resulted from changing the internal pressure by the syringe pump was measured using a system of three high-resolution, digital cameras (550D, Canon, Japan) with 100 mm (fixed-focal-length) macro lenses that were mounted on a support frame designed for both access and rigidity. The cameras were positioned around the equator of the eye globe facing perpendicular to the sides of the triangular support box. The camera images were initially used to obtain sufficient topography information to construct the eye-specific FE numerical model of the test specimen used in later inverse analysis. The cameras (Fig. 1) were further used during the test to obtain images of the deformed shape of the eye, which were then analysed using digital image correlation (DIC) software (Istra4D, Dantec Dynamics A/S, Denmark) to quantify the displacement distribution across the surface of the eye.
In addition, a laser displacement sensor (LK-2001, Keyence, UK) measured the displacement at the cornea’s apex with 1 μm resolution. This addition to the test rig helped (1) enable the conduct of the test while controlling the displacement at the apex, and (2) provided a direct measurement of an important displacement parameter that aided the validation of the deformation measurements obtained from analysis of the camera images.
Gelatine support system
The eye specimen was cast into a gelatine support, which was produced from Type-A gelatine flakes derived from porcine skin (Sigma-Aldrich cat. No. G2500). This particular gelatine was chosen from trials for its clarity and low spring stiffness and because it had been reported not to chemically alter collagen based tissue . In order to determine the effect of the gelatine on eye globe deformation, its spring stiffness k’ was measured using an Instron uniaxial machine fitted with a 10 N load cell and found to be linear of the value k ' = F/(δA) = 0.0001 N/mm 3, where F is the applied force, δ the resulting deformation, and A the cross-sectional area of gelatine material.
Finite element analysis (FEA) was conducted using the non-linear solver Abaqus/Standard 6.13 (Dassault Systèmes Simulia Corp., Rhode Island, USA) to determine the effect of this stiffness on the deformation of an inflating eye globe. The topography of the model was described by: a corneal radius (7.8 mm) and shape factor (0.82); a sclera radius (11.5 mm); a central corneal thickness, CCT (545 μm – reported as average value in a number of previous studies ) and a peripheral corneal thickness (695 μm, which was consistent with Gullstrand’s No. 1 schematic eye ). The intraocular void was represented by fluid cavity bound by the internal surface of the ocular vessel to which it provided a uniform hydrostatic pressure that varied in the analysis to inflate the eye globe model from 0–30 mmHg.
Deformation at discrete nodal locations from 0–30 mmHg
Control FEM (displacement)
FEM with simulation of gelatine support (displacement, difference)
201 μm, 0.49 %
69.5 μm, 0.52 %
58.2 μm, −0.05 %
86.3 μm, 0.52 %
Calibration of the laser displacement sensor
The laser displacement device uses an optical triangulation position sensor to determine the distance to a near object by the recorded change in angle between the output beam and the reflected beam. Altering the media the laser travels through consequently alters the refractive angle, which in turn affects the triangulation angle. A linear displacement calibrator with a resolution of 1 μm was used to recalibrate the laser displacement transducer through the glass box and gelatine by directly controlling and changing the distance measured by the laser. The relationship between the known relative distance and the relative distance measurements from the laser provided the calibration factor to correct the laser readings. The resulting root mean square (RMS) error after recalibration was 2.7 μm.
Deformation measurement by digital image correlation
DIC is an established technique that relies on the analysis of successive camera images taken for a test specimen to derive the deformation distribution across the specimen surface . This technique has already been successfully applied, in various forms, to ocular tests [21, 28, 32]. This paper describes the application of two-dimensional DIC to determine the planar deformation of the eye globes. For each camera, we consider the circumferential horizon of the eye to be the planar surface. DIC was performed using the software Istra4D. There were errors in the interpretation of initial images with all image processing techniques. This was particularly apparent from distortion, which usually occurs as a result of the lens geometry and due to diffraction in the glass box and the gelatine. The process of minimising the distortion and removing, to a reasonable degree, its influence on the data is described below.
Discrete speckles of cellulose-based paints were chosen for their stability, rapid solvent evaporation and organic nature. The speckles optimised the DIC output and did not provide another surface layer to the specimens. Uniaxial tests were performed on corneal tissue, both with and without applied media, to ensure no mechanical effects resulted from this process.
Referring back to the physical setup described in Fig. 1, the cameras were located orthogonal to the sides of the glass box. This minimises the distortion, which occurs as light was refracted while passing between the air and the glass and between the glass and the gelatine. The error in this system of measurement was quantified using the same setup as used to calibrate the laser beam. Multiple known positions within the gelatine were measured through DIC, and the resulting RMS error between known positions and their associated measurement was 1.7 μm.
Accurate measurements of deformation and loading, and accurate representation of boundary conditions are required to derive representative material properties from an inverse analysis process. Additionally, the reliability of the derived material properties is reliant on an accurate representation of the geometry of the specimen within the finite element simulation. This section describes the method used to recreate the geometry of the specimen in the form of an FEM and the method used to quantify and minimise errors.
where N is the number of pressure levels, M = 26 is the number of measurement locations and δ i,j is the deformation at each particular pressure, at each location. As the deformation was measured relative to the geometry at IOP = 2 mmHg, FEA was performed over two stages. The first stage elevated IOP to match that of the initial pressure of the experiment, and the second stage increased IOP to the final pressure. The deformation from the FEA (δ i numerical ) was recorded relative to the beginning of the second stage.
where λ k ' α are the deviatoric principal stretches equal to J − 1/3 × λ k (k = x, y, z); J = λ x ⋅ λ y ⋅ λ z where λ x , λ y , λ z are the principal stretches in the three main Cartesian directions. The material parameters denoting the strain hardening exponent and the shear modulus are α i and μ i (i = 1 … N) respectively, where N is the function order. The product of stretch in all three directions, J = λ x ⋅ λ y ⋅ λ z = 1, following the approximation that ocular tissue is an almost incompressible material [46, 47]. The values of material parameters α i and \( \mu \) i represented the output of the inverse modelling process described above. The use of a first order material model, N = 1, reduced computation time during the modelling procedure by reducing the number of variables requiring optimisation and was found to produce stable results.
In addition to the utilisation of Monte Carlo sampling to enhance uniqueness, uniqueness was ensured by repeating the inverse analysis twelve times. For each run, the baseline and boundaries of the solution space where altered within sensible but exploratory limits, (0 < p < 4s), where p is the optimized parameter and s is the best fit parameter value.
To illustrate the results of the test method, one human eye and one porcine eye were tested and the results analysed to produce estimates of material parameters for different regions of the eye globes. The human eye of a 69-year-old male with no known ocular diseases was obtained fresh from the Fondazione Banca degli Occhi del Veneto, Italy. The porcine eye was obtained within 6 h of slaughter. Both eye globes were stored in 6 % Dextran prior to the test. The aqueous and vitreous were removed through the needle and the intra-ocular void was filled with 6 % Dextran. Once contained within the gelatine and connected to the test rig, the IOP was altered by the control system described above. Both the human and porcine eyes were loaded and unloaded from 2 to 60 mmHg at a rate of 40 mmHg/min, with a rest period of 1 min following each of 10 loading cycles. The peak load of 60 mmHg, which was above the normal physiological range, was chosen to reach the higher IOP ranges associated with eye rubbing and tonometry.
Numerical parameters derived to represent the regional variation in material response for the human and porcine eye globes based on the Ogden model
Posterior pole region
Biomechanical representation can be considered by the resulting contour plots of stress and strain across the eye globe, as provided in Fig. 8 at 60 mmHg. In the model, the lamina cribrosa exhibits the greatest stress and strain due to its mechanical compliance and linear representation. Across the stroma, the cornea exhibits the greatest strain in the model. Unlike the derived material behavioural properties, there is no significant trend in the strain from anterior to posterior sclera. The stress in the cornea at IOP = 60 mmHg is similar to the majority of the sclera. However, there are areas of the sclera, which exhibit higher stress particularly in the temporal region. The region of highest stromal stress is in the limbus.
The inverse analysis for the porcine eye provided material parameter estimations with a RMS error at 26 % due to the significantly larger latitudinal variations in displacement within each eye region. The values of material parameters, α and μ, obtained for the tested eye are presented in Table 2.
A significant benefit of intact globe testing is the ability to estimate the regional variation of stiffness across the cornea, limbus and sclera. Previous experimental analysis of these variations has mostly been limited to separated corneas and scleras, obtained from different donors in most cases. The limbal region is commonly used for clamping the separated cornea and sclera specimens and hence is not usually characterised. Additionally, the clamps provide unrealistic boundary conditions, which are likely to affect the behaviour obtained experimentally in the adjacent areas. Therefore, the procedure described in this study is of particular benefit for obtaining material stiffness properties at the corneoscleral junction, providing the eye with physiologic loading and supporting conditions and correlating the behaviour in the cornea to that in the sclera.
The main objective of this study was to provide a method for determining the material stiffness of the eye globe, which more accurately represents the in vivo state. The microstructural effects of tissue manipulation (as is inherent in strip uniaxial testing), and the associated deviation from in vivo strain distributions, have led in other test methods to the necessity for preconditioning before a consistent, cyclic stress–strain behaviour is observed . The realignment of collagen fibrils during load cycles has been characterised by Quinn and Winkelstein  for ligaments, which also consist mainly of collagen fibrils. Quinn and Winkelstein’s study demonstrated a strong correlation between changes in collagen fibril alignment and changes in the mechanical response during preconditioning cycles as fibrils are able to change direction towards that of the principal strain. The pressure-displacement measurements of the eye globes tested suggest that no stiffening occurred between the 10 repetitive cycles, possibly as a result of no microstructural changes taking place. Subsequently, there is reduced need for pre-conditioning with the intact eye inflation procedure, which loads and supports the eye in a similar fashion to the in vivo conditions.
The limited experimental testing of one human and one porcine eye, included in this study, was presented mainly to illustrate the whole-globe inflation and the results that can be obtained from it. While the results presented show interesting trends that could only be obtained through whole-globe testing, they will need confirmation in studies involving larger numbers of test specimens.
In the test conducted on a human eye, inverse analysis allowed observation of high stiffness at the limbal region and the ring of scleral tissue immediately adjacent to it, reducing gradually towards the posterior pole. At the same time, the central cornea had a higher stiffness than the peripheral cornea, and both were much more compliant than the limbus. This behaviour is compatible with trends in material properties reported in previous studies [21, 27, 33]. While there was a decrease in the non-linear stiffness parameter, α, from the central to peripheral cornea, analysis revealed an increase in the initial shear modulus parameter, μ. This suggests a consistency with findings of previous studies which describe increased interweaving of collagen lamellae in the peripheral cornea . On the other hand, while the derived material stiffness parameters reduce from anterior to posterior sclera, observation of the stress and strain maps reveal relative consistency across these regions – suggesting that these changes are compensated by the variation in thickness (geometric contribution to stiffness).
The inverse analysis included within the method description and performed in this study assumed isotropy of material behaviour. This is despite the fact that the mechanical response of ocular tissue is known to be anisotropic and that a number of studies have utilised characteristic anisotropic microstructure arrangements [51, 52] in numerical representation of the cornea [53, 54]. However, since there was no data providing characteristic microstructure relationships covering the entire eye globe, it was not possible in this study to adopt an anisotropic numerical representation of the tissue.
The results obtained immediately adjacent to the posterior pole may have been compromised by the insertion of the inflation needle. Due to the process of reorienting the DIC measurement such that the origin was no longer located at the posterior pole, the potential inaccuracies in this area did not compound smaller errors in the derivation of material parameters elsewhere across the eye globe. The response to increases in IOP is a global response and errors in the representation of one region will always affect another region. However, it is suggested that this transfer of the origin to the limbal region has improved the reliability in the derivation of properties particularly at the limbus and across the cornea as this effect has been limited.
For the porcine eye, while the experimental procedure was applied successfully, the inverse analysis produced larger errors due to considerable latitudinal variation in behaviour within each of the model’s nine regions. This is most prominent in the mid cornea, limbal and equatorial regions. While the non-rotational representation of the specimen specific geometry was considered, the FEM was limited to rotationally symmetric material representation, which meant that the variation in stiffness could only arise due to the slightly irregular geometry around latitudinal regions. However, it is clear from the fitting of numerical to experimental data in the porcine eye that variations in geometry alone did not entirely represent the variation in stiffness, suggesting the need to consider possibly significant variations in material stiffness around these latitudinal regions when modelling porcine eye globes. This makes it more important to consider the anisotropic fibril distribution across the ocular surface once this information becomes available.
The study demonstrated that preconditioning cycles may not be required to obtain consistent load-deformation behaviour in globe inflation testing. This is a result of the test procedure providing a strain distribution that closely represents that of the in vivo eye and does not involve significant surgical manipulation as those factors would otherwise result in fibre re-orientation, and hence tissue stiffening, between successive preconditioning cycles. These findings are consistent with an earlier study on sections of non-human ocular tissue .
While the test procedure presented in this study represents a significant improvement in the ability to derive regional stiffness variations consistent with the in vivo eye, the method including the experimental testing and associated data analysis is substantially more demanding and more complex than uniaxial strip extensometry and inflation of separated corneas and scleras. However, despite these challenges and the difficulties in acquiring whole eye globes for research, the novelty of the technique and the more comprehensive data it offers, make the method invaluable in ocular material characterisation, especially as microstructure data covering the whole ocular surface becomes available.
The research was partially supported by the Engineering and Physical Sciences Research Council of the UK.
CW was the primary researcher in the development of the study, the analysis and interpretation of the data, and the originator of the manuscript. AJ participated in the development of the methodology and the data analysis. SJ contributed to the design of the test equipment and assisted drafting the manuscript. FB contributed to data analysis and interpretation. PR provided the human eye globe, advised on the methodology development. As group lead, AE supervised all aspects of the study and contributed significantly to the development of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis 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.
- Elsheikh A, Wang D, Kotecha A, Brown M, Garway-Heath D. Evaluation of goldmann applanation tonometry using a nonlinear finite element ocular model. Ann Biomed Eng. 2006;34:1628–40.View ArticlePubMedGoogle Scholar
- Kapnisis K, Van Doormaal M, Ross EC. Modeling aqueous humor collection from the human eye. J Biomech. 2009;42:2454–7.View ArticlePubMedGoogle Scholar
- Roy AS, Dupps WJ. Effects of altered corneal stiffness on native and postoperative LASIK corneal biomechanical behavior: a whole-eye finite element analysis. J Refract Surg. 2009;25:875–87.View ArticlePubMedGoogle Scholar
- Elsheikh A, Alhasso D, Gunvant P, Garway-Heath D. Multiparameter correction equation for goldmann applanation tonometry. Optom Vis Sci. 2011;88:102–12.View ArticleGoogle Scholar
- Norman RE, Flanagan JG, Sigal IA, Rausch SMK, 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:4–12.View ArticlePubMedGoogle Scholar
- Roy AS, Dupps WJ. Patient-specific modeling of corneal refractive surgery outcomes and inverse estimation of elastic property changes. J Biomech Eng. 2011;133:011002.View ArticlePubMedGoogle Scholar
- Abahussin M, Hayes S, Knox Cartwright NE, Kamma-Lorger CS, Khan Y, Marshall J, et al. 3D collagen orientation study of the human cornea using X-ray diffraction and femtosecond laser technology. Invest Ophthalmol Vis Sci. 2009;50:5159–64.View ArticlePubMedGoogle Scholar
- Komai Y, Ushiki T. The three-dimensional organization of collagen fibrils in the human cornea and sclera. Invest Ophthalmol Vis Sci. 1991;32:2244–58.PubMedGoogle Scholar
- Boyce BL, Jones RE, Nguyen TD, Grazier JM. Stress-controlled viscoelastic tensile response of bovine cornea. J Biomech. 2007;40:2367–76.View ArticlePubMedGoogle Scholar
- Downs JC, Suh JK, Thomas KA, Bellezza AJ, Hart RT, Burgoyne CF. Viscoelastic material properties of the peripapillary sclera in normal and early-glaucoma monkey eyes. Invest Ophthalmol Vis Sci. 2005;46:540–6.View ArticlePubMedGoogle Scholar
- Maurice DM. The structure and transparency of the cornea. J Physiol. 1957;136:263–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Cogan DG. Applied anatomy and physiology of the cornea. Trans Am Acad Ophthalmol Otolaryngol. 1951;55:329–59.PubMedGoogle Scholar
- Borcherding MS, Blacik LJ, Sittig RA, Bizzell JW, Breen M, Weinstein HG. Proteoglycans and collagen fibre organization in human corneoscleral tissue. Exp Eye Res. 1975;21:59–70.View ArticlePubMedGoogle Scholar
- Watson PG, Young RD. Scleral structure, organisation and disease. A review. Exp Eye Res. 2004;78:609–23.
- Ethier CR, Johnson M, Ruberti J. Ocular biomechanics and biotransport. Annu Rev Biomed Eng. 2004;6:249–73.View ArticlePubMedGoogle Scholar
- Girard MJ, Suh JK, Bottlang M, Burgoyne CF, Downs JC. Scleral biomechanics in the aging monkey eye. Invest Ophthalmol Vis Sci. 2009;50:5226–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Fratzl P, Misof K, Zizak I, Rapp G, Amenitsch H, Bernstorff S. Fibrillar structure and mechanical properties of collagen. J Struct Biol. 1998;122:119–22.View ArticlePubMedGoogle Scholar
- Govrin-Yehudain J, Levanon D, Hirshowitz B. Collagen fiber alignment in the human preauricular dermis before and after rhytidectomy--a scanning electron microscopy study. Ann Plast Surg. 1986;16:200–3.View ArticlePubMedGoogle Scholar
- Girton TS, Barocas VH, Tranquillo RT. Confined compression of a tissue-equivalent: collagen fibril and cell alignment in response to anisotropic strain. J Biomech Eng. 2002;124:568–75.View ArticlePubMedGoogle Scholar
- Downs JC, Suh JK, Thomas KA, Bellezza AJ, Burgoyne CF, Hart RT. Viscoelastic characterization of peripapillary sclera: material properties by quadrant in rabbit and monkey eyes. J Biomech Eng. 2003;125:124–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Elsheikh A, Geraghty B, Alhasso D, Knappett J, Campanelli M, Rama P. Regional variation in the biomechanical properties of the human sclera. Exp Eye Res. 2010;90:624–33.View ArticlePubMedGoogle Scholar
- Montiani-Ferreira F, Li T, Kiupel M, Howland H, Hocking P, Curtis R, et al. Clinical features of the retinopathy, globe enlarged (rge) chick phenotype. Vision Res. 2003;43:2009–18.View ArticlePubMedGoogle Scholar
- Elsheikh A, Brown M, Alhasso D, Rama P, Campanelli M, Garway-Heath D. Experimental assessment of corneal anisotropy. J Refract Surg. 2008;24:178–87.PubMedGoogle Scholar
- Hoeltzel DA, Altman P, Buzard K, Choe K. Strip extensiometry for comparison of the mechanical response of bovine, rabbit, and human corneas. J Biomech Eng. 1992;114:202–15.View ArticlePubMedGoogle Scholar
- Elsheikh A, Kassem W, Jones SW. Strain-rate sensitivity of porcine and ovine corneas. Acta Bioeng Biomech. 2011;13:25–36.PubMedGoogle Scholar
- Curtin BJ. Physiopathologic aspects of scleral stress–strain. Trans Am Ophthalmol Soc. 1969;67:417–61.PubMedPubMed CentralGoogle Scholar
- Elsheikh A, Alhasso D, Rama P. Biomechanical properties of human and porcine corneas. Exp Eye Res. 2008;86:783–90.View ArticlePubMedGoogle Scholar
- Boyce BL, Grazier JM, Jones RE, Nguyen TD. Full-field deformation of bovine cornea under constrained inflation conditions. Biomaterials. 2008;29:3896–904.View ArticlePubMedGoogle Scholar
- Bryant MR, McDonnell PJ. Constitutive laws for biomechanical modeling of refractive surgery. J Biomech Eng. 1996;118:473–81.View ArticlePubMedGoogle Scholar
- Elsheikh A, Wang D, Brown M, Rama P, Campanelli M, Pye D. Assessment of corneal biomechanical properties and their variation with age. Curr Eye Res. 2007;32:11–9.View ArticlePubMedGoogle 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:051011.View ArticlePubMedPubMed CentralGoogle Scholar
- Myers KM, Coudrillier B, Boyce BL, Nguyen TD. The inflation response of the posterior bovine sclera. Acta Biomater. 2010;6:4327–35.View ArticlePubMedGoogle Scholar
- Coudrillier B, Tian J, Alexander S, Myers KM, Quigley HA, Nguyen TD. Biomechanics of the human posterior sclera: age- and glaucoma-related changes measured using inflation testing. Invest Ophthalmol Vis Sci. 2012;53:1714–28.View ArticlePubMedPubMed CentralGoogle Scholar
- Elsheikh A, Anderson K. Comparative study of corneal strip extensometry and inflation tests. J R Soc Interface. 2005;2:177–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Lari DR, Schultz DS, Wang AS, Lee OT, Stewart JM. Scleral mechanics: comparing whole globe inflation and uniaxial testing. Exp Eye Res. 2012;94:128–35.View ArticlePubMedGoogle Scholar
- Mattson MS, Huynh J, Wiseman M, Coassin M, Kornfield JA, Schwartz DM. An in vitro intact globe expansion method for evaluation of cross-linking treatments. Invest Ophthalmol Vis Sci. 2010;51:3120–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Kling S, Remon L, Pérez-Escudero A, Merayo-Lloves J, Marcos S. Corneal biomechanical changes after collagen cross-linking from porcine eye inflation experiments. Invest Ophthalmol Vis Sci. 2010;51:3961–8.View ArticlePubMedGoogle Scholar
- Jaycock PD, Lobo L, Ibrahim J, Tyrer J, Marshall J. Interferometric technique to measure biomechanical changes in the cornea induced by refractive surgery. J Cataract Refract Surg. 2005;31:175–84.View ArticlePubMedGoogle Scholar
- Bayless KJ, Kwak HI, Su SC. Investigating endothelial invasion and sprouting behavior in three-dimensional collagen matrices. Nat Protoc. 2009;4:1888–98.View ArticlePubMedGoogle Scholar
- Bennett G, Rabbetts RB. Clinical Visual Optics. 4th ed. Edinburgh: Butterworth-Heinemann; 2007.Google Scholar
- Hjortdal JO, Jensen PK. In vitro measurement of corneal strain, thickness, and curvature using digital image processing. Acta Ophthalmol Scand. 1995;73:5–11.View ArticlePubMedGoogle Scholar
- Sigal IA, Flanagan JG, Ethier CR. Factors influencing optic nerve head biomechanics. Invest Ophthalmol Vis Sci. 2005;46:4189–99.View ArticlePubMedGoogle Scholar
- Yu JG, Bao FJ, Feng YF, Whitford C, Ye T, Huang YB, et al. Assessment of corneal biomechanical behavior under posterior and anterior pressure. J Refract Surg. 2013;29:64–70.View ArticlePubMedGoogle Scholar
- Groves RB, Coulman SA, Birchall JC, Evans SL. Quantifying the mechanical properties of human skin to optimise future microneedle device design. Comput Methods Biomech Biomed Engin. 2012;15:73–82.View ArticlePubMedGoogle Scholar
- Moran R, Smith JH, García JJ. Fitted hyperelastic parameters for Human brain tissue from reported tension, compression, and shear tests. J Biomech. 2014;47:3762–6.View ArticlePubMedGoogle Scholar
- Kampmeier J, Radt B, Birngruber R, Brinkmann R. Thermal and biomechanical parameters of porcine cornea. Cornea. 2000;19:355–63.View ArticlePubMedGoogle Scholar
- Zhou B, Xu F, Chen CQ, Lu TJ. Strain rate sensitivity of skin tissue under thermomechanical loading. Philos Trans A Math Phys Eng Sci. 2010;368:679–90.View ArticlePubMedGoogle Scholar
- Cheng S, Clarke EC, Bilston LE. The effects of preconditioning strain on measured tissue properties. J Biomech. 2009;42:1360–2.View ArticlePubMedGoogle Scholar
- Quinn KP, Winkelstein BA. Preconditioning is correlated with altered collagen fiber alignment in ligament. J Biomech Eng. 2011;133(6):064506.View ArticlePubMedGoogle Scholar
- Abass A, Hayes S, White N, Sorensen T, Meek KM. Transverse depth-dependent changes in corneal collagen lamellar orientation and distribution. J R Soc Interface. 2015;12(104):20140717.View ArticlePubMedPubMed CentralGoogle Scholar
- Aghamohammadzadeh H, Newton RH, Meek KM. X-ray scattering used to map the preferred collagen orientation in the human cornea and limbus. Structure. 2004;12:249–56.View ArticlePubMedGoogle Scholar
- Boote C, Hayes S, Abahussin M, Meek KM. Mapping collagen organization in the human cornea: left and right eyes are structurally distinct. Invest Ophthalmol Vis Sci. 2006;47:901–8.View ArticlePubMedGoogle Scholar
- Studer H, Larrea X, Riedwyl H, Büchler P. Biomechanical model of human cornea based on stromal microstructure. J Biomech. 2010;43:836–42.View ArticlePubMedGoogle Scholar
- Whitford C, Studer H, Boote C, Meek KM, Elsheikh A. Biomechanical model of the human cornea: considering shear stiffness and regional variation of collagen anisotropy and density. J Mech Behav Biomed Mater. 2015;42:76–87.View ArticlePubMedGoogle Scholar
- Tonge TK, Murienne BJ, Coudrillier B, Alexander S, Rothkopf W, Nguyen TD. Minimal Preconditioning Effects Observed for Inflation Tests of Planar Tissues. J Biomech Eng. 2013;135(11):114502.View ArticlePubMedGoogle Scholar