# Patent

## Title
Method for stabilizing contact lenses

## Abstract

A method for stabilizing contact lenses includes providing a lens design with a nominal set of stabilization zone parameters, applying a merit function to the lens design based on balancing moments of momentum, and creating a contact lens design with improved stabilization based on the application of the merit functions to the lens design with a nominal set of stabilization zone parameters.

## Background
BACKGROUND OF THE INVENTION
Correction of certain optical defects can be accomplished by imparting non-spherical corrective aspects to one or more surfaces of a contact lens such as cylindrical, bifocal, or multifocal characteristics. These lenses must generally be maintained at a specific orientation while on the eye to be effective. Maintenance of the on-eye orientation of a lens typically is accomplished by altering the mechanical characteristics of the lens. Prism stabilization including decentering of the lens' front surface relative to the back surface, thickening of the inferior lens periphery, forming depressions or elevations on the lens' surface, and truncating the lens edge are examples of stabilization approaches. Additionally, dynamic stabilization has been used in which the lens is stabilized by the use of thin zones, or areas in which the thickness of the lens' periphery is reduced. Typically, the thin zones are located at two regions that are symmetric about either the vertical or horizontal axis of the lens from the vantage point of its on-eye placement.
Evaluating lens design involves making judgments concerning the performance of the lens on-eye and then optimizing the design if necessary and possible. This process is typically done by clinically evaluating the test design in patients. However, this process is time consuming and expensive because it requires a significant number of patients to be tested since patient to patient variability must be accounted for.
There is a continuing need for improving the stabilization of certain contact lenses.
SUMMARY OF THE INVENTION
The invention is a contact lens designed with improved stabilization relative to a nominal stabilized design.
In another aspect of the invention, a method for stabilizing contact lenses incorporates a lens design with a nominal set of stabilization zone parameters, evaluating the on-eye performance of the lens design, calculating a merit function based on this performance, and optimizing the stabilization zone parameters by applying the merit function. This process may be performed iteratively via a virtual model (e.g., software based) that simulates the effects of eye mechanics such as blinking and adjusts the stabilization scheme accordingly.
In yet another aspect of the invention, contact lenses are stabilized according to a scheme in which the moment of momentum of torques acting on the lens on-eye are balanced.
In yet another aspect of the invention, contact lenses are stabilized by the formation of one or more zones having a different thickness than the remainder of the lens and in which these zones are located on the lenses such that the moment of momentum of torques acting on the lens when the lens is on-eye will be balanced.
In yet another aspect of the invention a contact lens has a stabilization zone with the bulk of its length lying beneath the horizontal axis of the lens.
In yet another aspect of the invention a contact lens has a stabilization zone having a differing rate of change of slope (from its peak) in one direction relative to the other.
In yet another aspect of the invention a contact lens has a different in height profile above the horizontal axis than it does below the horizontal axis.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a front or object view of a stabilized contact lens.
FIGS. 2A-C are schematic representations of an eye with an inserted lens that identify the axis of rotation and various torques acting at the lens.
FIG. 3 is a flowchart showing a stabilization optimization process according to the invention.
FIGS. 4A-C is a front view of a stabilized lens with stabilization zones and circumferential and radial thickness graphs corresponding to example 1.
FIGS. 5A-C is a front view of a stabilized lens with stabilization zones and circumferential and radial thickness graphs corresponding to example 2.
FIGS. 6A-C is a front view of a stabilized lens with stabilization zones and circumferential and radial thickness graphs corresponding to example 3.
FIGS. 7A-C is a front view of a stabilized lens with stabilization zones and circumferential and radial thickness graphs corresponding to example 4.
FIG. 8 is a graph showing rotation speed measurements.



