A viscoelastic biomechanical model of the cornea describing the effect of viscosity and elasticity on hysteresis

PURPOSE. To develop a method for evaluating viscosity and elasticity of the cornea and to examine the effect that both properties have on hysteresis. METHODS. A three-component spring and dashpot model was created in Simulink in Matlab to represent the purely elastic and viscoelastic behavior of the cornea during a measurement using device called an ocular response analyzer (ORA). Values for elasticity and viscosity were varied while sinusoidal stress was applied to the model. The simulated stresses were used to determine how hysteresis is affected by the individual components of elasticity, viscosity, and maximum stress. To validate the model, high-speed photography was used to measure induced strain in a corneal phantom during ORA measurement. This measured strain was compared with the strains simulated by the model. RESULTS. When the spring in the viscoelastic portion of the model was stiffened, hysteresis decreased. When the spring in the purely elastic element was stiffened, hysteresis increased. If both springs were stiffened together, hysteresis peaked strongly as a function of the viscosity of the viscoelastic element. Below the peak value, lower elasticity was associated with higher hysteresis. Above the peak value, higher elasticity was associated with higher hysteresis. In addition, hysteresis increased as the air maximum pressure was increased. Measurements from phantom corresponded to predictions from the model. CONCLUSIONS. A viscoelastic model is presented to illustrate how changing viscosity and elasticity may affect hysteresis. Low hysteresis can be associated with either high elasticity or low elasticity, depending on the viscosity, a finding consistent with clinical reports.