## 1.

## Introduction

It is much easier to construct a numerical model of the human eyeball nowadays than it was
$20\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}30\phantom{\rule{0.3em}{0ex}}\text{years}$
ago due to advances in professional software programs for construction design based on finite element analysis (FEA). Such programs allow the construction of models of a complete eye globe with complicated geometry that possess both nonlinear and anisotropic materials.^{1, 2} However, most of the work done with such numerical tools has been dedicated to the specific elements of the eye globe, and in particular to the cornea and its mechanical properties, before and after corneal refractive surgery.^{3, 4, 5} Also, there has been a continuous effort to model and analytically describe the mechanical properties of the eye using closed-form expressions.^{6, 7}

Thus, it would appear that biomechanical modeling of the eyeball should have been resolved and that the accuracy of such a model would only be a matter of performing a larger number of iterations to find the numerical solution. However, notwithstanding the efforts that are being made toward development of accurate algorithms, current eye models are still far from describing the behavior of a real human eye. The reason for this is rather prosaic. Despite the eye’s apparent structural simplicity, the *in vivo* mechanical properties of the human eye have not been categorically described. Using the words of Shin
^{8} in referring to the cornea, “experimental data is [often] meagre and flawed.” Detailed knowledge of the structure of all eyeball tissue components (i.e., cornea, sclera, iris, lens, choroid, and retina) as well as their material properties (nonlinear and anisotropic) is needed. Two of these eyeball components, the cornea and the sclera, are fairly well described when it comes to their geometries, but their material properties have been disputed.
^{8, 9, 10, 11, 12, 13, 14} On the other hand, knowledge of the *in vivo* mechanical properties of the other components of the eye is very limited. As a result, most reported biomechanical models of the human eyeball contain the cornea and the sclera, and sometimes the optic nerve,^{15} but omit the other components.

Another deficiency in the majority of reported biomechanical models of the eyeball is the developers’ apparent low interest in the optics of the eye. Most authors develop their models as if the eye was a strictly mechanical structure while its optics was just a simple consequence of the mechanical parameters, which are often chosen *ad hoc*. Subsequently, this often leads to problems with accurate mechanical identification. To date, optical properties of the eye are hardly used to identify the mechanical parameters of the eyeball model.^{16}

In this work we postulate that the optical functions of the eye have strongly influenced its mechanical structure; that it is the optics of the eyeball that took, in the evolution process, control over the mechanics (in terms of geometry and material parameters) and subjected it toward its needs.^{17} It is well known from studies on a variety of species that there is a direct role of vision in the material properties of the eye and in regulation of its growth.^{18, 19} Therefore, it is highly likely that the material properties of the eye must have been selected by nature to optimize the optics and visual function in the first place. Thus, biomechanical modeling of the eyeball should include some relationships between the mechanical and optical parameters of the eye so the goal of developing an accurate mechanical eyeball model can be achieved.

Using this concept, we develop a numerical model of an eyeball in which the structural parameters (in terms of the eyeball geometry and its material properties) are linked to the optical function of the eye. This is achieved by controlling in the model all essential optical functions of the eye while still assuming material properties from the range of those available in the literature experimental results.^{20, 21, 22, 23} By combining the mechanical and optical functions of the eye in the model, the number of free model parameters significantly decreases, which makes the model more robust and the comparison between any two model conditions more tractable.

## 2.

## Methods

A model of the eyeball is defined as a set of geometrical and material parameters in the design of a mechanical construction that can be solved using the finite element method (FEM). In this view, the type of the elements used or the particular solution parameters of the FEM is of some importance, as they may have a significant impact on the accuracy of the solution.

## 2.1.

### The FEM Model

The model that we use to link the mechanical properties of the eye to its optical functions is built from a 2-D quadrilateral 8-node body of revolution elements. More than 500 of such elements were used. The solution parameters include both the nonlinear material effects as well as the changes in the construction configuration under loading.