DETAILED DESCRIPTION
Contact lenses of this invention have designs that optimize stabilization based on balancing various forces that act on the lenses. This involves the application of a design process that balances torques operating on the eye, components of the eye, and ultimately the stabilized lens that is placed on the eye. Preferably, improved stabilization is attained by starting the improvement process with a nominal design that includes stabilization elements. For example, a lens design that has two stabilization zones that are symmetric about both the horizontal and vertical axes running through the center is a convenient reference from which to optimize stabilization of the lens according the inventive methods. By “stabilization zone” is meant an area of the peripheral zone of the lens which has thickness values larger than the average thickness of the remaining areas of the peripheral zone. By “peripheral zone” is meant the area of the lens surface circumferentially surrounding the optic zone of the lens, and extending up to but not including the edge of the lens. Another stabilization design that is a useful starting point is described in US Patent Publication 20050237482 which is incorporated herein by reference but any stabilization design can be used as the nominal design that is then optimized according to this invention. The stabilization design improvement process can also incorporate testing the improvement with the eye model described below, evaluating the results of the testing, and continuing iteratively with the improvement process until a desirable level of stabilization is achieved.
FIG. 1 depicts the front, or object side, surface of a stabilized lens. Lens 10 has an optical zone 11. The lens' periphery surrounds optic zone 11. Two thick regions 12 are located in the periphery and are stabilization zones.
The model that is preferably used in the process to produce the new designs incorporates various factors and assumptions that simulate the mechanical operation and their effects on lens stability. Preferably, this model is reduced to software using standard programming and coding techniques according to well-known programming techniques. In broad overview, the model is used in the process for designing stabilized lenses by simulating the application of the forces described below in a prescribed number of blinks of the eye. The degree to which the lens rotates and decenters is determined accordingly. The design is then altered in a way that is directed to bringing rotation and/or centration to more desirable levels. It is then subjected to the model again to determine translation upon blinking after the pre-determined number of blinks. The alteration of the design is accomplished by the application of merit functions described in more detail below.
The model assumes that the eye preferably consists of at least two spherical surface parts representing the cornea and the sclera and that the origin of the x-y-z coordinate axes is in the center of the sphere representing the cornea. More complex surfaces such as aspheric surfaces may also be used. The base shape of the lens consists of spherical surface parts but the base curve radius of the lens is allowed to change from the center of the lens towards the edge. More than one base curve may be used to describe the back surface. It is assumed that a lens positioned on the eye assumes the same shape as that of the eye. The thickness distribution of the lens need not necessarily be rotationally symmetric and indeed is not symmetric according to some preferred embodiments of the inventive lenses. Thick zones at the edge of the lens may be used to control the position and orientation behavior of the lens. A uniform thin film of liquid (tear film) exists between the lens and the eye, with a typical thickness of 5 μm. This tear film is referred to as the post-lens tear film. At the lens edge the thickness of the liquid film between the lens and eye is much smaller and is referred to as the mucin tear film. A uniform thin film of liquid (also, tear film) with a typical thickness of 5.0 μm exists between the lens and the lower and upper eyelids and these are referred to as the pre-lens tear films. The boundaries of both the lower and the upper eyelids lie in planes having unit normal vectors in the x-y plane. Hence, the projections of these boundaries on the plane perpendicular to the z-axis are straight lines. This assumption is also made during the motion of the eyelids. The upper eyelid exerts a uniform pressure on the contact lens. This uniform pressure is exerted on the whole area of the contact lens covered by the upper eyelid or on a part of this area near the boundary of the upper eyelid with uniform width (measured in the direction perpendicular to the plane through the curve describing the edge of the eyelid). The lower eyelid exerts a uniform pressure on the contact lens. This pressure is exerted on the whole area of the contact lens covered by the lower eyelid. The pressure exerted by the eyelids on the contact lens contributes to the torque acting at the lens through a non-uniform thickness distribution (thick zone) of the contact lens, especially near the edge. The effect of this pressure on the torques acting at the contact lens is referred to as the melon seed effect. Viscous friction exists in the post-lens tear film if the lens moves with respect to the eye. Viscous friction also exists in the mucin tear film between lens edge and the eye if the lens moves with respect to the eye. Additionally, viscous friction exists in the pre-lens tear film if the lens moves and/or the eyelids move. Strains and stresses in the lens occur due to the deformation of the lens. These strains and stresses result in an elastic energy content of the lens. As the lens moves with respect to the eye and the deformation of the lens changes, the elastic energy content changes. The lens tends toward a position at which the elastic energy content is minimal.
The parameters describing the geometry of the eye (cornea and sclera), base shape of the lens and the motion of the eyelids the eyelids are shown in FIG. 2. The motion of the lens follows from the balance of moment of momentum acting at the lens. Inertial effects are neglected. Then the sum of all moments acting at the lens is zero. Hence,

{right arrow over (0)}={right arrow over (M)}l,cor+{right arrow over (M)}l,muc+{right arrow over (M)}l,low+{right arrow over (M)}l,upp+{right arrow over (M)}l,Ulow+{right arrow over (M)}l,Uupp+{right arrow over (M)}l,Vupp+{right arrow over (M)}ms,low+{right arrow over (M)}ms,upp+{right arrow over (M)}elast+{right arrow over (M)}grav 

The first 4 moments are resisting torques and are linearly dependent on the lens motion. The remaining torques are driving torques. This balance of moment of momentum results in a non-linear first order differential equation for the position β of the lens




  
    
      A
      ⁡
      
        (
        
          
            β
            _
          
          ,
          t
        
        )
      
    
    ⁢
    
      
        d
        
          β
          →
        
      
      
        d
        t
      
    
  
  =
  
    
      
        M
        →
      
      total
      driving
    
    ⁡
    
      (
      
        
          β
          _
        
        ,
        t
      
      )
    
  




This equation is solved with a fourth order Runge-Kutta integration scheme. The positions of points on the contact lens follow from a rotation around the rotation vector β(t). The rotation matrix R(t) transforming the old position of points to the current position follows from Rodrigues's formula



  
    
      
        x
        →
      
      new
    
    =
    
      
        
          x
          →
        
        old
      
      +
      
        sin
        ⁢
        
          
        
        ⁢
        
          β
          ⁡
          
            (
            
              
                n
                →
              
              ×
              
                
                  x
                  →
                
                old
              
            
            )
          
        
      
      +
      
        
          (
          
            1
            -
            
              cos
              ⁢
              
                
              
              ⁢
              β
            
          
          )
        
        ⁢
        
          (
          
            
              n
              →
            
            ×
            
              (
              
                
                  n
                  →
                
                ×
                
                  
                    x
                    →
                  
                  old
                
              
              )
            
          
          )
        
      
    
  




  
    
      
        x
        →
      
      new
    
    =
    
      
        R
        ⁡
        
          (
          t
          )
        
      
      ⁢
      
        
          x
          →
        
        old
      
    
  




  
    
      where
      ⁢
      
        
      
      ⁢
      
        n
        →
      
    
    =
    
      
        
          
            β
            →
          
          
            
            
              β
              →
            
            
          
        
        ⁢
        
          
        
        ⁢
        and
        ⁢
        
          
        
        ⁢
        β
      
      =
      
        
          
          
            β
            →
          
          
        
        .
      
    
  



In the numerical integration method a time-discretization is used. Then the motion of the lens can be seen as a number of subsequent rotations, hence at the next time step tn+1 the rotation matrix is

Rn+1=RΔtRn 

where RΔt is the rotation during the time step Δt.

The rotation matrix is decomposed into a rotation Ra and a decentration Rθ of the lens

R(t)=Rθ(t)Ra(t)

The rotation of the lens is a rotation around the centerline of the lens. The decentration is a rotation around a line in the (x, y) plane. Hence, the position of the lens is seen as a rotation α of the lens around its centerline followed by a decentration {right arrow over (θ)}.
In a preferred method of the invention, merit functions (MFs) based on these relationships are fashioned to adjust and thereby improve stabilization schemes of nominal designs. These merit functions are defined based on lens on-eye performance requirements. In a preferred embodiment the merit functions are defined but not limited to: a) lens rotation and centration performance (Equation 1), b) lens stability around a resting position (Equation 2), or c) Lens rotation and centration performance and stability around resting position (Equation 3).




  
    
      
        
          MF
          1
        
        =
        
          
            
              
                
                  W
                  R
                
                ⁡
                
                  (
                  
                    Rot
                    
                      R
                      REF
                    
                  
                  )
                
              
              2
            
            +
            
              
                
                  W
                  C
                
                ⁡
                
                  (
                  
                    Cent
                    
                      C
                      REF
                    
                  
                  )
                
              
              2
            
          
        
      
    
    
      
        (
        
          Equation
          ⁢
          
            
          
          ⁢
          1
        
        )
      
    
  




By lens rotation is meant the angular movement of the lens around its z-axis occurring during and in between blinks. Rotation may be clockwise or counter clockwise depending on the initial position of the lens on the eye or the lens behavior when modeled on the eye.
By lens centration is meant the distance between the lens geometric center and the corneal apex. Centration is recorded in an x-y coordinate system in the plane of the corneal apex.
By lens stability is meant the amount of maximum lens movement in the horizontal direction (x axis) and the vertical direction (y axis) and amount of lens rotation during the blink period. Lens stability is preferably recorded with no lens misorientation and decentration after the lens reaches its final position.
Using Equation 1 as exemplary of the purpose and application of the merit function, Rot and Cent respectively describes lens performance in rotation and centration of the lens design to be optimized. RREF and CREF are variables describing the lens performance in rotation and centration of the initial lens design. WR and WC are two weighting factors allowing the adjustment of the contribution of one factor relative to the other and can take values between 0 and 1. When applied, as exemplified below, these functions are best solved numerically. The weighting factors are applied so that components of interest are given appropriate consideration. They may be equal or one component may be of more interest than another. Thus, for example, if one is more concerned with optimizing rotation than centration they would select a WR that is greater than WC. A stabilized design is improved when its merit function is decreased relative to the design that preceded it under this construct. Further, it is optimized when the merit function is minimized in such a case. Of course, one lens design may be preferable to another for reasons other than stabilization so an improved stabilization may still be undertaken according to the invention without necessarily optimizing the stabilization aspect of the design.




  
    
      
        
          MF
          2
        
        =
        
          
            
              
                
                  W
                  X
                
                ⁡
                
                  (
                  
                    
                      X
                      Range
                    
                    
                      X
                      REF
                    
                  
                  )
                
              
              2
            
            +
            
              
                
                  W
                  Y
                
                ⁡
                
                  (
                  
                    
                      Y
                      Range
                    
                    
                      Y
                      REF
                    
                  
                  )
                
              
              2
            
            +
            
              
                
                  W
                  θ
                
                ⁡
                
                  (
                  
                    
                      θ
                      Range
                    
                    
                      θ
                      REF
                    
                  
                  )
                
              
              2
            
          
        
      
    
    
      
        (
        
          Equation
          ⁢
          
            
          
          ⁢
          2
        
        )
      
    
  




