%0 Journal Article
%A Xu, Mengchen
%A Paul, Funkenbusch D.
%A Lerner, Amy L.
%A Yoon, Geunyoung
%T Estimating the material properties of collagen fibril using optical behavior of the human cornea
%B Investigative Ophthalmology & Visual Science
%D 2017
%J Investigative Ophthalmology & Visual Science
%V 58
%N 8
%P 4332-4332
%@ 1552-5783
%X To improve the accuracy of estimating collagen fibril properties by using optical behavior of the cornea with an inverse analysis method. A three-dimensional anisotropic corneal model with regionally varying fibril distribution (Pandolfi, 2008) was generated in finite element software (ABAQUS 6.14). The model used two known average material parameters, depth-dependent matrix stiffness C10 (Sloan, 2014) and fibril dispersion parameter κ (Pandolfi, 2008), and two unknown parameters for the model, fibril stiffness k1 and fibril nonlinearity k2 were estimated with an inverse algorithm developed based on design of experiment method. This iterative method estimates k1 and k2 by minimizing the average difference between the predicted and the theoretical corneal responses simulated under multiple intraocular pressure (IOP) levels ranging from 5-30mmHg. To test the accuracy, three different combinations of k1 and k2 were assumed: k1 = 30kPa, k2 = 800 (case 1); k1 = 80kPa, k2 = 800 (case 2); k1 = 80kPa, k2 = 1200 (case 3). We compared estimated k1 and k2 based on apical displacement and spherical aberration (SA). The iterative algorithm was considered converged with a threshold of 0.01mm for apical displacement and 0.1μm for SA. Accuracy of the estimation was judged by prediction error of apical displacement and SA, as well as refractive power as an independent corneal response. Estimation of k1 and k2 using apical displacement led to non-unique solutions for all three cases. Different combinations of k1 and k2 predicted similar apical displacement and SA within the thresholds. However, average estimation error of k1 and k2 relative to their assumed magnitudes for the three cases was 28.8% and 18.8%, which led to an average error of 0.25D in refractive power compared to the theoretical data. In contrast, the algorithm based on SA resulted in a unique solution of k1 and k2 and the error was reduced to 8.5% and 2.8%. In addition, the resulting corneal responses for the three cases all matched the theoretical data, with an average difference of 0.006mm in apical displacement, 0.006μm in SA and 0.13D in refractive power. For this model with regionally varying fibril distribution, an inverse analysis based on optical responses of the cornea at multiple IOPs improved the accuracy in estimating corneal collagen fibril properties compared to the result based on apical displacement. This is an abstract that was submitted for the 2017 ARVO Annual Meeting, held in Baltimore, MD, May 7-11, 2017.
%[ 8/15/2020