Revealing and predicting the long-term biomechanical response of orthokeratology by developing a patient-specific computational model

作     者：Yifeng Li Zhuoran Yang Ziming Yan Huibin Shi Zhanli Liu Kaijie Wang 

作者机构：Department of Engineering Mechanics School of Aerospace Engineering Tsinghua University Beijing Tongren Eye Center Beijing Tongren Hospital Capital Medical University Beijing Ophthalmology & Visual Sciences Key Lab 

出 版 物：《Science China(Physics,Mechanics & Astronomy)》 (中国科学:物理学 力学 天文学(英文版))

年 卷 期：2025年第Mechanics & Astronomy) .卷第1期

页      面：156-168页

核心收录：

学科分类：0831[工学-生物医学工程（可授工学、理学、医学学位）] 1002[医学-临床医学] 08[工学] 100212[医学-眼科学] 0836[工学-生物工程] 0801[工学-力学（可授工学、理学学位）] 10[医学] 

基　　金：supported by the National Natural Science Foundation of China (Grant Nos. 82371038, and 11921002) the R&D Program of Beijing Municipal Education Commission (Grant No. KZ20231002503) the Chinese Institutes for Medical Research,Beijing (Grant No. CX24PY17) Beijing Municipal Administration of Hospitals Incubating Program (Grant No. PX2024011) 

主　　题：orthokeratology computational biomechanical model long-term viscoelasticity myopia correction 

摘      要：Orthokeratology(OK) is widely used for effective myopia correction and control. However, the incomplete understanding of its biomechanical mechanisms makes OK lens fitting rely heavily on clinician judgment, complicating accurate predictions of treatment outcomes. In this paper, we performed clinical experiments and numerical analysis to study corneal deformation modes and long-term changes in central corneal thickness. Clinical experiments were conducted on 194 Chinese myopic patients under OK treatment for 3 months. Based on the experimental data, a patient-specific computational biomechanical model for OK was established and validated. Specifically, the anisotropic mechanical properties of the cornea were incorporated into the model to describe the significant difference between its shear modulus(29.5 kPa) and tensile modulus(768.4 kPa). Additionally, a viscohyperelastic material model with a prolonged corneal relaxation time of 5.6 h was developed to capture the long-term deformation response. The results show that corneal thickness reduction in OK is primarily due to out-of-plane shear deformation,influenced by the cornea s low shear resistance. Modeling the extended corneal relaxation time is crucial for predicting long-term biomechanical responses. The computational model effectively captures long-term changes in central corneal thickness, potentially improving OK lens fitting accuracy.