In Equation 2, XRange, YRange and θRange describe lens performance in stability in horizontal direction, vertical direction and rotation of the lens design to be optimized, XREF, YREF and θREF the lens performance in stability in horizontal direction, vertical direction and rotation of the initial lens design, and WX, WY and Wθ the weighting factors allowing the adjustment of the contribution of the factors relative to each others.




  
    
      
        
          MF
          3
        
        =
        
          
            
              
                
                  W
                  R
                
                ⁡
                
                  (
                  
                    Rot
                    
                      R
                      REF
                    
                  
                  )
                
              
              2
            
            +
            
              
                
                  W
                  C
                
                ⁡
                
                  (
                  
                    Cent
                    
                      C
                      REF
                    
                  
                  )
                
              
              2
            
            +
            
              
                
                  W
                  S
                
                ⁡
                
                  (
                  
                    Stab
                    
                      S
                      REF
                    
                  
                  )
                
              
              2
            
          
        
      
    
    
      
        (
        
          Equation
          ⁢
          
            
          
          ⁢
          3
        
        )
      
    
  




In Equation 3 Rot, Cent and Stab describe lens performance in rotation, centration and stability of the lens design to be optimized, RREF, CREF and SREF the lens performance in rotation, centration and stability of the initial lens design, and WR, WC, and WS are the weighting factors allowing the adjustment of the contribution of the factors relative to each others.
In another embodiment merit functions include wearing comfort and can also include stabilization zone volume, stabilization zone surface area, soft contact lens wearer awareness to stabilization zone or any other relevant criteria.
In further preferred embodiments the merit functions are defined from the following parameters in the same fashion as those set forth above:

    
    
        Rotation performance:
        
            Surface area below the rotation curve response
            Time to reach resting position in rotation within +1-5.0 degrees
            Initial rotation speed
        
        
        Centration performance:
        
            Surface area below the centration curve response
            Time to reach resting position in centration
            First time to reach the final resting position
            Centration speed
        
        
        Stability performance:
        
            Magnitude of movement in horizontal direction
            Magnitude of movement in vertical direction
            Magnitude of rotation
            Duration of horizontal movement
            Duration of vertical movement
            Duration of rotation.
        
        
        Wearing comfort:
        
            Volume of material in excess for building stabilization zone
            Surface area covered by stabilization zone
            Lens wearer awareness to stabilization zone
        
        
    
    


There is no limitation in the type of stabilization that can be produced by the method. Stabilization zones can be of the following type:

    
    
        Symmetrical with respect to X and Y axis
        Symmetrical with respect to X or Y-axis
        Asymmetric with respect to both X and Y axis
        Constant radial distance
        Variable radial distance
    
    


Various stabilization zone parameters may be evaluated during the optimization, including without limitation the following: zone length, peak thickness location, ramp angles on either side of the peak, circumferential tilt of the zone, and zone width. The optimization parameters can also include lens diameter, base curve, thickness, optic zone diameter, peripheral zone width, material properties, other parameters describing lens features.
In a preferred embodiment of the invention, two types of improvement approaches are disclosed. In the first, a full optimization is conducted in which the model of on-eye behavior with a given iteration of the stabilization adjustment induced by the MF requires several blink cycles until the lens reaches its resting position. In another embodiment, the design is improved during a predetermined number of blink cycles. Three blink cycles are generally the least that will be effective in providing meaningful stabilization improvements. In either case, the process is conducted iteratively with the application of the MF to a nominal design. In the case in which three blink cycles are used, an initial blink orients the lens at an angle α from the horizontal, intermediate blink where the lens is oriented at an angle β from the horizontal, and final blink where the lens is positioned at the resting location. In the most preferred embodiment the angle α is set to 45 degrees and the angle β is set to 22 degrees (but both angles are not limited to these values). In another embodiment the optimization process is a combination of both approaches where a reduced number of blink cycles is preliminary used to reach an intermediate solution then several blink cycles are used to validate that the optimization has been conducted to an acceptable degree.
FIG. 3 shows a flow chart of this improvement process. The initial stabilization zones design can be either a pre-existing design or a new design. Stabilization zone parameters from these designs are determined. These parameters are obtained from calculating the design performance when the parameters are modified around their initial values. Parameters that give the most variations in the lens performance are preferably selected for the optimization process. In step 1, stabilization zone parameters are selected for consideration. These can include, for example, magnitude of the stabilization zone (Z0), peak location along 0-180 deg. meridians (r0), peak location angularly around the 0-180 deg. meridian (θ0), slope above and below the peak location, the angular length of the stabilization zone (σθ), stabilization zones rotated around the peak location, and width of the stabilization zone (θR) among many others.
In step 2, the lens is defined mathematically in terms of the stabilization zones parameters to arrive at an initial or nominal design. There is no limitation in the type of mathematic function describing the stabilization zones. Stabilization zones can also be described using computer







TABLE 1





Performance indicators obtained from designs of 


examples 1, 2, 3, and 4 applied to the merit 


functions defined by equations (1) and (2).




















Weight
WR
WC
WX
WY
Wθ





Example 1
1.00
1.00
0.50
0.50
1.70


Example 2
1.00
1.00
0.50
0.50
1.70


Example 3
1.00
1.00
0.50
0.50
1.70


Example 4
1.00
1.00
0.50
0.50
1.70





