Research article

Biomechanical responses of the cornea after small incision lenticule extraction (SMILE) refractive surgery based on a finite element model of the human eye


  • Received: 20 March 2021 Accepted: 12 May 2021 Published: 17 May 2021
  • Purpose 

    To investigate the biomechanical responses of the human cornea after small incision lenticule extraction (SMILE) procedures, especially their effects of SMILE surgery on stress and strain.

    Methods 

    Based on finite element analysis, a three-dimensional (3D) model of the human eye was established to simulate SMILE refractive surgery procedures. Stress and strain values were calculated by inputting the intraocular pressure (IOP).

    Results 

    After SMILE refractive surgery procedures, the stress and strain of the anterior and posterior corneal surfaces were significantly increased. The equivalent stress and strain on the anterior and posterior corneal surfaces increased with increasing diopter and were concentrated in the central area, whereas the values of stress and strain at the incision site on the anterior surface of the cornea were approximately 0. Compared with the anterior corneal surface, the stress and strain of the posterior surface were larger. Increasing IOP caused an approximately linear change in stress and a nonlinear increase in corneal strain. In addition, we found that the incision sizes and direction had less of an influence on stress and strain. In summary, SMILE surgery increased the equivalent stress and strain on the human cornea.

    Conclusions 

    The equivalent stress and strain of the anterior and posterior human corneal surfaces increased after SMILE refractive surgery; these increases were particularly noticeable on the posterior surface of the cornea.

    Citation: Yinyu Song, Lihua Fang, Qinyue Zhu, Ruirui Du, Binhui Guo, Jiahui Gong, Jixia Huang. Biomechanical responses of the cornea after small incision lenticule extraction (SMILE) refractive surgery based on a finite element model of the human eye[J]. Mathematical Biosciences and Engineering, 2021, 18(4): 4212-4225. doi: 10.3934/mbe.2021211

    Related Papers:

    [1] Ruirui Du, Lihua Fang, Binhui Guo, Yinyu Song, Huirong Xiao, Xinliang Xu, Xingdao He . Simulated biomechanical effect of aspheric transition zone ablation profiles after conventional hyperopia refractive surgery. Mathematical Biosciences and Engineering, 2021, 18(3): 2442-2454. doi: 10.3934/mbe.2021124
    [2] Yan Li, Kuan Wang, Lejun Wang, Tongbo Chang, Shengnian Zhang, Wenxin Niu . Biomechanical analysis of the meniscus and cartilage of the knee during a typical Tai Chi movement—brush-knee and twist-step. Mathematical Biosciences and Engineering, 2019, 16(2): 898-908. doi: 10.3934/mbe.2019042
    [3] Fang Wang, Jiaming Wang, Mingxin Li, Jun Hu, Kehua Song, Jianguo Zhang, Yubo Fan . Biomechanical study of the effect of traction on elbow joint capsule contracture. Mathematical Biosciences and Engineering, 2023, 20(12): 21451-21466. doi: 10.3934/mbe.2023949
    [4] Niksa Praljak, Shawn D. Ryan, Andrew Resnick . Pulsatile flow through idealized renal tubules: Fluid-structure interaction and dynamic pathologies. Mathematical Biosciences and Engineering, 2020, 17(2): 1787-1807. doi: 10.3934/mbe.2020094
    [5] Debao Guan, Yuqian Mei, Lijian Xu, Li Cai, Xiaoyu Luo, Hao Gao . Effects of dispersed fibres in myocardial mechanics, Part I: passive response. Mathematical Biosciences and Engineering, 2022, 19(4): 3972-3993. doi: 10.3934/mbe.2022183
    [6] Xu Bie, Yuanyuan Tang, Ming Zhao, Yingxi Liu, Shen Yu, Dong Sun, Jing Liu, Ying Wang, Jianing Zhang, Xiuzhen Sun . Pilot study of pressure-flow properties in a numerical model of the middle ear. Mathematical Biosciences and Engineering, 2020, 17(3): 2418-2431. doi: 10.3934/mbe.2020131
    [7] Christophe Prud'homme, Lorenzo Sala, Marcela Szopos . Uncertainty propagation and sensitivity analysis: results from the Ocular Mathematical Virtual Simulator. Mathematical Biosciences and Engineering, 2021, 18(3): 2010-2032. doi: 10.3934/mbe.2021105
    [8] Xiaoyu Du, Yijun Zhou, Lingzhen Li, Cecilia Persson, Stephen J. Ferguson . The porous cantilever beam as a model for spinal implants: Experimental, analytical and finite element analysis of dynamic properties. Mathematical Biosciences and Engineering, 2023, 20(4): 6273-6293. doi: 10.3934/mbe.2023270
    [9] Longting Ding, Yuan Li, Zhanchuang Han, Mengyuan Zhang, Xuancang Wang, Lu He . The effect of RAP content on fatigue damage property of hot reclaimed asphalt mixtures. Mathematical Biosciences and Engineering, 2024, 21(2): 3037-3062. doi: 10.3934/mbe.2024135
    [10] Rui Zhu, Yuhang Chen, Qianqian Yu, Siqing Liu, Jianjie Wang, Zhili Zeng, Liming Cheng . Effects of contusion load on cervical spinal cord: A finite element study. Mathematical Biosciences and Engineering, 2020, 17(3): 2272-2283. doi: 10.3934/mbe.2020120
  • Purpose 

    To investigate the biomechanical responses of the human cornea after small incision lenticule extraction (SMILE) procedures, especially their effects of SMILE surgery on stress and strain.

    Methods 

    Based on finite element analysis, a three-dimensional (3D) model of the human eye was established to simulate SMILE refractive surgery procedures. Stress and strain values were calculated by inputting the intraocular pressure (IOP).

    Results 

    After SMILE refractive surgery procedures, the stress and strain of the anterior and posterior corneal surfaces were significantly increased. The equivalent stress and strain on the anterior and posterior corneal surfaces increased with increasing diopter and were concentrated in the central area, whereas the values of stress and strain at the incision site on the anterior surface of the cornea were approximately 0. Compared with the anterior corneal surface, the stress and strain of the posterior surface were larger. Increasing IOP caused an approximately linear change in stress and a nonlinear increase in corneal strain. In addition, we found that the incision sizes and direction had less of an influence on stress and strain. In summary, SMILE surgery increased the equivalent stress and strain on the human cornea.

    Conclusions 

    The equivalent stress and strain of the anterior and posterior human corneal surfaces increased after SMILE refractive surgery; these increases were particularly noticeable on the posterior surface of the cornea.



    The incidence of myopia is increasing yearly and has become a global public health problem [1]. Corneal refractive surgery is an one of the important way of treating refractive errors. To change the refractive ability of the eye, many refractive surgery techniques have been developed [2]. SMILE refractive surgery has emerged as a new surgical approach. It differs from Laser-assisted in situ keratectomy (LASIK) surgery because no corneal flap is created, and the corneal lenticule is extracted through a small incision. Because it is flapless and creates only a small incision, SMILE surgery is safer. The most important advantage of SMILE over LASIK is that it leaves the cornea with greater tensile strength [3], causing less damage to corneal biomechanics [4]. Therefore, SMILE has become the new choice for myopia correction procedures in traditional refractive surgery.

    Corneal refractive surgery cuts corneal tissue to correct vision, and it has influence on the biomechanics of the corneal. Gyldenkerne et al. [4] confirmed that the curvature of the anterior corneal surface changed significantly after SMILE refractive surgery. Wei et al. [5] found that refractive surgery changed corneal structure and biomechanics. The stress and strain on cornea also increased after refractive surgery. Scarcelli et al. [6] showed that tensile tests were strongly correlated with the biomechanical properties of the cornea and the depth of the matrix. Knox et al. [7] found that the corneal stress changed after the incision was created in LASIK surgery. Franus [8] also studied the stress and strain of the cornea after LASIK. In addition, the main physiological function of the cornea is refraction, accounting for 70% of the refractive power of the eye [9]. And visual quality is closely related to corneal biomechanics after refractive surgery. Therefore, it is necessary to explore the postoperative characteristics of corneal stress and strain.

    The finite element analysis method is a useful tool for calculating geometric, biomechanical and structural characteristics of biological [10]. Since the mid-20th century, its practicability has been verified in many fields. The method has been widely used in structural and fluid research and has been used to solve complex research problems [11]. Ariza-Gracia et al. [12] used the Gasser-Holzapfel-Ogden (GHO) constitutive model for finite element simulation of the human cornea. Kaliske et al. [13] described the viscoelastic material model. Anderson et al. [14] proposed the anisotropy model of the corneal shell. Pinsky et al. [15] applied the numerical model of anisotropy of corneal and scleral mechanics in surgery. In brief, the finite element analysis could accurately assess postoperative corneal biomechanics. Based on previous studies, we established a three-dimensional (3D) finite element whole-eye model to simulate SMILE refractive surgery to study postoperative changes in corneal stress and strain and provide valuable theoretical data for the diagnosis of ophthalmic diseases.

    Corneal refractive surgery is predictable and stable [16]. However, the biomechanics become weaker and the vision problems, such as corneal thinning, myopia aggravation, and visual acuity decline still occur [17]. John et al. [18] found that the biomechanical properties of corneal reduced after SMILE and LASIK surgery, and the lower myopic had more tensile strength in SMILE procedure. Osman et al. [19] also confirmed that a retrospective study showed less reduction in corneal biomechanical properties of SMILE compared with the LASIK group. Peinado et al. [20] indicated that corneal distention caused by biomechanical impacted the whole corneal structure. Randleman et al. [17] showed that refractive surgery in which the cornea was cut might reduce the biomechanical stability of the cornea, even causing the complication of iatrogenic keratectasia. Cordero-Mendieta showed that corneal biomechanics had a great impact on the ophthalmological studies. For this reason, we used finite element analysis to study the biomechanical characteristics of cornea after SMILE surgery.

    Our study aims to evaluate postoperative corneal equivalent stress and strain. A 3D finite element model was constructed to simulate SMILE surgery. By loading IOP, the postoperative change in corneal strain and stress were explored through diopter, IOP, and incision sizes to explore the mechanism of how SMILE influences corneal biomechanics. It was not only important for preoperative diagnosis of eye disease, but also provided effective preventive measures for postoperative complications such as keratoconus.

    We used Siemens NX (Siemens AG, Munich, Germany) to build a 3D full-eye model of the human eye. According to the anterior surface of cornea and the thickness and size of the cornea cap, the point cloud data on the subsurface of the cornea cap were obtained. Then, the point cloud data were used to subtract the thickness of lenticule calculated by the Munnerlyn equation to simulate the removal of a corneal lenticule. After the corneal lenticule was removed, the lenticule interfaces were fitted the residual cornea by modifying the upper and lower surfaces of cornea cap and the postoperative SMILE eye model was constructed. The eye model consisted of a corneal cap, cornea and sclera. And the small incision was located at the junction of cornea cap and cornea, with a size of about 3 mm.

    The IOP was loaded inner surface the eyeball, which was a fluid cavity with the eye as an enclosed inner surface. And the model was restrained in posterior scleral nerves that were surrounded by optic nerve bundle and biological tissues. In order to prevent rigid-body movement in the eye model, the force of fixation from the posterior scleral nerve was used as the the boundary conditions of model. In addition, the cornea and sclera, corneal cap and cornea were set as bond-connected, respectively. The lenticule interfaces and residual cornea were connected in frictionless way. Figure 1 was the diagram of the human eye model.

    Figure 1.  (a) The 3D whole-eye model. (b) The cross-sectional view of the model mesh. (c) The diagram of small incision model, the location of the small incision was the arc between the dots A and B, and the size was about 3 mm. The radii of curvature of the anterior and posterior corneal surfaces were 7.7 mm, 6.8 mm, respectively. The central corneal thickness was 0.5 mm, and the refractive index of the cornea was 1.376.

    Meshing a model is a critical step in finite element analysis. The mesh quality of the model affects the calculation rate and the accuracy of the experimental data. In our study, the cornea cap, cornea and sclera were divided by tetrahedral and hexahedral meshes, respectively, and the mesh mass of the eye model was about 0.7. For the models of different diopters, the number of elements and nodes were slightly different, but the changes were not significant. The numbers of corneal nodes, corneal mesh cells, scleral nodes, and scleral mesh cells were 175,765, 59,104,112,038, and 31,346, respectively.

    Grytz et al. [21] showed that the cornea and sclera were nonlinear material properties. In the establishment of the mathematical model of the eye, the Ogden strain-energy function [22] was used to represent the isotropy, hyperelasticity and incompressibility of the cornea and sclera. The specific formula can be expressed as follows:

    W=Ni=1μiβi(ˉλβi1+ˉλβi2+ˉλβi33)+Nk=11dk(J1)2k (1)

    Where βi, μi and dk are material constants, the determinant of the elastic deformation gradient is J = λ1 λ2 λ3, the principal stretch are λ1, λ2, λ3, respectively. And the deviatoric principal stretches is λi=J1  3λi.

    The shear modulus µ is defined by the following equation:

    μ=Ni=1βiμi2 (2)

    The bulk modulus K is defined by:

    K=2d1 (3)

    Due to the near incompressibility of the cornea and sclera, d1 was about 0. According to the results of previous studies and material fitting of our study, N = 2 and N = 1 were independently selected as corneal and scleral fitting parameters, respectively, where corneal fitting parameters were as follows: µ1 = 0.0034801 MPa, a1 = 104.06, µ2 = 0.0034801 MPa, and a2 = 103.94. The scleral fitting parameters were as follows: µ1 = 0.030224 MPa, and a1 = 182.73.

    SMILE creates a corneal lenticule that changes the surface curvature of the cornea to correct refractive errors. The thickness of the corneal lenticule was calculated based on the Munnerlyn equation in our study [23], the formula is defined by:

    ls(x,y)=(R2ixx2+RiyRix)2y2R2fx2y2+R2f(Oz/2)2+RixRiyR2ix(Oz/2)2 (4)

    Rix and Riy are the preoperative radii of curvature of two major meridians of the anterior corneal surface, and their values are both 7.7 mm. Rf is the postoperative radius of curvature. OZ represents the X-axis diameter of the optical zone. Rf can be expressed as:

    Rf=1000(n1)Rix(n1)+DsRix (5)

    where Ds and n are the myopic refraction in diopters and the refractive index of the cornea, respectively.

    Figure 2 showed the principles of the SMILE procedure. Firstly, the corneal lenticule was made through incisions 1 and 2; then the corneal cap was made through incision 3; and finally, 3-mm small incision 4 was created, and the lenticule was extracted via the incision. In our study, the corrected refractions were related to the ablated corneal stromal lenticule, which was given by the ablation formula. The thickness of the corneal cap was 110 μm, and the diameters of the corneal cap and optical zone were 7.6 and 6.6 mm, respectively. After inputting the IOP, the equivalent stress and strain of the anterior and posterior surfaces of the cornea were calculated using the finite element software Ansys (Ansys, Canonsburg, PA, USA). We could investigate the effects of corrected diopters, IOP, and small-incision sizes on equivalent stress and strain.

    Figure 2.  The simplified schematic diagram of the small incision lenticule extraction procedure. The lenticule cut was 1. The lenticule side cuts 2 was made. The cap interface was 3. The lenticule interfaces were created. Then the lenticule was removed through the small incision 4.

    As the cornea thins, the stress distribution becomes less uniform, and the cornea deforms by increasing corneal curvature in the thinner higher stress region to reduce local stress [24]. Based on finite element analysis, we used Ansys software to analyze the stress and strain of the cornea equivalent to that occurring after SMILE refractive surgery. Corneal deformation was observed after SMILE, and the equivalent stress could represent the stress inside the eyeball. To make the calculation more accurate, the stress on the surface of the cornea in the X, Y, and Z directions were denoted as σx, σy, and σz, respectively. Then the stress can be expressed by the following formula:

    (σxσy)2+(σxσz)2+(σyσz)2=2σs=6k2 (6)

    It can be rewritten as:

    σs=22(σxσy)2+(σxσz)2+(σyσz)2 (7)

    The equivalent stress is that the complex stress from the model of cornea equivalent to the stress in unidirectional stretching, the σe can be expressed as :

    σe=22(σxσy)2+(σxσz)2+(σyσz)2+6(τ2xy+τ2yz+τ2zx) (8)

    From the Mises yield conditions, σe = σs, we obtain:

    σe=σs=22(σxσy)2+(σxσz)2+(σyσz)2 (9)

    Then, its strain increment dεe is computed as

    dεe=(dεxdεy)2+(dεydεz)2+(dεxdεz)2+6(dγ2xy+dγ2xz+dγ2yz)2(1+v) (10)

    The equivalent strain can be obtained by the equation:

    εe=11+v12[(εxεy)2+(εyεz)2+(εxεz)2] (11)

    where v' is Poisson's ratio of the eyeball material. Poisson ratio values for the human cornea ranges from 0.42 to 0.5, which indicates that the cornea tissue is nearly incompressible. A value of 0.49 was assumed in this study.

    After setting material parameters and the mesh of the eyeball model, the IOP was loaded to simulate SMILE surgery. Considering the small incision was located in the Y-axis direction of the eye model, the equivalent stress in the Y-axis direction on the anterior and posterior corneal surfaces was calculated. In this part of the study, the influence of corrected diopters on corneal stress was explored. The optical zone diameter was 6.6 mm, and the diopter ranged from −1D to −15D, with an interval of 1 diopter. Figure 3 showed the relationship between corneal stress and myopic diopter.

    Figure 3.  (a) Equivalent stress on the anterior surface of cornea under different corrected diopters. (b) Equivalent stress on the posterior surface of cornea under different corrected diopters. The corneal cap thickness was assumed to be 110 μm and the size was 7.6 mm, IOP was 15 mmHg. The diopter ranged from 0D to −15D.

    Overall, compared with preoperative stress, the stress on the corneal anterior and posterior surfaces increased with increasing diopter, showing a regional distribution in the central region. As observed in Figure 3a, a small peak in stress occurred at −3.85 mm on the Y-axis. When the diopter was −15D, the stress reached its peak value of 0.019 MPa. The reason for this was that the corneal cap and cornea were in separated forms during construction, and stress concentration occured at their junctions. And at 3.8 mm on the Y-axis, the stress value dropped to 0.007 MPa. From Figure 3b, we concluded that the stress on the posterior surface of the cornea also increased with an increase in the corrected diopter and was concentrated in the central area. However, it showed a symmetrical distribution. When the corrected diopters were 0D and −15D, the stresses at the center of the cornea were 0.016 and 0.036 MPa, respectively. The stress value of −15D was about twice that of 0D. In addition, the stress on the posterior surface of the cornea showed no stress concentration phenomenon and changed more than that on the anterior surface.

    By loading IOP, we utilized the same method to obtain the equivalent strain on the anterior and posterior corneal surfaces. Figure 4 showed the relationship curves between strain and myopic diopter, where the diopter ranged from 0D to −15D.

    Figure 4.  (a) Equivalent strain on the anterior surface of cornea under different corrected diopters. (b) Equivalent strain on the posterior surface of cornea under different corrected diopters. The corneal cap thickness was assumed to be 110 μm and the size was 7.6 mm, IOP was 15 mmHg. The diopter ranged from 0D to −15D.

    According to Figure 4, the variation in corneal strain showed a similar pattern to that of stress. With increasing diopter, the strain increased and became more concentrated in the central area. As indicated by Figure 4a, the strain reached a small peak value at −3.85 mm with a maximum value of 0.18. At 3.8 mm, a sudden drop in strain to 0.009 occurred due to the small incision. Overall, the strain of the cornea increased with increased diopter and changed obviously at the site of incision. Figure 4b showed that the posterior surface strain of the cornea increased with a change of diopter. However, no sudden change in the strain occurred at the location of the incision. In addition, it changed more than that of the anterior surface.

    The 3D finite element eye model with a diopter of −6D was loaded different IOPs ranging from 13 to 35 mmHg with an interval of 2 mmHg. To explore the effect of IOP on anterior and posterior corneal stress, the value of stress variation on the corneal surfaces under different IOPs was calculated. The results were showed in Figure 5.

    Figure 5.  (a) The equivalent stress on the anterior surface of the cornea under different IOP. (b) The equivalent stress on the posterior surface of the cornea under different IOP. IOP ranged from 13 to 35 mmHg with intervals of 2 mmHg. The diopter was −6D.

    In general, stress on the anterior and posterior surface of the cornea increased linearly with increasing IOP, whereas stress on the anterior surface of the cornea suddenly dropped to 0.001 MPa at the time of incision. According to Figure 5a, when the IOP was 13, 15 and 17 mmHg, the maximum stress were about 0.013, 0.015 and 0.017 MPa, respectively. Until the IOP was 35 mmHg, the maximum stress was 0.035 MPa. The result showed that the stress was increased by about 0.002 MPa at an interval of 2 mmHg. As observed in Figure 5b, the stress on the posterior surface of the cornea also increased linearly. When the IOP was 13, 15 and 17 mmHg, the maximum stress were about 0.017, 0.021 and 0.024 MPa, respectively. In conclusion, the stress on the posterior of the cornea increased by about 0.003 MPa with every 2 mmHg increase of IOP.

    Different IOPs were applied to the eye model with a diopter of −6D and the equivalent strain on the corneal surface obtained. The IOP ranged from 13 to 35 mmHg with intervals of 2 mmHg. The results were as follows.

    Figure 6 indicated that the strain on the anterior and posterior surfaces of the cornea increased nonlinearly with increasing IOP. Different from the effect of IOP on corneal stress, the strain on the anterior surface of the cornea decreased slowly to 0.007 at the incision site. Figure 6a showed that the strain difference at the center decreased under the same IOP interval. When the IOP increased from 13 to 15 mmHg, the strain at the center increased by about 0.0025. The strain at the center increased by about 0.0022 from 15 to 17 mmHg. Between 17 and 19 mmHg, the strain at the center increased by about 0.0014. According to Figure 6b, when the IOP increased from 13 to 15 mmHg, the posterior cornea surface of strain at the center increased by about 0.0028. The strain at the center increased by about 0.0026 between 15 to 17 mmHg. In summary, the change of the strain on the posterior surface of the cornea was not linear with increasing IOP.

    Figure 6.  (a) The equivalent strain on the anterior surface of the cornea under different IOP. (b) The equivalent strain on the posterior surface of the cornea under different IOP. IOP ranged from 13 to 35 mmHg with intervals of 2 mmHg. The diopter was −6D.

    The biomechanical response of the cornea is affected by eye diseases and surgical stability [25]. The flapless and small incision of SMILE refractive surgery has improved surgical safety [26]. Clinical studies have shown the biomechanical of corneal changed after surgery. Wu et al. [27] showed that corneal hysteresis and the corneal resistance factor both decreased, obviously changing the biomechanics of the corneal. Balidis et al. [28] showed that corneal biomechanics changed after SMILE. In addition, many studies showed that the stress and strain increased after surgery. Shih et al. [29] found that the central area of the cornea was the most deformed after SMILE. Franus et al. [8] also suggested that the postoperative cornea was more deformed than the normal cornea. The finite element analysis has been used to research corneal biomechanics after refractive surgery. And the corneal biomechanics were also explored by modeling research. Shih et al. [29] used finite element method to simulate the changes of corneal shape and stress distribution after refractive surgical, which showed the changes of biomechanical cornea. Cordero showed that using a corneal model to analyze biomechanics were useful tools could quickly diagnose eye diseases [30]. NT Mohammad et al. [31] also used finite element model to assessing the mechanical behavior of the cornea. Therefore, previous studies have confirmed that corneal biomechanics changed after corneal refractive surgery. Our results suggested that postoperative corneal biomechanics might also change, the stress and strain increased. And the increases were concentrated in the central region after SMILE surgery. By constructing a 3D finite element model of the human eye and using computer algorithm to analyze the corneal biomechanics, our study had a certain theoretical significance for refractive surgery.

    The sizes of incisions ranged from 2 to 5 mm with intervals of 0.5 mm. After applying 15 mmHg IOP, the stress and strain on the anterior surface of the cornea for different incision sizes were calculated. The results were showed in Figure 7.

    Figure 7.  (a) Relationship between the sizes of incision and equivalent stress. (b) Relationship between the sizes of incision and equivalent strain. The incision size ranged from 2 to 5 mm, with an interval of 0.5 mm, and the diopter was −6D. IOP was 15 mmHg.

    It was concluded from Figure 7 that the stress and strain on the anterior surface of the cornea increased with increasing incision size, and the changes were obviously at the incision site. Compared with diopter and IOP, incision size had little influence on corneal stress and strain. As observed in Figure 7a, the maximum stress in the central region was 0.016 MPa, and the change in value at the incision site was 0.001 MPa. According to Figure 7b, how strain and stress varies was consistent, and strain increased with increasing incision size. Overall, the incision size had only a slight effect on stress and strain in SMILE surgery, which might be associated with the small changes of corneal tension.

    Our results only considered stress in the Y-axis direction, but corneal stress in the X-axis direction should also be considered. Therefore, we calculated the X-axis direction stress changes on the anterior and posterior surfaces of the cornea under different diopters, ranging from 0D to −15D. The results were shown in Figure 8.

    Figure 8.  (a) The X-axis direction equivalent stress on the anterior surface of cornea under different corrected diopters. (b) The X-axis direction equivalent stress on the posterior surface of cornea under different corrected diopters. IOP was 15 mmHg.

    Overall, stress also increased on the X-axis direction with increasing diopter. However, it was clear from Figure 8a that the stress on the anterior surface of the cornea was symmetrical, and no descent was observed at the incision location. The stress peaked at −3.8 and 3.8 mm due to the construction of the corneal cap and corneal dissociation. Figure 8b showed that the stress on the posterior surface of the cornea changed in a similar way to that on the Y-axis direction, and the stress on the posterior surface was greater than that on the anterior surface.

    Our study assessed the biomechanical responses of the cornea after SMILE surgery, which could provide references for the SMILE refractive surgery. But there were some imperfections that need to be improved. Firstly, SMILE surgery had a greater impact on the cornea, so our model was composed of the corneal cap, the cornea, and the sclera. However, in order to make the study results more accurate, the impact of other components should also be considered, which needed to be improved in future jobs. Secondly, the material properties of the cornea and sclera had also an important effect on the experimental results. Because of the nonlinear and hyperelastic properties of the cornea and sclera, the Ogden function model was selected. But there were individual differences in the human eye in clinical practice. Finally, we only considered the effects of the diopter, IOP, and incision size on stress and strain. In practice, the factors influencing stress and strain may included treatment decentration, material parameters, corneal ablation depth, and others. How all these factors influence corneal stress and strain required further exploration.

    The biomechanics of ophthalmological procedures can be accurately analyzed by finite element analysis. Our next target is dynamic simulation of the clinical SMILE surgical procedure, so that our research findings will be closer to those of clinical surgery. Moreover, there are many factors that influence the outcome of surgery. We will study the postoperative biomechanics of the human eye from perspectives such as eye model, treatment decentration, and material parameters. Therefore, the finite element analysis used to study postoperative corneal biomechanics can provide a theoretical basis for making surgical plans and diagnosing postoperative ocular diseases after refractive surgery.

    A 3D human eye model was used to simulate the SMILE surgery procedure based on the finite element analysis method, and we studied postoperative corneal stress and strain. Our results showed that with increasing corrected diopter, the equivalent stress and strain on the anterior and posterior surfaces of the cornea increased and were concentrated in the central area. However, compared with the anterior surface of the cornea, the changes in stress and strain on the posterior surface were significantly larger. With increasing IOP, corneal stress increased approximately linearly, whereas the strain increased nonlinearly. In addition, we found that the sizes and directions of incisions had a slight influence on stress and strain. In summary, our study indicated SMILE refractive surgery increased stress and strain, having certain influence on corneal biomechanics. And the result could guide the design of clinical operation and the avoidance of postoperation complications.

    Supported by Natural National Science Foundation of China (NSFC) (61465010 & 81873684); National Key Research and Development Program of China (2018YFE0115700); Jiangxi Nature Science Foundation (20192BAB207035).

    The authors thank the participants of this multicentre study for their valuable contributions.

    The authors assume full responsibility for the analyses and interpretation of these data.

    All authors declare no conflicts of interest in this paper.



    [1] J. Wang, X. Xun, D. O. Ophthalmology, S. G. Hospital, Environmental risk factors for myopia and advances in its control and prevention, Shanghai Med. Pharm. J., 38 (2017), 3-7.
    [2] M. Kohlhaas, Corneal sensation after cataract and refractive surgery, J. Cataract Refract. Surg., 24 (1998), 1399-1409. doi: 10.1016/S0886-3350(98)80237-X
    [3] D. Z. Reinstein, T. J. Archer, M. Gobbe, The Key characteristics of corneal refractive surgery: biomechanics, spherical aberration, and corneal sensitivity after SMILE, in Small Incision Lenticule Extraction (SMILE), Springer, Cham, (2015), 123-142.
    [4] A. Gyldenkerne, A. Ivarsen, J. Hjortdal, Comparison of corneal shape changes and aberrations induced By FS-LASIK and SMILE for myopia, J. Refractive Surg., 31 (2015), 223-229. doi: 10.3928/1081597X-20150303-01
    [5] P. Wei, G. P. Cheng, J. Zhang, A.L. Ng, Y. Wang, Changes in corneal volume at different areas and its correlation with corneal biomechanics after SMILE and FS-LASIK surgery, J. Ophthalmol., 2020 (2020), 1-7.
    [6] G. Scarcelli, S. H. Yun, Brillouin scanning microscopy in keratoconus, in Keratoconus: Recent Advances in Diagnosis and Treatment, Springer International Publishing, (2017), 167-173.
    [7] N. Cartwright, J. R. Tyrer, P. D. Jaycock, J. Marshall, Effects of variation in depth and side cut angulations in LASIK and thin-flap LASIK using a femtosecond laser: a biomechanical study, J. Refractive Surg., 28 (2012), 419-425. doi: 10.3928/1081597X-20120518-07
    [8] D. V. Franus, Change in the stress-strain state of the cornea after refractive surgery, in 2015 International Conference on Mechanics-Seventh Polyakhov's Reading, IEEE, (2015), 1-4.
    [9] J. C. He, J. Gwiazda, F. Thorn, R. Held, Wave-front aberrations in the anterior corneal surface and the whole eye, J. Opt. Soc. Am. A, 20 (2003), 1155-1163. doi: 10.1364/JOSAA.20.001155
    [10] A. S. Roy, W. J. Dupps, Effects of altered corneal stiffness on native and postoperative LASIK corneal biomechanical behavior: a whole-eye finite element analysis, J. Refractive Surge., 25 (2009), 875-887. doi: 10.3928/1081597X-20090917-09
    [11] K. J. Bathe, D. Chapelle, Computational Fluid And Solid Mechanics: the Finite Element Analysis of Shells- Fundamentals, Springer, 2003.
    [12] M. Á. Ariza-Gracia, J. F. Zurita, D. P. Pinero, J. F. Rodriguez, Coupled biomechanical response of the cornea assessed by non-contact tonometry, a simulation study, PLOS ONE, 10 (2015), e0121486.
    [13] M. Kaliske, A formulation of elasticity and viscoelasticity for fibre reinforced material at small and finite strains, Comput. Methods Appl. Mech. Eng., 185 (2000), 225-243. doi: 10.1016/S0045-7825(99)00261-3
    [14] K. Anderson, A. El-Sheikh, T. Newson, Application of structural analysis to the mechanical behaviour of the cornea, J. R. Soc. Interface, 1 (2004), 1742-5662.
    [15] P. M. Pinsky, D. V. Datye, A microstructurally-based finite element model of the incised human cornea, J. Biomech., 24 (1991), 907-922. doi: 10.1016/0021-9290(91)90169-N
    [16] R. Shah, S. Shah, S. Sengupta, Results of small incision lenticule extraction: all-in-one femtosecond laser refractive surgery, J. Cataract Refractive Surg., 37 (2011), 127-137. doi: 10.1016/j.jcrs.2010.07.033
    [17] J. B. Randleman, B. Russell, M. A. Ward, K. P. Thompson, R. D. Stulting, Risk factors and prognosis for corneal ectasia after LASIK, Ophthalmology, 110 (2003), 267-275. doi: 10.1016/S0161-6420(02)01727-X
    [18] K. A. John, Comparison of corneal biomechanics after myopic small-incision lenticule extraction compared to LASIK: an ex vivo study, Clin. Ophthalmol., 12 (2018), 237-245. doi: 10.2147/OPTH.S153509
    [19] I. M. Osman, H. A. Helaly, M. Abdalla, M. A. Shousha, Corneal biomechanical changes in eyes with small incision lenticule extraction and laser assisted in situ keratomileusis, BMC Ophthalmol., 16 (2016), 123. doi: 10.1186/s12886-016-0304-3
    [20] T. F. Peinado, D. P. Piñ ero, I. A. López, J. L. Alio, Correlation of both corneal surfaces in corneal ectasia after myopic LASIK, Optom. Vision Sci., 88 (2011), E539-E542.
    [21] R. Grytz, K. Krishnan, R. Whitley, V. Libertiaux, J. C. Downs, A mesh-free approach to incorporate complex anisotropic and heterogeneous material properties into eye-specific finite element models, Comput. Methods Applied Mech. Eng., 358 (2020), 112654. doi: 10.1016/j.cma.2019.112654
    [22] J. G. Yu, F. J. Bao, Y. F. Feng, C. Whitford, A. Elsheikh, Assessment of corneal biomechanical behavior under posterior and anterior pressure, J. Refractive Surg., 29 (2013), 64-70. doi: 10.3928/1081597X-20121228-05
    [23] L. Fang, The influence of the aspheric profiles for transition zone on optical performance of human eye after conventional ablation, J. Eur. Opt. Soc. Rapid Publ., 9 (1990).
    [24] C. J. Roberts, W. J. Dupps, Biomechanics of corneal ectasia and biomechanical treatments, J. Cataract Refractive Surg., 40 (2014), 991-998. doi: 10.1016/j.jcrs.2014.04.013
    [25] M. Girard, W. J. Dupps, M. Baskaran, G. Scarcelli, S. H. Yun, H. A. Quigley, et al., Translating ocular biomechanics into clinical practice: current state and future prospects, Curr. Eye Res., 40 (2015), 1-18. doi: 10.3109/02713683.2014.914543
    [26] S. Ganesh, S. Brar, R. R. Arra, Refractive lenticule extraction small incision lenticule extraction: a new refractive surgery paradigm, Indian J. Ophthalmol., 66 (2018), 10-19. doi: 10.4103/ijo.IJO_761_17
    [27] W. Di, Y. Wang, L. Zhang, S. Wei, T. Xin, Corneal biomechanical effects: Small-incision lenticule extraction versus femtosecond laser-assisted laser in situ keratomileusis, J. Cataract Refractive Surg., 40 (2014), 954-962. doi: 10.1016/j.jcrs.2013.07.056
    [28] M. Balidis, Biomechanical profile of refractive surgery procedures, Acta Ophthalmol., 97 (2019).
    [29] P. J. Shih, I. J. Wang, W. F. Cai, J. Y. Yen, Biomechanical simulation of stress concentration and intraocular pressure in corneas subjected to myopic refractive surgical procedures, Sci. Rep., 7 (2017), 13906. doi: 10.1038/s41598-017-14293-0
    [30] M. I. Cordero-Mendieta, E. Pinos-Vélez, R. Coronel-Berrezueta, Study of corneal biomechanics and modeling of Young's module in healthy and pathological corneas, in International Conference on Applied Human Factors and Ergonomics, Springer, Cham, (2020), 99-105.
    [31] N. T. Mohammad, F. Craig, G. Dipika, Finite element modelling of cornea mechanics: a review, Arq. Bras. Oftalmol., 77 (2014), 60-65.
  • This article has been cited by:

    1. Francisco J. Perdomo-Argüello, Estelina Ortega-Gómez, Purificación Galindo-Villardón, Víctor Leiva, Purificación Vicente-Galindo, STATIS multivariate three-way method for evaluating quality of life after corneal surgery: Methodology and case study in Costa Rica, 2023, 20, 1551-0018, 6110, 10.3934/mbe.2023264
    2. Gongpu Lan, Qun Shi, Yicheng Wang, Guoqin Ma, Jing Cai, Jinping Feng, Yanping Huang, Boyu Gu, Lin An, Jingjiang Xu, Jia Qin, Michael D. Twa, Spatial Assessment of Heterogeneous Tissue Natural Frequency Using Micro-Force Optical Coherence Elastography, 2022, 10, 2296-4185, 10.3389/fbioe.2022.851094
    3. Yan Wang, Lixin Xie, Ke Yao, Walter Sekundo, Jorge L. Alió, Jod S. Mehta, Sanjay Goel, Ahmed Elmassry, Julie Schallhorn, Tatiana Shilova, Huazheng Cao, Lulu Xu, Xuan Chen, Fengju Zhang, Ji Bai, Wenfang Zhang, Qin Liu, Xingtao Zhou, Yueguo Chen, Zheng Wang, Vishal Jhanji, Kehu Yang, Evidence-based Guidelines for Keratorefractive Lenticule Extraction Surgery, 2024, 01616420, 10.1016/j.ophtha.2024.11.016
    4. Yifeng Li, Zhuoran Yang, Ziming Yan, Huibin Shi, Zhanli Liu, Kaijie Wang, Revealing and predicting the long-term biomechanical response of orthokeratology by developing a patient-specific computational model, 2025, 68, 1674-7348, 10.1007/s11433-024-2487-6
    5. 郭济西 Guo Jixi, 方利华 Fang Lihua, 石刚 Shi Gang, 于瑞轩 Yu Ruixuan, 张毅 Zhang Yi, 汪新超 Wang Xinchao, 张旭阳 Zhang Xuyang, 光学相干弹性成像评估SMILE术后角膜弹性的有限元分析, 2025, 52, 0258-7025, 0307207, 10.3788/CJL241231
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3378) PDF downloads(177) Cited by(5)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog