Verification of a virtual fields method to extract the mechanical properties of human optic nerve head tissues in vivo

We aimed to verify a custom virtual fields method (VFM) to estimate the patient-specific biomechanical properties of human optic nerve head (ONH) tissues, given their full-field deformations induced by intraocular pressure (IOP). To verify the accuracy of VFM, we first generated 'artificial' ONH displacements from predetermined (known) ONH tissue biomechanical properties using finite element analysis. Using such deformations, if we are able to match back the known biomechanical properties, it would indicate that our VFM technique is accurate. The peripapillary sclera was assumed anisotropic hyperelastic, while all other ONH tissues were considered isotropic. The simulated ONH displacements were fed into the VFM algorithm to extract back the biomechanical properties. The robustness of VFM was also tested against rigid body motions and noise added to the simulated displacements. Then, the computational speed of VFM was compared to that of a gold-standard stiffness measurement method (inverse finite element method or IFEM). Finally, as proof of principle, VFM was applied to IOP-induced ONH deformation data (obtained from one subject's eye imaged with OCT), and the biomechanical properties of the prelamina and lamina cribrosa (LC) were extracted. From given ONH displacements, VFM successfully matched back the biomechanical properties of ONH tissues with high accuracy and efficiency. For all parameters, the percentage errors were less than 0.05%. Our method was insensitive to rigid body motions and was also able to recover the material parameters in the presence of noise. VFM was also found 125 times faster than the gold-standard IFEM. Finally, the estimated shear modulus for the prelamina and the LC of the studied subject's eye were 33.7 and 63.5 kPa, respectively. VFM may be capable of measuring the biomechanical properties of ONH tissues with high speed and accuracy. It has potential in identifying patient-specific ONH biomechanical properties in the clinic if combined with optical coherence tomography.