Performance indicator
RREF 
CREF
Xrange
Yrange
θrange





Example 1
505.110
1.100
1.03
2.65
1.88


Example 2
218.91
0.416
1.02
2.67
0.52


Example 3
277.22
0.356
1.03
2.68
0.67


Example 4
349.32
0.780
1.02
2.67
0.55













Merit functions
% Improvement












Equ. (1)
Equ. (2)
Equ. (1)
Equ.(2)





Example 1
1.414
1.643
N/A
N/A


Example 2
0.575
1.062
59.32
35.35


Example 3
0.637
1.106
54.96
32.68


Example 4
0.990
1.070
29.97
34.88









generated software such as CAD applications. The mathematically described design (with defined parameters) is entered into the eye model in Step 3 and rotation, centration, and stability data is generated as shown in Table 1). This data can then be used to modify one or more of the stabilization parameters in optional step 4.

The stabilization zones are modified either by reshaping, scaling, rotating, shifting, or using any other techniques to modify the current design. In steps 5a-5d the modified stabilization parameters are again run through the eye model to generate rotation, centration, and stability data for each of now modified designs. In each case of corresponding steps 6a-6d merit functions are created and applied to each new design to generate new rotation, centration, and stability data in steps 7 and 8 as the lens is maneuvered (preferably through rotation). Again, in each iteration merit functions are calculated in step 9 and checked to see if they are decreasing in step 10. A decrease is an improvement over the previous iteration. If the merit function did not decrease then a stabilization parameter can then be modified again in optional step 11 and the resulting modified lens design is then put back into the selection and data generation steps 7 and 8. If the merit function did decrease it shows an improvement in stabilization and the lens design is determined to be the final design (step 12) or other zones are again improved in optional step 13. The invention may find its greatest utility in toric and multifocal lenses. Additionally, the designs may be useful in lenses customized to a specific individual's corneal topography, lenses incorporating high order wave-front aberration correction, or both. Preferably, the invention is used to stabilize toric lenses or toric multifocal lenses as, for example, disclosed in U.S. Pat. Nos. 5,652,638, 5,805,260 and 6,183,082 which are incorporated herein by reference in their entireties.
As yet another alternative, the lenses of the invention may incorporate correction for higher order ocular aberrations, corneal topographic data, or both. Examples of such lenses are found in U.S. Pat. Nos. 6,305,802 and 6,554,425 incorporated herein by reference in their entireties.
The lenses of the invention may be made from any suitable lens forming materials for manufacturing ophthalmic lenses including, without limitation, spectacle, contact, and intraocular lenses. Illustrative materials for formation of soft contact lenses include, without limitation silicone elastomers, silicone-containing macromers including, without limitation, those disclosed in U.S. Pat. Nos. 5,371,147, 5,314,960, and 5,057,578 incorporated in their entireties herein by reference, hydrogels, silicone-containing hydrogels, and the like and combinations thereof. More preferably, the surface is a siloxane, or contains a siloxane functionality, including, without limitation, polydimethyl siloxane macromers, methacryloxypropyl polyalkyl siloxanes, and mixtures thereof, silicone hydrogel or a hydrogel, such as etafilcon A.
Curing of the lens material may be carried out by any convenient method. For example, the material may be deposited within a mold and cured by thermal, irradiation, chemical, electromagnetic radiation curing and the like and combinations thereof. Preferably, for contact lens embodiments, molding is carried out using ultraviolet light or using the full spectrum of visible light. More specifically, the precise conditions suitable for curing the lens material will depend on the material selected and the lens to be formed. Suitable processes are disclosed in U.S. Pat. No. 5,540,410 incorporated herein in its entirety by reference.
The contact lenses of the invention may be produced by any convenient method. One such method uses an OPTOFORMTM lathe with a VARIFORMTM attachment to produce mold inserts. The mold inserts in turn are used to form molds. Subsequently, a suitable liquid resin is placed between the molds followed by compression and curing of the resin to form the lenses of the invention. One ordinarily skilled in the art will recognize that any number of known methods may be used to produce the lenses of the invention.
The invention will now be further described with respect to the following non-limiting examples.
EXAMPLE 1
A contact lens having a known design for correcting the vision of astigmatic patients is shown in FIG. 6. It was designed using conventional lens design software with the following input design parameters:

    
    
        Sphere power: −3.00D
        Cylinder Power: −0.75D
        Cylinder Axis: 180 deg
        Lens diameter: 14.50 mm
        Front optical zone diameter of 8.50 mm
        Back optical zone diameter of 11.35 mm
        Lens base curve: 8.50 mm
        Center Thickness: 0.08 mm
        Eye model parameters used are listed in Table 2A and 2B.
    
    