A typical distribution of the density of elements used in our model is shown in Fig. 1
. The boundary conditions (model mounting) limit the movements of the construction as a rigid body; otherwise the model is allowed to move freely. The model includes three kinds of material properties corresponding to the cornea, sclera, and limbus (the circled area in Fig. 1). The only material that can remain linear in a model and satisfy optical self-adjustment is that of the limbus. Its Young’s modulus is set to
$E=1\phantom{\rule{0.3em}{0ex}}\mathrm{MPa}$
. All numerical calculations have been performed using Cosmos/M™.^{24}

## 2.2.

### Assumption of Optical Self-Adjustment

The phenomenon of optical self-adjustment of the eye can be observed during an increase in intraocular pressure (IOP). One would expect that an increase in IOP could result in a slight expansion of the eyeball and the following sequence of actions:

1. Corneal apex moving away from the retina,

2. Slight change in corneal curvature,

3. Slight axial movement of the limbus (that is internally attached to crystalline lens), and

4. Increase in the axial distance between the cornea and the crystalline lens (see parameter $d$ in Fig. 1).

^{17}and we used it in our earlier work where a solution for a linear model of the eyeball was devised.

^{16}Recent clinical studies of McMonnies and Boneham provide further support for the postulate, though it has not been proved.

^{25, 26}Under the assumptions given in the postulate, our preliminary study indicated that the following structural parameters of the eyeball have the largest impact on the property of self adjustment:

## 2.3.

### Model Geometry

We approximate the anterior and posterior surfaces of the cornea by spheres similar to the Gullstrand—Le Grand model,^{27} making sure that the data follow the available knowledge based on experimental measurements. The nominal parameter values for the considered eyeball are shown in Table 1
for a model under no load. Changes of these parameters with the increasing level of IOP will be shown later. The sclera is modeled as a
$12.5\text{-}\mathrm{mm}$
anterior radius sphere with varying thicknesses:
$0.8\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$
near the limbus,
$0.6\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$
in the equatorial area, and about
$1\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$
in the foveal region (see Fig. 1).

## Table 1

Main parameters of the considered model of the eyeball.

Parameter | Source (e.g.) | Value |
---|---|---|

Central anterior radius of curvature | Kiely
^{28} | $R=7.84\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ |

Central posterior radius of curvature | Edmund^{20} | $r=6.75\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ |

Central corneal thickness | Liu
^{22} | $\mathrm{CCT}=0.52\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ |

Peripheral corneal thickness adjacent to the limbus | Dubbelman
^{29} | $\mathrm{PCT}=0.65\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ |

Diameter of the cornea | Baumeister
^{30} | ${D}_{c}=11.5\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ |

Average corneal refractive index | Patel
^{21} | ${n}_{c}=1.3771$ |

Average refractive index of aqueous humour and vitreous body | Le Grand and El Hage^{27} | ${n}_{a}=1.336$ |

Refractive power of the crystalline lens | Le Grand and El Hage^{27} | ${F}_{\mathit{\text{lens}}}=22.07\phantom{\rule{0.3em}{0ex}}\mathrm{D}$ |

Poisson ratio | (see text) | $\nu =0.49$ |

Nominal intraocular pressure | Eysteinsson
^{23} | $\mathrm{IOP}=16\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}$ |

The corneal surfaces are better modeled by aspheric surfaces, such as an ellipsoid,^{31, 32} because those shapes correspond better to keratometric measurements of real corneas than a sphere. In our previous study on optical self-adjustment,^{16} we did consider an elliptic shape with eccentricity
$e=0.5$
. This result, however, was achieved for a linear model. Here, we decided to use the classical spherical model of the cornea so we could first study the behavior of a model with nonlinear materials.

## 2.4.

### Boundary Conditions

There is no agreement in the literature on how the model of the eyeball should be mounted despite the fact that there are many studies devoted to orbital mechanics of the eye.^{33, 34, 35, 36} For obvious reasons, those analytical models that consider the corneal surface only have either a roller or a fixed support.^{12, 13, 14} However, such simplified mountings are also considered in numerical models.^{3, 8, 37} As we will show later, this could lead to misinterpretation of the numerical results, particularly when dealing with the optical function of the eye. A recent *in vivo* study on corneal rigidity^{38} supports the assumption of free mounting the eyeball, which we will also adopt here.

## 2.5.

### Materials

