Skip to main content

Human cornea thermo-viscoelastic behavior modelling using standard linear solid model

Abstract

Background

Corneal biomechanics is of great interest to researchers recently. Clinical findings relate them to corneal diseases and to outcomes of refractive surgery. To have a solid understanding of corneal diseases’ progression, it is important to understand corneal biomechanics. Also, they are essential for better explaining outcomes of refractive surgeries and their undesired consequences. There is a difficulty for studying corneal biomechanics in-vivo and multiple limitations arise for ex-vivo studies. Hence mathematical modelling is considered as a proper solution to overcome such obstacles. Mathematical modelling of cornea in-vivo allows studying corneal viscoelasticity with taking into consideration all boundary conditions existing in real in-vivo situation.

Methods

Three mathematical models are used to simulate corneal viscoelasticity and thermal behavior in two different loading situations: constant and transient loading. Two models of the three are used for viscoelasticity simulation which are Kelvin-Voigt and standard linear solid models. Also, temperature rise due to the ultrasound pressure push is calculated using bioheat transfer model for both the axial direction and as a 2D spatial map using the third model (standard linear solid model).

Results

Viscoelasticity simulation results show that standard linear solid model is efficient for describing the viscoelastic behavior of human cornea in both loading conditions. Results show also that the deformation amplitude obtained from standard linear solid model is more reasonable for corneal soft-tissue deformation with respect to corresponding clinical findings than that obtained from Kelvin-Voigt model. Thermal behavior results estimated corneal temperature rise to be roughly 0.2 °C, which conforms with FDA regulations for soft tissue safety.

Conclusion

Standard Linear Solid (SLS) model is better describing the human corneal behavior in response to constant and transient load more efficiently. Temperature rise (TR) for the corneal tissue of about 0.2 °C is conforming with FDA regulations and even less than the FDA regulations for soft tissue safety.

Peer Review reports

Introduction

Deformation behavior of soft tissues due to stress is investigated by measuring their biomechanics properties. Recently, research is giving attention to biomechanical properties of ocular tissues so as to understand the progression of ocular diseases such as: keratoconus, and post-refractive ectasia which are characterized by corneal tissue deformation. Also, elasticity changes before and after the refractive surgeries are of their interest as they affect the developing of post-refractive ectasia. Biomechanical properties are reported to be changing due to keratoconous corneas and refractive surgeries [1], [2]. A comprehensive review about the impact of post-refractive surgeries on corneal biomechanics and how can they be measured is presented by our previous paper [3].

Understanding of cornea biomechanics leads to better understanding of these corneal alterations and how to treat them. Many researches have focused on investigation of the corneal biomechanics ex-vivo such as the work reported in [4]. However, ex-vivo cornea experience swelling and loss of its tear film leading to changes in its biomechanical properties [1]. Hence, there is a great need for in-vivo cornea studies that allows for normal and diseased cornea biomechanics. Mathematical modelling of corneal biomechanics is the first step in understanding its behavior in normal and diseased states. A suitable model should describe the elastic and viscous components of the cornea as its tissue exhibits viscoelastic behavior when loaded under certain forces transiently [2].

Human cornea experiences a viscoelastic behavior when subjected to a transient stress [5]. This behavior can be represented by two components; elastic component and viscous component. The elastic component gives an instantaneous deformation while the viscous component gives a damping deformation. Elastic and viscous components can be modelled by spring; with elasticity of (E); and dashpot; with viscosity of (\(\eta\)); system components respectively. The spring represents the pure elastic behavior to an applied load, while the dashpot represents the time-dependent viscous resistance to that load.

There are multiple models that are used to describe the viscoelasticity behavior, amongst them Maxwell model, Kelvin-Voigt model and the standard linear solid model [6]. The three models are presented in Fig. 1. Maxwell’s model consists of a spring and a dashpot connected in series to each other, while Kelvin-Voigt model consists of the same two components but connected in parallel to each other. Standard linear solid model in Kelvin-Voigt representation has an extra spring element connected in series to the parallel spring-dashpot connection.

Fig. 1
figure 1

(a) Maxwell Model, (b) Kelvin-Voigt and (c) Standard linear solid model in Kelvin-Voigt representation