The stabilization zone is an extra thick zone added to the thickness profile of that lens. The initial stabilization zone is constructed using a combination of normalized Gauss functions describing the radial and angular changes in thickness. The mathematical expression describing the Sag of the stabilization zone in polar coordinates is:




  
    Z
    ⁡
    
      (
      
        R
        ,
        θ
      
      )
    
  
  =
  
    
      Z
      0
    
    ·
    
      Exp
      ⁡
      
        (
        
          
            -
            0.5
          
          ·
          
            
              (
              
                
                  r
                  -
                  
                    r
                    0
                  
                
                
                  σ
                  R
                
              
              )
            
            2
          
        
        )
      
    
    ·
    
      Exp
      ⁡
      
        (
        
          
            -
            0.5
          
          ·
          
            
              (
              
                
                  θ
                  -
                  
                    θ
                    0
                  
                
                
                  σ
                  θ
                
              
              )
            
            2
          
        
        )
      
    
  




Where Z0 is the maximum magnitude of the stabilization zone, r0 and θ0 the radial and angular location of the peak, and σR and σθ are parameters controlling the profile of the thickness change in the radial and angular directions.
The change in slope along the radial and angular direction is obtained using log-normal Gauss distributions. The equation becomes:




  
    Z
    ⁡
    
      (
      
        R
        ,
        θ
      
      )
    
  
  =
  
    
      Z
      0
    
    ·
    
      Exp
      ⁡
      
        (
        
          
            -
            0.5
          
          ·
          
            
              (
              
                
                  
                    Log
                    ⁡
                    
                      (
                      r
                      )
                    
                  
                  -
                  
                    r
                    0
                  
                
                
                  σ
                  R
                
              
              )
            
            2
          
        
        )
      
    
    ·
    
      Exp
      ⁡
      
        (
        
          
            -
            0.5
          
          ·
          
            
              (
              
                
                  
                    Log
                    ⁡
                    
                      (
                      θ
                      )
                    
                  
                  -
                  
                    θ
                    0
                  
                
                
                  σ
                  θ
                
              
              )
            
            2
          
        
        )
      
    
  




The design parameters controlling the stabilization zones are:

    
    
        Change in magnitude of the stabilization zone (Z0).
        Peak location changes along 0-180 deg. meridians (r0).
        Peak location changes angularly around the 0-180 deg. meridian (θ0).
        Slope changes above and below the peak location.
        Change of the angular length of the stabilization zone (σθ).
        Stabilization zones rotated around the peak location.
        Change in width of the stabilization zone (σR) along the 0-180 deg. meridian.
    
    


The values from which the initial stabilization zone were built were:

    
    
        Z0=0.25 mm
        r0=5.75 mm
        σR=0.50 mm
        θ0=180 degrees and 0 degrees for left and right stabilization zones, respectively
        σθ=25.0 degrees
    
    


The stabilization zone was then added to the original lens thickness profile. The final maximum lens thickness was 0.38 mm. A graphic illustration of the profile is shown in FIG. 4. Stabilization zones are symmetric about both horizontal and vertical axes with slopes that uniformly descend from peak height







TABLE 2A





Initial parameters provided to the eye model.


















Tear film
Aqueous layer viscosity
8.30E−04
[Pa.s]



Mucin layer viscosity
1.50E−03
[Pa.s]



Mucin layer thickness
3.50E−07
[m]



Pre-lens tear film thickness
5.00E−06
[m]



Post-lens tear film thickness
5.00E−06
[m]


Eye
Cornea Radius
7.95E−03
[m]


geometry
Sclera Radius
1.15E−02
[m]



Visible Cornea Radius
5.82E−03
[m]



File (eye geometry)

[mm]


Lens
Lens base curve radii
8.50E−03
[m]


properties
Lens transition radii
5.50E−03
[m]



File (lens back geometry)

[mm]



Contact angle edge
−5.00
[degree]



Contact area lens edge
2.40E−05
[m2]



Density of lens material
1000
[Kg/m3]



Young's modulus
280000
[N/m2]



Poisson's ratio
0.48
[-]



File (lens normal thickness profile)

[mm]


Eyelid
Lower eyelid lateral displacement
4.00E−03
[m]


geometry
Upper eyelid lateral displacement
3.50E−03
[m]


and blink
Time to complete downward 
0.082
[s]


properties
motion upper eyelid





Time to complete blink
0.258
[s]



Time between two blinks
3
[s]



Eyelid pressure
200
[N/m2]



Lower eyelid position at start of blink
6.35E−03
[m]



Upper eyelid position at start of blink
4.70E−03
[m]



Width pressure band at 
5.00E−04
[m]



upper eyelid edge





Angle upper eyelid at start of blink
−4.47
[degree]



Angle lower eyelid at start of blink
−2.07
[degree]



Eye velocity
2
—


Gaze
Eye gaze (choice in predefined 
0
—


direction
transient gaze direction)





Amplitude gaze
20
[degree]



Frequency gaze
2.78
[Hz]


Initial
Initial rotation angle lens
0.00
[degree]


position
Initial decentration in X-direction
0.00
[m]



Initial decentration in Y-direction
0.00
[m]


Gravity
Gravity
9.80
[m/s2]


Simulation
Number of cycles that 
5
—


parameters
has to be simulated





Number of time steps in [0, Tdownblink]
−400
—



(If <0 then specified time step is used)





Specified time step
0.005
[second]



Lens discretization in radial direction
20
—



Lens discretization in 
90
—



circumferential direction


















TABLE 2B





Initial parameters provided to the eye model.







Eye velocity










1
Velocity eyelid constant



2
Velocity eyelid matching human eye velocity







Eye gaze (choice in predefined transient gaze direction)










0
No change in gaze direction



1
Horizontal motion



2
Vertical motion



3
Circular motion (counter clockwise)



4
Circular motion (clockwise)



5
Constant gaze in horizontal direction



6
Constant gaze in vertical direction










The contact lens rotation and centration characteristics were determined using the eye model described above with initial parameters provided in Table 2. Rotation of the lens steadily decreased from about 45 deg to less than 10 degrees as the number of blinks modeled went from 0 to 20. Over the course of blinks 1-20, centration remained relatively steady from about 0.06 mm to just over 0.08 mm. The resulting value of the merit function defined by equation 1 applied on the prior art lens was 1.414 with WR=WC=1.0. This example shows the rotation, centration and stability achieved by the lens of these parameters where maintenance of the on-eye orientation is achieved using depressions or elevations on the periphery of the front surface.
EXAMPLE 2
A new stabilization zone was designed using the eye model and optimization method described above and the initial design described in Example. The merit function was defined using

    
    
        Surface area below the response in rotation.
        Surface area below the response in centration.
        Identical weight for rotation and centration, WR=WC=1.0.
    
    


The values from which the initial stabilization zone were built were:

    
    
        Z0=0.25 mm
        r0=5.75 mm
        σR=0.50 mm
        θ0=180 degrees and 0 degrees for left and right stabilization zones, respectively
        σθ=25.0 degrees
    
    


The stabilization zone was then added to the original lens thickness profile.
The stabilization zone was rotated around the peak location until the lens performance characteristics represented a significant improvement over the initial design. The rotation was obtained by applying a coordinate transformation (rotation around the peak location) on the original stabilization zone coordinates:




  
    (
    
      x
      ,
      y
    
    )
  
  =
  
    
      [
      
        
          
            
              Cos
              ⁡
              
                (
                α
                )
              
            
          
          
            
              Sin
              ⁡
              
                (
                α
                )
              
            
          
        
        
          
            
              Sin
              ⁡
              
                (
                α
                )
              
            
          
          
            
              Cos
              ⁡
              
                (
                α
                )
              
            
          
        
      
      ]
    
    ⁢
    
      (
      
        
          x
          0
        
        ,
        
          y
          0
        
      
      )
    
  




Where (x0,y0) were the original coordinates and (x,y) were the new coordinates, and α the angle of rotation.
An improved stabilization design was obtained in which the final orientation of the stabilization zone was 10.0 degrees off from the vertical with the upper portion of the stabilization oriented toward the center of the lens as shown in FIG. 5. Additionally, stabilization zones are not symmetric about the horizontal axis. In this case, the bulk of the long dimension of each zone lies above the horizontal axis. The final value of the merit function was 0.58. Improvement on the merit function was of the order of 59%. Rotation declined sharply relative to the initial stabilization design. Beginning at blink 4 less than 30 degrees rotation were seen with no rotation from blink 12 onward compared to about 40-25 degrees of rotation seen in the initial design over the same range of blinks. Centration remained steady with less than 0.04 mm on blink 1 and less than 0.03 thereafter in the improved design relative to 0.06 to greater than 0.08 for the initial design over the same number of blink cycles. This example shows improved rotation, centration and stability as compared with the lens of Example 1.
EXAMPLE 3
A new stabilization zone was designed using the eye model and optimization method described above and the initial design described in Example 1. The merit function was defined using

    
    
        Surface area below the response in rotation.
        Surface area below the response in centration.
        Identical weight for rotation and centration, WR=WC=1.0.
        The values from which the initial stabilization zone were built were:
        Z0=0.25 mm
        r0=5.75 mm
        σR=0.50 mm
        θ0=180 degrees and 0 degrees for left and right stabilization zones, respectively
        σθ=25.0 degrees
    
    


The stabilization zone was added to the original lens thickness profile.
An improved stabilization design was obtained in which the final orientation of the stabilization zone was such that the peak location of the stabilization zone was changed angularly around the 0-180 deg. meridian from the geometrical center of the lens as shown in FIG. 6. Stabilization zones are no longer symmetric about the horizontal axis and the rate of change of the slope of those zones differs in direction away from the 0-180 meridian. The final value of the merit function was 0.64. Improvement on the merit function was of the order of 55%. Rotation declined sharply relative to the initial stabilization design. Beginning at blink 4 less than 30 degrees rotation were seen with about 10 degrees of rotation seen at blink 10 and no rotation from blink 16 onward compared to about 40-30-15 degrees of rotation seen in the initial design over the same range of blinks. Centration was less than 0.06 mm on blink 1 and less than 0.04 at blink 4. Thereafter it dropped sharply, being less than 0.02 at 8 blinks and zero at blink 16 relative to greater than 0.06 to greater than 0.07 and greater than 0.08 for the initial design over the same number of blink cycles. This example shows improved rotation, centration and stability as compared with the lens of Example 1.
EXAMPLE 4
A new stabilization zone was designed using the eye model and optimization method described above and the initial design described in Example 1. The merit function was defined using

    
    
        Surface area below the response in rotation.
        Surface area below the response in centration.
        Weight for rotation WR=0.84, weight for centration WC=1.14.
        The values from which the initial stabilization zone were built were:
        Z0=0.25 mm
        r0=5.75 mm
        σR=0.50 mm
        θ0=1.954
        σθ=0.14
    
    