The anterior surfaces of cornea and sclera have been well documented and their physical structures are known.^{39} However, the knowledge of their mechanical properties is less complete. It is known that the materials are nonlinear, anisotropic, and that under a large load, their rheological properties must be taken into account. Fortunately, not all of these aspects need always to be considered. In particular, for a model of a normal living eye, one can assume an isotropic and elastic material. However, this material must be *nonlinear*. Here we also restrict ourselves to only one layer of corneal tissue, and from this point onward we identify the cornea with the stroma. This simplification does not significantly reduce the accuracy of the solution,^{2} since it appears that the Young’s modulus of the Descemet’s membrane is much smaller than that of the stroma.^{40}

The stress-strain functional relationship,
$\sigma =f\left(\epsilon \right)$
, that can be used to characterize the eye surface materials can be expressed in the known exponential form^{39, 41}

## Eq. 1

$$\sigma =A[\mathrm{exp}\left(\alpha \epsilon \right)-1],\phantom{\rule{1em}{0ex}}\epsilon \u2a7e0,$$^{39}The constitutive Eqs. 1, 3 approximate mechanical functions of a real human cornea, particularly in a clinically relevant range of IOP. The secant modulus of elasticity is often expressed in a form ${E}_{\text{Secant}}=\sigma \u2215\epsilon $ . In our case, the secant modulus of elasticity is given byDespite the fact that expression 4 is a function, we will still call it a modulus. It should, however, be distinguished from the Young’s modulus, which is constant and used for linear materials only.

In the analysis of the stress state present in ocular surfaces under IOP, we also need to consider the value of another variable: the Poisson’s ratio. The limit of the Poisson’s ratio for an isotropic incompressible material (e.g., water) is
$\nu =0.5$
, and in some works it has been used to describe corneal tissue.^{3, 42} However, this limit value cannot be used in our modeling, because the stress-strain relationship in this case is no longer uniquely described. For this reason, we consider Poisson’s ratio values of the eye constructions values that are just below this limit. The problem is that Poisson’s ratios ranging from
$\nu =0.15$
through
$\nu =0.3$
to
$\nu =0.49$
are being reported.^{2, 8, 37, 43, 44} However, this wide range of values is being disputed, because some of the considered values are supported by real measurements.^{1, 45} In some cases, the use of real measurements and the belief that they are correct may lead to an unrealistic situation. See, for example, the work of Pandolfi and Manganiello,^{37} where the chosen values of the secant modulus of elasticity and the Poisson’s ratio led to unrealistically low values of the bulk modulus. In our opinion, such a discussion is of no serious significance because the model solution only marginally depends on the Poisson’s ratio. If we assume that the bulk modulus of elasticity of considered soft tissues is of the same order as that of water (
$K=2200\phantom{\rule{0.3em}{0ex}}\mathrm{MPa}$
), that is
${10}^{3}\phantom{\rule{0.3em}{0ex}}\mathrm{MPa}$
, while their secant (Young’s) modulus of elasticity is about three orders smaller, i.e.,
${10}^{0}\phantom{\rule{0.3em}{0ex}}\mathrm{MPa}$
, then the Poisson’s ratio is

The material parameters of the considered biomechanical model of the eyeball have been chosen in a way so that the model results match those of the Goldmann applanation tonometry.^{46} Despite its deficiencies it is still the most widely accepted method of determining intraocular pressure.^{47, 48} In this technique, the cornea is modeled as a membrane filled with fluid under certain pressure. Then, using a specially calibrated probe, a small force is applied to the central
$3.06\text{-}\mathrm{mm}$
-diameter flattened corneal area, which in turn exerts a force directly related to the internal pressure as per the so-called Imbert-Fick law.

Identification of the corneal material parameters with the Goldmann’s test showed that, to adequately represent the performance of a real eye in our modeling, the following criteria must be met:

1. The material [function $\sigma =f\left(\epsilon \right)$ ] must be strongly nonlinear (even under the physiological load).

2. The modulus ${E}_{0}$ from Eq. 2 must be relatively small (of an order of magnitude smaller than the average secant modulus of the cornea, ${E}_{\text{Secant}}$ , for $\mathrm{IOP}=16\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}$ .