Several research proposed mathematical models for studying corneal biomechanics. One of them is the model reported in [5], where two Kelvin-Voigt models in series with a spring is proposed to compare the effect of fast versus slow viscoelastic behavior. In [7], a Kelvin-Voigt model is proposed to study the effect of elasticity and viscosity separately on corneal hysteresis. In [8], a mathematical model based simulation is proposed to investigate induced corneal vibrations due to air-puff and its relationship to intra-ocular pressure and viscoelasticity of the cornea. Other mathematical models are proposed to study the heating effect of different wavelengths of UV lasers used in corneal ablation procedures as well [9,10,11]. Temperature rise calculation due to microwave radiation formed the first step towards modelling the thermal behavior of the human eye [12], [13]. While the first FEM model developed for describing the human eye thermal behavior while undergoing LASER radiation is introduced by [14].

This paper presents a mathematical model to investigate the corneal tissue viscoelastic behavior and thermal behavior, i.e., temperature rise (TR); when subjected to a transient force under the effect of intra-ocular pressure and tear film pressure. This study represents a continuation of our project that is presented in [15].

The paper is organized as follows, Sect. 1 presents an introduction about the importance of corneal biomechanics to corneal diseases and refractive surgeries outcomes, different mathematical models that are used to describe soft-tissue biomechanics and thermal behavior and advantages of each one with respect to the others. The proposed methodology for the simulation procedure using both models is introduced in Sect. 2. Results are presented and discussed in Sect. 3 and the paper is concluded in Sect. 4.

Methodology

The viscoelastic and the thermal behaviors of the corneal tissue are simulated using mathematical equations to investigate the corneal behavior in response to ultrasonic transient load. The viscoelastic behavior is simulated using a constant and a transient load simulation. While the thermal behavior is simulated in case of a transient load only. The methodology is divided into two major sections, one for simulating the viscoelastic behavior and the other for simulating the thermal behavior. The first major section is furthermore divided into two subsections, one for simulating the corneal viscoelastic behavior due to constant load and the other for simulating its behavior due to the transient load. Afterwards; the corneal viscoelastic behavior is simulated for a wide range of Young’s moduli starting from 150 KPa till 3 MPa [16], [17] in order to predict whether the corneal behavior is linear or not. This is in correspondence with Fung et al. [18], [19] and Hoeltzel et al. [20], where they recommend that for nonlinear behavior tissues, the study should not be carried out for a single Young’s modulus.

Viscoelastic behavior

In this section, both the Kelvin-Voigt and the standard linear solid models are used to simulate the human cornea viscoelastic behavior. Maxwell’s model is discarded due to its indefinite creep behavior and because it does not fully recover from deformation due to the pure viscous component.

The complete vertical cross-section of human eye with cornea represented as standard linear solid model is shown in Fig. 2. The corneal tissue is modelled using Kelvin-Voigt model and standard linear solid model. Its behavior using both models is simulated under two scenarios; while applying a constant load and a transient load respectively. The constant load is simulated as a unit step function, and the transient load is simulated as a sinusoidal burst of 3 cycles applied for 40% of the simulation time. The transient load is chosen to be a sinusoidal burst as it simulates the acoustic radiation force burst as an internal actuator generated from an ultrasound transducer during elastography procedure [3], [15], [21].

Fig. 2
figure 2

Complete vertical cross-section of human eye, and corneal tissue represented with standard linear solid model

The corneal tissue elasticity is set to 150 KPa to simulate a normal corneal tissue [22]. Either the constant load or the transient load is applied at the corneal apex. The simulation time is 5 msec for both Kelvin-Voigt and standard linear solid model to study their strain behavior. Constant load is applied for the full 5 msec, while the transient load is applied for only 2 msec. The complete simulation parameters are listed in Table 1.

For a Kelvin-Voigt material; the stress-strain equation as a function of time is given by Eq. (1):

$$\sigma \left(t\right)=E\epsilon \left(t\right)+ \eta \frac{d\epsilon \left(t\right)}{dt}$$
(1)

where \(\sigma\) is the applied load in (Pa), \(E\) is the modulus of elasticity, \(\epsilon\) is the strain in (mm) and \(\eta\) is the viscosity in (Pa.sec.).

The strain as a function of time due to Kelvin-Voigt model will follow Eq. (2) when subjected to a sudden stress:

$$\epsilon \left(t\right)=\frac{{\sigma }_{o}}{E} (1- {e}^{\frac{-t}{{\tau }_{R}}})$$
(2)

where \({\tau }_{R}= \frac{\eta }{E}\) is the retardation time and t is the simulation time in seconds.

Strain as a function of time due to standard linear solid model is given in Eq. (3):

$$\epsilon \left(t\right)=\frac{{\sigma }_{o}}{E} \left(1- {e}^{\frac{-t}{{\tau }_{R}}}\right)\left({e}^{\frac{-t}{{\tau }_{R}}}\right)$$
(3)

When applying either a constant or a transient load Papplied; the total load experienced by the cornea is given by the following Eq. (4):

$${P}_{resultant}={P}_{applied} + {P}_{tear film}- IOP$$
(4)

where Ptear film is the tear film pressure working in the same direction of the applied load and equals 4.15 mmHg [23], and IOP is the intra-ocular pressure. Tear film pressure is responsible of maintaining corneal tissue constituents in place together.

Table 1 Simulation parameters for Kelvin-Voigt and Standard Linear Solid models

Thermal behavior

Temperature rise (TR)

In this section, we will consider a focused ultrasound transducer placed at 28 mm above the corneal tissue of the human eye. A saline cup is placed in the distance between the transducer face till the corneal tissue. Thermal behavior is simulated using the mathematical model of bioheat transfer equation formulated by Pennes [24], [25]. Bioheat equation is used in a former work presented by [24] to study the thermal behavior of human eye due to InfraRed (IR) LASER exposure in ophthalmic refractive procedures. Bioheat transfer for corneal tissue is given by Eq. (5):

$${\rho }_{t}{C}_{t}\frac{\partial T}{\partial t}= \nabla \left(k\nabla T\right)+ {W}_{b}{C}_{pb}\left({T}_{a}-T\right)+ {Q}_{m}+H$$
(5)

Where \({\rho }_{t}\) represents tissue density, \({C}_{t}\) represents tissue specific heat, \({W}_{b}\) represents the blood perfusion rate, \({C}_{pb}\) represents the specific heat of blood, \({T}_{a}\) represents the arterial blood temperature, \(k\) represents the specific tissue thermal conductivity, and \(\nabla T\) represents the temperature spatial gradient. The heat generated inside the tissue in response to an external heat source due to temperature rise (TR) is represented by the term H. we will model H mathematically in the next section. The external heat source in this study is the pressure push due to ultrasonic transient load. The boundary condition equation of the cornea is given by Eq. (6) [24], [26]:

$$-k\frac{\partial T}{\partial x}= {h}_{c}\left(T- {T}_{amb}\right)+ \sigma \left({T}^{4}- {T}_{amb}^{4}\right)$$
(6)

Where hc represents the heat transfer coefficient of the cornea. Stefan-Boltzmann constant is represented by σ. The emissivity of the corneal surface is represented by ϵ. \({T}_{amb}\) represents the ambient temperature surrounding the anterior cornea. Equation 6 represents the convection-type thermal exchange between the cornea and the ambient. Bioheat equation parameters’ values are given in Table 2.

Table 2 Bioheat parameters’ values for Eqs. (5) & (6) [24]

The temperature rise (TR) rate of change is represented as the sum of all heat effects, such as, the metabolic heat effect, the cooling down / heating up effect due to arterial blood flow, and the net heat conduction-type effect in the human tissue.

Pressure wave modelling

Hence, we need to model the ultrasound transient push that is imparted to the corneal tissue mathematically in terms of its pressure function. Then by substituting with this imparted pressure to the medium in Eq. (5), we solve for the temperature rise (TR) for the corneal tissue due to this imparted transient push. The ultrasound pressure push is modelled mathematically for a circular planar piston by Rayleigh-Sommerfeld integral [27], [28] and given by Eq. (7):

$$p(z, w)=\rho c{v}_{0}({e}^{\left(ikz\right)}- {e}^{ik\sqrt{({z}^{2}- {a}^{2})}})$$
(7)

Where a is the radius of the circular piston transducer, v0 is the uniform velocity field over the entire aperture, k is the wave number, and \(\rho\) is the density. The equation states that the pressure wave generated is due to two components the direct wave component which is equal to \({e}^{\left(ikz\right)}\), and the edge wave component which is equal to \({e}^{ik\sqrt{({z}^{2}- {a}^{2})}}\). Both components form the on-axis pressure wave. However, for a focused ultrasound transducer the equation can be reduced to Eq. (8). Focused ultrasound transducer can be visualized by means of Fig. 3, where the spherical aperture is denoted by the surface S, the circular aperture is of radius a, and the spherical surface is of radius R0. Now the integral can be reduced to the following equation:

$$p(z, w)=\frac{\rho c{v}_{0}}{{q}_{0}}({e}^{\left(ikz\right)}- {e}^{ik{r}_{e}})$$
(8)

Where q0 equals to 1-z/R0, and re equals to \(\sqrt{({(z-h)}^{2}- {a}^{2})}\).

Fig. 3
figure 3

Spherically focused ultrasound transducer

The resulting pressure wave that is imparted to the corneal tissue is represented graphically by Fig. 4. The resulting pressure wave is fed into Eq. (5) in order to obtain the temperature rise (TR) due to its effect. The pressure wave distribution presented by Fig. 4 is for both the axial line in front of the transducer and the average pressure wave for lateral positions around the transducer’s center.

Fig. 4
figure 4

Resulting on-axis pressure wave

Results and discussion

Simulation is performed using MATLAB program (R2020b) and using a machine of 8 GB RAM and core i5 processor with a best value of execution time of 3.3 msec for transient load scenario and 4 msec for constant load scenario.

Constant load

For the Kelvin-Voigt model, when a constant load is applied, it starts to deform instantaneously till it reaches an asymptote. The model describes the instantaneous (elastic) deformation due to sudden load, but it fails in describing the viscous damping deformation effect of the viscoelastic material.

The applied constant load is shown in Fig. 5(a). Kelvin-Voigt model time dependent strain (deformation) behavior due to applying constant load is presented in Fig. 5(b). The material does not recover from its deformation, but it continues to creep constantly. The deformation amplitude due to the applied constant load is about 4.4 mm. The material reaches its peak deformation amplitude at 3 msec approximately.

Fig. 5
figure 5

(a) Constant load, (b) Kelvin-Voigt model time dependent strain behavior and (c) Standard linear solid model time dependent strain behavior

The standard linear solid model is able to describe both the instantaneous (elastic) deformation and the viscous damping deformation of the material. The time dependent strain behavior due to constant load is shown in Fig. 5(c). The deformation amplitude due to this constant load is about 1.1 mm. The material reaches its peak deformation amplitude at 0.45 msec approximately. The material returns to its original state recovering from its deformation.

From the constant load simulation, it is clear that the standard linear solid model is able to describe the corneal tissue viscoelasticity more effectively. This is because it models the material recovery from deformation and because it predicts more reasonable deformation amplitude value compared to Kelvin-Voigt model.

Transient load

The applied transient sinusoidal load burst is presented in Fig. 6(a). Kelvin-Voigt model yields a sinusoidal strain behavior that increases with time. The model describes again the instantaneous elastic deformation effectively. However, it fails in describing the viscous damping behavior. The deformation amplitude reaches 0.007 m approximately and continues to increase with each new cycle of excitation which is not logical for corneal tissue. This behavior is presented in Fig. 6(b).

Fig. 6
figure 6

(a) Sinusoidal load burst, (b) Kelvin-Voigt time dependent strain behavior to transient load and (c) standard linear solid model in Kelvin-Voigt time dependent strain behavior to transient load

Standard linear solid model behavior is presented in Fig. 6(c). The model predicts the material to have an instantaneous elastic deformation combined with exponential damping for deformation behavior. The deformation amplitude is approximately 1.2 mm which decreases significantly with each cycle of load excitation till reaches original state. The deformation amplitude is reasonable for cornea tissue in response to the applied transient load compared to the deformation amplitude predicted by Kelvin-Voigt model.

In terms of temporal localization, where the deformation is supposed to be highly localized within short temporal duration, the predicted behavior due to Kelvin-Voigt model is not suitable as it does not fulfill the temporal localization for the corneal deformation. On the other hand, standard linear solid model predicts properly the temporal localization of the corneal deformation.

The deformation behavior is calculated for different multiple elastic moduli starting from 150KPa to 3 MPa to investigate the mathematical model behavior for a wide range of elastic moduli using SLS model. This is done in order to understand the behavior in both pre-refractive and post-refractive states of the cornea. The complete behavior of the mathematical model for transient load is presented in Fig. 7. It is observed that the behavior is a non-linear behavior.

Fig. 7
figure 7

Generalized mathematical model behavior for transient load

Temperature rise (TR)

The temperature rise (TR) along the axial direction of the transducer is presented by Fig. 8. Temperature rise (TR) is at its maximum near the ultrasound transducer with a ∆T of 0.9 °C. while at the corneal tissue the temperature rise is about 0.17 °C. In addition, corneal thermal behavior is observed to be exponentially decaying with two decaying rates. The first decaying rate is found in the region between the transducer and the surface of the corneal tissue, where this decay rate is slow. The second decaying rate is faster than the first one and found in the region where the corneal tissue is involved and the rest of the human eye structure. Second thermal decaying rate is observed to be faster than the first one because the corneal tissue absorbs much of the imparted thermal energy than the saline cup does. Also, faster decaying rate is due to the fluid-like structures in the human eye such as; aqueous humor in the anterior chamber and the vitreous humor in the posterior chamber, which both act as heat sinks for the imparted thermal energy.

Fig. 8
figure 8

Temperature Rise (TR) in Celsius for the axial line extending from the transducer face till the posterior part of the human eye

The 2D distribution for the pressure wave generated by the ultrasound transducer is calculated as well and fed to the bioheat equation in order to calculate the 2D temperature rise around both the corneal tissue; anteriorly and posteriorly; and inside the corneal tissue itself. The 2D spatial map for the temperature rise due to the pushing is calculated using bio-heat transfer equation [29] and is presented in Fig. 9.

Fig. 9
figure 9

Temperature rise 2D spatial map due to pushing beam

In terms of temperature rise due to pushing beam, the increase in temperature is concentrated at the proximal area of the transducer forming a delta shaped area of temperature of about 0.9 °C. Hence, the temperature rise is concentrated in front of the transducer directly. From literature, there is a debate about the acceptable temperature rise due to ultrasound waves. As in [30], a temperature rise from 1 to 2 °C is acceptable for all body organs. However, tissue temperature may increase by as much as 10 °C, raising safety concerns even though the acoustic output was still within the Food and Drug Administration’s recommended maximum output exposure level for diagnostic ultrasound [31]. Moreover, The FDA regulates the temperature rise (TR) of tissue to be no more than 6 °C as stated by [32]. In general, the anterior chamber contains a watery-like fluid called aqueous humor, and the posterior chamber contains a gel-like fluid called vitreous humor. Fluids inside the human eye either in the anterior chamber or the posterior chamber act as a heat sink for any generated heat. By consequence, the thermal effect induced due to the acoustic radiation force impulse is reduced by these fluids.

Conclusion

In this paper, two mathematical models are proposed to simulate corneal viscoelastic behavior. Their performance is compared when applying a constant and transient load respectively. The simulation time is 5 msec for both models. The constant load is applied for the whole simulation time, while the transient load is applied for about 2 msec.

From the obtained results of the two proposed models and with comparison of these two models with each other, it is concluded that standard linear solid model is describing the human corneal behavior in response to constant and transient load more efficiently. This is because of the lack of Kelvin-Voigt model in describing the viscous damping behavior of corneal biomechanics. Also, because the deformation amplitude in Kelvin-Voigt is greater than the reasonable values for the corneal tissue for the same applied load either constantly or transiently.

On the other hand, standard linear solid model predicts a reasonable deformation amplitude for corneal tissue for the same applied load either constantly or transiently. This model also describes the viscous damping behavior of the corneal tissue efficiently.

Also, the thermal behavior of the corneal tissue is calculated using bioheat transfer equation. Temperature rise for corneal due ultrasound pressure push is calculated.

Both the axial and 2D spatial distributions of the temperature rise is calculated for the corneal tissue. Maximum temperature rise (TR) of about 0.9 °C is found to be near to the transducer face. However, the temperature rise for the corneal tissue is found to be about 0.2 °C, which is less the FDA regulations for soft tissue safety.

Study limitations

We used the fundamental models in biomechanics as a first step in modelling corneal biomechanics. We did not include other advanced models such as Prony and Kevin series with more model parameters for this study. However, we intend to include them in our future improvement of the presented mathematical model to improve the results.

In addition, we calculated the temperature rise for one certain ultrasonic push. Yet, in the future improvement we intend to calculate the temperature rise for a range of ultrasonic pushes in order to determine the maximum push after which the temperature rise inside the cornea will exceed the threshold recommended by FDA. This will guide the manufacturing specifications of the ultrasound probes.

Data availability

The datasets presented and analyzed during the current study are available from the corresponding author upon reasonable request.

References

  1. Edmund C. Corneal topography and elasticity in normal and keratoconic eyes. A methodological study concerning the pathogenesis of keratoconus. ” Acta Ophthalmol Suppl (Oxf). 1989;193:1–36.

    CAS  Google Scholar 

  2. Edmund C. Corneal elasticity and ocular rigidity in normal and keratoconic eyes. Acta Ophthalmol. 1988;66(2):134–40.

    Article  CAS  Google Scholar 

  3. Ahmed HM, Salem NM, Al-Atabany W. Impact of post–refractive surgeries on corneal Biomechanics—A review. J Clin Eng. 2021;46(4):177–88.

    Article  Google Scholar 

  4. Tanter M, Touboul D, Gennisson J, Bercoff J, Fink M. “High-Resolution Quantitative Imaging of Cornea Elasticity Using Supersonic Shear Imaging,” IEEE Trans. Med. Imaging, vol. 28, no. 12, pp. 1881–1893, Dec. 2009, doi: https://doi.org/10.1109/TMI.2009.2021471.

  5. Kobayashi AS, Staberg LG, Schlegel WA. Viscoelastic properties of human cornea. Exp Mech. 1973;13(12):497–503.

    Article  Google Scholar 

  6. Lakes RS. Viscoelastic solids. Volume 9. CRC press; 1998.

  7. Glass DH, Roberts CJ, Litsky AS, Weber PA. A viscoelastic biomechanical model of the cornea describing the effect of viscosity and elasticity on hysteresis. Investig Ophthalmol \& Vis Sci. 2008;49(9):3919–26.

    Article  Google Scholar 

  8. Han Z, et al. Air puff induced corneal vibrations: theoretical simulations and clinical observations. J Refract Surg. 2014;30(3):208–13.

    Article  PubMed  Google Scholar 

  9. Abdelhalim I, Hamdy O, Hassan AA, Elnaby SH. Dependence of the heating effect on tissue absorption coefficient during corneal reshaping using different UV lasers: a numerical study. Phys Eng Sci Med. 2021;44(1):221–7.

    Article  PubMed  Google Scholar 

  10. Rahbar S, Shokooh-Saremi M. Mathematical modeling of laser linear thermal effects on the anterior layer of the human eye. Opt \& Laser Technol. 2018;99:72–80.

    Article  CAS  Google Scholar 

  11. Gokul KC, Gurung DB, Adhikary PR. “Mathematical Model: comparative study of Thermal Effects of laser in corneal refractive surgeries.,” Appl \& Appl Math, vol. 10, no. 1, 2015.

  12. Emery AF, Kramar P, Guy AW, Lin JC. “Microwave induced temperature rises in rabbit eyes in cataract research,” 1975.

  13. Taflove A, Brodwin ME. Computation of the electromagnetic fields and induced temperatures within a model of the microwave-irradiated human eye. IEEE Trans Microw Theory Tech. 1975;23(11):888–96.

    Article  Google Scholar 

  14. Amara EH. Numerical investigations on thermal effects of laser-ocular media interaction. Int J Heat Mass Transf. 1995;38(13):2479–88.

    Article  Google Scholar 

  15. Ahmed HM, Al-Atabany W, Salem NM. Corneal Biomechanics Assessment using high frequency Ultrasound B-Mode Imaging. IEEE Access. 2021;9:106014–28. https://doi.org/10.1109/ACCESS.2021.3100098.

    Article  Google Scholar 

  16. 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(1):11–9. https://doi.org/10.1080/02713680601077145.

    Article  PubMed  Google Scholar 

  17. Elsheikh A, Wang D, Pye D. “Determination of the modulus of elasticity of the human cornea,” J. Refract. Surg, vol. 23, no. 8, pp. 808–818, 2007, doi: https://doi.org/10.3928/1081-597x-20071001-11.

  18. Fung YC, Skalak R. “Biomechanics: mechanical properties of living tissues,” 1981.

  19. Fung YC. Elongation of Soft Biological Tissues in simple elongation. Am J Physiol. 1987;213(6):1532–43.

    Article  Google Scholar 

  20. Hoeltzel DA, Altman P, Buzard K, Il Choe K. Strip extensiometry for comparison of the mechanical response of bovine, rabbit, and human corneas. J Biomech Eng. 1992;114(2):202–15. https://doi.org/10.1115/1.2891373.

    Article  CAS  PubMed  Google Scholar 

  21. H. M., N. M., A. F., and M. I., “On Shear Wave Speed Estimation for Agar-Gelatine Phantom,” Int J Adv Comput Sci Appl, vol. 7, no. 2, pp. 401–9, 2016, doi: https://doi.org/10.14569/ijacsa.2016.070255.

  22. Ambrósio R, Alonso RS, Luz A, Coca Velarde LG. Corneal-thickness spatial profile and corneal-volume distribution: tomographic indices to detect keratoconus. J Cataract Refract Surg. 2006;32(11):1851–9. https://doi.org/10.1016/j.jcrs.2006.06.025.

    Article  PubMed  Google Scholar 

  23. Damji KF, Muni RH, Munger RM. Influence of corneal variables on accuracy of intraocular pressure measurement. J Glaucoma. 2003;12(1):69–80.

    Article  PubMed  Google Scholar 

  24. Cvetkovic M, Peratta A, Poljak D. Thermal modelling of the human eye exposed to infrared radiation of 1064 nm nd: YAG and 2090 nm Ho: YAG lasers. Environ Heal Risk. 2009;14:221–31.

    CAS  Google Scholar 

  25. Pennes HH. Analysis of tissue and arterial blood temperatures in the resting human forearm. J Appl Physiol. 1948;1(2):93–122.

    Article  CAS  PubMed  Google Scholar 

  26. Cvetkovi M, Poljak D, Peratta A, others. 3D FEM temperature distribution analysis of the human eye exposed to laser radiation. WIT Trans Eng Sci. 2010;68:303–12.

    Article  Google Scholar 

  27. Farjat AE, Etcheverry JI. “Mathematical modeling of the radiated acoustic field of ultrasonic transducers,” J. Nondestruct. Test. Ultrason, vol. 2007, no. 11, pp. 1–12, 2007.

  28. Cawley P, Cegla FB. “Ultrasonic non-destructive testing.” Google Pat, 2016.

  29. Scott JA. A finite element model of heat transport in the human eye. Phys Med \& Biol. 1988;33(2):227.

    Article  CAS  Google Scholar 

  30. Soltan A, et al. A head mounted device stimulator for optogenetic retinal prosthesis. J Neural Eng. 2018;15(6):65002.

    Article  Google Scholar 

  31. Liu Y, Herman BA, Soneson JE, Harris GR. Thermal safety simulations of transient temperature rise during acoustic radiation force-based ultrasound elastography. Ultrasound Med \& Biol. 2014;40(5):1001–14.

    Article  Google Scholar 

  32. Song P, Zhao H, Manduca A, Urban MW, Greenleaf JF, Chen S. Comb-push ultrasound shear elastography (CUSE): a novel method for two-dimensional shear elasticity imaging of soft tissues. IEEE Trans Med Imaging. 2012;31(9):1821–32.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not Applicable.

Funding

Not Applicable.

Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB).

Author information

Authors and Affiliations

Authors

Contributions

Hassan M. Ahmed was involved in conceptualization, investigation and writing original draft preparation.Nancy M. Salem was involved in data curation, visualization, reviewing and editing.Walid Al-Atabany was involved in data curation, visualization, reviewing and editing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Hassan M. Ahmed.

Ethics declarations

Ethics approval and consent to participate

Not Applicable.

Consent for publication

Not Applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ahmed, H.M., Salem, N.M. & Al-Atabany, W. Human cornea thermo-viscoelastic behavior modelling using standard linear solid model. BMC Ophthalmol 23, 250 (2023). https://doi.org/10.1186/s12886-023-02985-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12886-023-02985-3

Keywords