The stabilization zone was added to the original lens thickness profile. The stabilization zone was adjusted to change the slope around the peak location. Peak locations remain on the 0-180 degrees meridian as shown in FIG. 7. Stabilization zones are not symmetric about the horizontal axis and the rate of change of the slope of those zones differs in direction away from peak height. This is accentuated in this case with a much more gradual decline in slope toward the bottom portion of the lens. Slope changes were obtained using a log-normal Gauss distribution function for describing angularly the thickness change. The final value of the merit function was 0.86. Improvement on the merit function was of the order of 30%. Rotation declined moderately relative to the initial stabilization design. Beginning at blink 6 less than 30 degrees rotation were seen with about 10 degrees of rotation seen at blink 12 and no rotation from blink 16 onward compared to about 38-30-15 degrees of rotation seen in the initial design over the same range of blinks. Centration was less than 0.08 mm on blink 1 and less than 0.07 at blink 4. Thereafter it dropped sharply, being less than 0.05 at 8 blinks and 0.04 at blink 16 relative to 0.06 to greater than 0.07 and 0.08 for the initial design over the same number of blink cycles. This example shows improved rotation, centration and stability as compared with the lens of Example 1.
FIG. 8 summarizes the rotation speed versus lens orientation on eye for the example 1, 2, 3, and 4. The initial design described in example 1 has an average rotation speed of about −0.55°/sec. in the 45°-0° miss-orientation range, while designs given in examples 2, 3 and 4 have an average rotation speed above −0.70°/sec. within the same miss-orientation range. Example 2 and 4 have higher rotation speed for miss-orientations below 15°. Both designs are more adequate for lenses requiring single orientation on eye such as soft contact lenses designed for high order aberration correction. These designs might require different fitting method requiring special fiducials on the front surface to help the patient for lens insertion. As the lens orientation on eye is unique due to the asymmetry of the stabilization and because of the markings on the front surface, the orientation of the lens during insertion should be very close to the final orientation of the lens after it reaches its resting position. High rotation speed for small miss-orientation at insertion will provide faster full vision correction. Those designs also present better centration performance over the design of example 3. Lens centration gets stable over a smaller number of blinks.

## Claims

1. A method for stabilizing contact lenses comprising: a) providing a lens design with a initial set of stabilization zone parameters, b) applying a merit function to the lens design based on balancing moments of momentum, and c) creating a contact lens design with improved stabilization based on the application of said merit function to the lens design with a initial set of stabilization zone parameters; wherein said merit function is selected from the group consisting of


MF1=√{square root over (WR(Rot/RREF)2+WC(Cent/CREF)2)}{square root over (WR(Rot/RREF)2+WC(Cent/CREF)2)}  (Equation 1)



MF2=√{square root over (WX(Xrange/XREF)2+WY(Yrange/YREF)2+Wθ(θrange/θREF)2)}{square root over (WX(Xrange/XREF)2+WY(Yrange/YREF)2+Wθ(θrange/θREF)2)}{square root over (WX(Xrange/XREF)2+WY(Yrange/YREF)2+Wθ(θrange/θREF)2)}  (Equation 2)



and



MF3=√{square root over ((Rod/RREF)2+WC(Cent/CREF)2+WS(Stab/SREF)2)}{square root over ((Rod/RREF)2+WC(Cent/CREF)2+WS(Stab/SREF)2)}{square root over ((Rod/RREF)2+WC(Cent/CREF)2+WS(Stab/SREF)2)}  (Equation 3)


and wherein; WR and WC are weighting factors and can take values between 0 and 1; XRange, YRange and θRange are lens performance variables in stability in horizontal direction vertical direction and rotation of the completed lens design; XREF, YREF and θREF are lens REF and REF are lens performance variables in stability in horizontal direction, vertical direction and rotation of the initial lens design; WX, WY, WS and Wθ are weighting factors. Rot Cent and Stab are lens performance variables in rotation centration and stability of the completed lens design; RREF, CREF and SREF are lens performance variables in rotation centration and stability of the initial lens design, and RREF, CREF and SREF are weighting factors.

2. The method of claim 1 wherein steps b and c are conducted iteratively.

3. The method of claim 2 wherein a virtual model that simulates the effects of eye mechanics is used to determine the contact lens designs.

4. The method of claim 3 wherein blinking is one of the eye mechanics and the simulated effects of blinking is used to adjust the stabiliztaion.

5. A contact lens made according to the method of claim 1.

6. A contact lens made according to the method of claim 2.

7. A contact lens made according to the method of claim 3.