3. The material elasticity ${E}_{\text{Secant}}$ for $\epsilon <0$ must be much smaller than that for stretching, as described in Eq. 3.

## Eq. 5

$${M}_{s}:\phantom{\rule{0.3em}{0ex}}A=500\phantom{\rule{0.3em}{0ex}}\mathrm{Pa},\phantom{\rule{1em}{0ex}}\alpha =55,\phantom{\rule{1em}{0ex}}\text{for}\phantom{\rule{0.3em}{0ex}}\epsilon \u2a7e0,\text{and}\phantom{\rule{0.3em}{0ex}}{E}_{0}=0.0275\phantom{\rule{0.3em}{0ex}}\mathrm{MPa}\phantom{\rule{0.3em}{0ex}}\text{for}\phantom{\rule{0.3em}{0ex}}\epsilon <0;$$## Eq. 6

$${M}_{h}:\phantom{\rule{0.3em}{0ex}}A=200\phantom{\rule{0.3em}{0ex}}\mathrm{Pa},\phantom{\rule{1em}{0ex}}\alpha =130,\phantom{\rule{1em}{0ex}}\text{for}\phantom{\rule{0.3em}{0ex}}\epsilon \u2a7e0,\text{and}\phantom{\rule{0.3em}{0ex}}{E}_{0}=0.0260\phantom{\rule{0.3em}{0ex}}\mathrm{MPa}\phantom{\rule{0.3em}{0ex}}\text{for}\phantom{\rule{0.3em}{0ex}}\epsilon <0.$$The assumption that the model satisfies the principle of optical self-adjustment leads to another structural condition. Specifically, for a given stress value, the secant modulus of the sclera must be several times larger than that of the cornea. By denoting this relationship by

## Eq. 7

$$m=\frac{{E}_{\text{Secant}}\left(\text{sclera}\right)}{{E}_{\text{Secant}}\left(\text{cornea}\right)},$$## Eq. 8

$$\sigma =A[\mathrm{exp}\left(m\alpha \epsilon \right)-1],\phantom{\rule{1em}{0ex}}\epsilon \u2a7e0,$$^{49}and has its confirmation in experimental studies.

^{9}

We have numerically validated this proposed nonlinear model by comparing its stiffness to that of a real eye.^{38, 50} In this particular aspect, material
${M}_{h}$
showed much better performance than that given by material
${M}_{s}$
. That is why in the remaining part of this paper we will only consider the corneal model realized from material
${M}_{h}$
. The material for the sclera, on the other hand, satisfies relationship 8 with
$m=5$
.

Finally, the main role of the limbus in the model is to “optically couple” the cornea and sclera so that optical self-adjustment is satisfied. The limbus is defined as an annular transition zone between the cornea and sclera. Histologically, at the limbus, the epithelium gradually thickens toward the sclera where it is replaced by conjunctival tissue. In our modeling, we simplify this limbus area to a limbus ring with mechanical properties that are between those of the cornea and sclera. Since the secant modulus of the cornea is close to $0.5\phantom{\rule{0.3em}{0ex}}\mathrm{MPa}$ and that of the sclera is $m=5$ times larger, i.e. $2.5\phantom{\rule{0.3em}{0ex}}\mathrm{MPa}$ , the modulus for the limbus must be set between these two values, say $E=1\phantom{\rule{0.3em}{0ex}}\mathrm{MPa}$ , so that the structure is not perturbed.

## 2.6.

### Optical System of the Eyeball

The total optical power of the eye depends on the respective powers of the cornea and crystalline lens and their mutual positions (i.e., distance $d$ ). Figure 1 shows the optical system and the notation used. The parameters of the eyeball are given in Table 1. Changes in the pressure $p$ influence the optical power of the cornea (by changing its axial curvature $R$ ) as well as the distance $d$ between the cornea and the crystalline lens.

The optical power of the eyeball is given by.

## Eq. 9

$${F}_{\text{eye}}={F}_{\text{cornea}}+{F}_{\text{lens}}-\frac{d\cdot {F}_{\text{cornea}}\cdot {F}_{\text{lens}}}{n},$$The displacement of the focal point with respect to the back of the eye is influenced not only by the changes in optical power of the cornea-crystalline-lens interface, but also by the “optical” elongation of the eyeball, which is understood here as a shift of the secondary principal plane ${H}_{e}^{\prime}$ with respect to the back of the eye. Thus, displacement of the focal point is composed of the axial elongation of the sclera (i.e., axial translation of the crystalline lens ${a}_{l}$ ) and the change in the distance ${S}_{H}^{\prime}$ between the crystalline lens and the second principal plane ${H}_{e}^{\prime}$ (see Fig. 1),

The distance ${S}_{H}^{\prime}$ is given bywhile the axial translation of the crystalline lens ${a}_{l}$ is obtained from the FEM solution.Finally, the displacement in focus (positive sign toward the cornea) is given by

The function $B$ can be treated as a measure of optical blur on the retina. If its change is due to the intraocular pressure $p$ , then the criterion of optical self-adjustment can be analytically expressed asThe deformation of the sclera has a crucial role in the above criterion. Thus, studying optical self-adjustment in models that contain only the corneal surface that is fixed at its perimeter is an oversimplification.The property of optical self-adjustment determines the ratio of the corneal longitudinal modulus of elasticity to that of the sclera (parameter $m$ ). Together with the material properties of the cornea, they “automatically” determine the parameters of the sclera and limbus, hence reducing the number of free parameters used in the model.

## 3.

## Results

## 3.1.

### Full Model

First, we consider a full model of the eyeball with standard geometry (see Table 1), material ${M}_{h}$ , and free movement mounting under internal load $p$ . The free model parameters include the thickness of the sclera ( $0.8\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ near the limbus, $0.6\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ in the equatorial area, and about $1\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ in the foveal region), the central corneal thickness (CCT), peripheral corneal thickness adjacent to the limbus (PCT), central anterior radius of curvature, parameters $A$ and $\alpha $ of the corneal material, and the ratio $m$ .

Figure 4 shows the solutions of the FEM for several functions dependent on $p$ . The load was applied in steps and the FEM solutions were fast-converging. The solutions include $B$ , the displacement in the focus; $R$ , the central anterior corneal radius of curvature; ${V}_{\text{Limbus}}$ and ${U}_{\text{Limbus}}$ , the axial and radial displacements of the limbus; and ${V}_{\text{Apex}}$ , the axial displacement of the corneal apex. Comparing the displacement in the focus $B\left(p\right)$ with the displacements in the corneal curvature and limbus, we note that they are all of the same order, say about $0.1\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ . However, function $B\left(p\right)$ shows that for the pressure larger than $10\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}$ , the changes in displacement of focus are much smaller; thus, the model fullfills the criterion of optical self-adjustment. Note also that because the cornea has a flexible mount at the limbus, its radius of curvature $R$ increases with the pressure $p$ .

Figure 5 shows the radial profile of the anterior cornel surface for the same seven loads considered in Fig. 4, i.e., from $0\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}40\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}$ in steps of $8\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}$ , and additionally for $p=3.2\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}$ , for which the displacements are most rapid.

## 3.2.

### Fixed Cornea Model

In many applications, it is a quite common practice to reduce the full model of the eyeball to a model of the cornea that is fixed at the limbus.^{5, 13, 47} However, such an approach significantly alters the results of the calculated optical power of the model because it significantly changes the location of the focal point with respect to the cornea. It also changes other optical functions of the eye that are important in biomechanical studies of the eyeball.

The FEA solution for the fixed cornea model is similar to that of a model with a limbal ring that has infinitely high stiffness. However, two such models are not identical, because in the fixed corneal model the sclera expands. Figure 6 indicates that the optical mechanism of this model is different than that of the full model described earlier (see Fig. 4). The most important difference is in the shape of function $\Delta R$ . That is, unlike in the full model, the central anterior corneal radius of curvature decreases with an increase in $p$ . As a consequence, the focal distance in the fixed cornea model shortens and its amplitude gets larger than in the case of the full model. Also, the displacement of the focal point differs with changes in $p$ —it moves forward much more while previously it moved slightly backward. The reason for all of these differences between the two considered models is the dynamic behavior of the anterior corneal surface when it is fixed at the limbus (see Fig. 7 ). However, it is counterintuitive that these radically different boundary conditions lead also to a model that is optically self-adjusted in the range from $16\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}32\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}$ [see function $B\left(p\right)$ in Fig. 6]. It should be noted that it is not easy to achieve optical self-adjustment for the model. For the material ${M}_{h}$ , a deviation from criterion 7 and ${E}_{\text{Limbus}}=1\phantom{\rule{0.3em}{0ex}}\mathrm{MPa}$ quickly puts the model “out of tune” in which function $B\left(p\right)$ changes strongly in the whole range of considered pressures.

## 4.

## Discussion

We have considered biomechanical models of a standard eyeball (surgically unaltered). The FEA solutions of the two considered models clearly indicate that the boundary conditions (in particular the way the cornea is fixed to the sclera) strongly influence the optics of the eyeball. The investigated case of the fixed cornea model is somewhat extreme. However, even for the model with the roller support of the cornea^{3, 13} at
$40\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}45\phantom{\rule{0.3em}{0ex}}\mathrm{deg}$
, it is not possible to achieve optical self-adjustment because this type of support does not fullfill condition 7.

Thus, the most important feature of our proposed optically inspired model is its optical self-adjustment for intraocular pressure above $10\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}$ . The principal condition for such behavior of the eye model is the constraint on the secant modulus of the sclera to be several times larger than that of the cornea. The Young’s modulus of the limbal ring, ${E}_{\text{Limbus}}=1\phantom{\rule{0.3em}{0ex}}\mathrm{MPa}$ , also had to be chosen carefully. It is quite remarkable that these purely mechanical conditions were imposed simply by the optics.

These modeling results have indicated that the increase of focal length $\Delta f$ , as well as the corresponding change in eye length $\Delta l$ , are functions of $p$ and are closely related for pressures above $10\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}$ since their difference (equal to $B$ ) is almost constant. For example, for $p=16\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}$ , the focal length changed by $\Delta f=0.264\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ (from the clinically unseen in a real eye but biomechanically relevant case where $p=0$ ). If the length of the eye had not changed at the same time, one would expect that the focal point had moved by this value. However, in the optical self-adjustment model, the focus moved only by about $B=-0.121\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ , a difference of about $0.7\phantom{\rule{0.3em}{0ex}}\text{to}\phantom{\rule{0.3em}{0ex}}0.3\phantom{\rule{0.3em}{0ex}}\mathrm{D}$ (again, from the nominal point of $p=0$ ). Note that in a real eye, where the intraocular pressure may range between 10 and $20\phantom{\rule{0.3em}{0ex}}\mathrm{mm}\phantom{\rule{0.2em}{0ex}}\mathrm{Hg}$ , the condition of optical self-adjustment results in changes of focal point that are not noticeable by a human.

The solutions for the fixed cornea model (shown in Figs. 6 and 7 indicate that it is a drastically different model. The displacement of the second principal plane $\Delta l$ , despite originating mainly from the displacement of the crystalline lens $\Delta {S}_{H}^{\prime}$ [see Eq. 12], is much larger than in the full model. However, the largest difference can be seen in changes of the focal point $\Delta f$ with respect to the initial value of $f=22.426$ (for $p=0$ ). Not only has the change radically increased, but the sign also changed. Instead of moving by $-0.121\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ (backward), it moved by $0.894\phantom{\rule{0.3em}{0ex}}\mathrm{mm}$ (forward). This leads to significant optical power changes:

The displacements of the limbus, although small, have a crucial impact on the optical power ${F}_{\text{eye}}$ of the whole system being a function of the intraocular pressure $p$ . This functional relationship has an important role in the numerical studies of the eyeball. On the other hand, the displacements of the limbus caused by changes in $p$ depend on the material properties of both the sclera and the limbal ring. Therefore, it is important that these displacements are carefully considered in numerical models of the eyeball.

## Acknowledgments

The authors wish to thank Henryk Kasprzak for his encouragement, comments, and helpful suggestions. Part of this research was supported by the Australian Research Council 2007 Linkage Grant, LX0775988.

## References

*In vivo*determination of Young’s modulus for the human cornea,” Bull. Math. Biol., 49 (2), 217 –232 (1987). 0092-8240 Google Scholar