Technical Brief

Appropriate Objective Functions for Quantifying Iris Mechanical Properties Using Inverse Finite Element Modeling

[+] Author and Article Information
Anup D. Pant

Department of Biomedical Engineering,
The University of Akron,
Akron, OH 44325
e-mail: adp63@zips.uakron.edu

Syril K. Dorairaj

Department of Ophthalmology,
Mayo Clinic,
Jacksonville, FL 32224
e-mail: dorairaj.syril@mayo.edu

Rouzbeh Amini

Department of Biomedical Engineering,
The University of Akron,
Akron, OH 44325
e-mail: ramini@uakron.edu

1Corresponding author.

Manuscript received July 7, 2017; final manuscript received February 27, 2018; published online April 30, 2018. Assoc. Editor: Thao (Vicky) Nguyen.

J Biomech Eng 140(7), 074502 (Apr 30, 2018) (6 pages) Paper No: BIO-17-1298; doi: 10.1115/1.4039679 History: Received July 07, 2017; Revised February 27, 2018

Quantifying the mechanical properties of the iris is important, as it provides insight into the pathophysiology of glaucoma. Recent ex vivo studies have shown that the mechanical properties of the iris are different in glaucomatous eyes as compared to normal ones. Notwithstanding the importance of the ex vivo studies, such measurements are severely limited for diagnosis and preclude development of treatment strategies. With the advent of detailed imaging modalities, it is possible to determine the in vivo mechanical properties using inverse finite element (FE) modeling. An inverse modeling approach requires an appropriate objective function for reliable estimation of parameters. In the case of the iris, numerous measurements such as iris chord length (CL) and iris concavity (CV) are made routinely in clinical practice. In this study, we have evaluated five different objective functions chosen based on the iris biometrics (in the presence and absence of clinical measurement errors) to determine the appropriate criterion for inverse modeling. Our results showed that in the absence of experimental measurement error, a combination of iris CL and CV can be used as the objective function. However, with the addition of measurement errors, the objective functions that employ a large number of local displacement values provide more reliable outcomes.

Copyright © 2018 by ASME
Your Session has timed out. Please sign back in to continue.


Quigley, H. A. , and Broman, A. T. , 2006, “ The Number of People With Glaucoma Worldwide in 2010 and 2020,” Br. J. Ophthalmol., 90(3), pp. 262–267. [CrossRef] [PubMed]
Quigley, H. A. , 1996, “ Number of People With Glaucoma Worldwide,” Br. J. Ophthalmol., 80(5), pp. 389–393. [CrossRef] [PubMed]
Sigal, I. A. , Grimm, J. L. , Jan, N.-J. , Reid, K. , Minckler, D. S. , and Brown, D. J. , 2014, “ Eye-Specific IOP-Induced Displacements and Deformations of Human Lamina Cribrosa direct Measurement of Lamina Cribrosa Biomechanics,” Invest. Ophthalmol. Visual Sci., 55(1), pp. 1–15. [CrossRef]
AGIS Investigators, 2000, “ The Advanced Glaucoma Intervention Study (AGIS)—7: The Relationship Between Control of Intraocular Pressure and Visual Field Deterioration,” Am. J. Ophthalmol., 130, pp. 429–440. [CrossRef] [PubMed]
Gordon, M. O. , Beiser, J. A. , Brandt, J. D. , Heuer, D. K. , Higginbotham, E. J. , Johnson, C. A. , Keltner, J. L. , Miller, J. P. , Parrish, R. K. , Wilson, M. R. , and Kass, M. A. , 2002, “ The Ocular Hypertension Treatment Study: Baseline Factors That Predict the Onset of Primary Open-Angle Glaucoma,” Arch. Ophthalmol., 120(6), pp. 714–720. [CrossRef] [PubMed]
Herschler, J. , 1977, “ Trabecular Damage Due to Blunt Anterior Segment Injury and Its Relationship to Traumatic Glaucoma,” Trans. Sect. Ophthalmol. Am. Acad. Ophthalmol. Otolaryngol., 83(2), pp. 239–248. [PubMed]
Gsellman, L. , and Amini, R. , 2016, “ Patients With Intravitreal Gas Bubbles at Risk of High Intraocular Pressure Without Exceeding Elevation of Surgery: Theoretical Analysis,” Invest. Ophthalmol. Visual Sci., 57(7), pp. 3340–3347. [CrossRef]
Amini, R. , Barocas, V. H. , Kavehpour, H. P. , and Hubschman, J. P. , 2011, “ Computational Simulation of Altitude Change-Induced Intraocular Pressure Alteration in Patients With Intravitreal Gas Bubbles,” Retina, 31(8), pp. 1656–1663. [CrossRef] [PubMed]
Jouzdani, S. , Amini, R. , and Barocas, V. H. , 2013, “ Contribution of Different Anatomical and Physiologic Factors to Iris Contour and Anterior Chamber Angle Changes During Pupil Dilation: Theoretical Analysis,” Invest. Ophthalmol. Visual Sci., 54(4), pp. 2977–2984. [CrossRef]
Barocas, V. H. , and Taravella, M. J. , 2001, “ Modeling Passive Mechanical Interaction Between Aqueous Humor and Iris,” ASME J. Biomech. Eng., 123(6), pp. 540–547. [CrossRef]
Narayanaswamy, A. K. , Hoon, N. M. , Lim, C. T. , and Aung, T. , 2015, “ Young's Modulus Determination of Normal and Glaucomatous Human Irides,” Invest. Ophthalmol. Visual Sci., 56(7), pp. 6139–6139.
Quigley, H. A. , Silver, D. M. , Friedman, D. S. , He, M. , Plyler, R. J. , Eberhart, C. G. , Jampel, H. D. , and Ramulu, P. , 2009, “ Iris Cross-Sectional Area Decreases With Pupil Dilation and Its Dynamic Behavior is a Risk Factor in Angle Closure,” J. Glaucoma, 18(3), pp. 173–179. [CrossRef] [PubMed]
Aptel, F. , and Denis, P. , 2010, “ Optical Coherence Tomography Quantitative Analysis of Iris Volume Changes After Pharmacologic Mydriasis,” Ophthalmology, 117(1), pp. 3–10. [CrossRef] [PubMed]
Quigley, H. A. , 2010, “ The Iris is a Sponge: A Cause of Angle Closure,” Ophthalmology, 117, pp. 1–2. [CrossRef] [PubMed]
Brouwer, I. , Ustin, J. , Bentiey, L. , Dhruv, A. , and Tendick, F. , 2001, “ Measuring In Vivo Animal Soft Tissue Properties for Haptic Modeling in Surgical Simulation,” Medicine Meets Virtual Reality, Vol. 81, IOS Press, Amsterdam, The Netherlands, pp. 69–74.
Kagemann, L. , Wang, B. , Wollstein, G. , Ishikawa, H. , Nevins, J. E. , Nadler, Z. , Sigal, I. A. , Bilonick, R. A. , and Schuman, J. S. , 2014, “ IOP Elevation Reduces Schlemm's Canal Cross-Sectional Area,” Invest. Ophthalmol. Visual Sci., 55(3), pp. 1805–1809. [CrossRef]
Pant, A. D. , Kagemann, L. , Schuman, J. S. , Sigal, I. A. , and Amini, R. , 2017, “ An Imaged-Based Inverse Finite Element Method to Determine In-Vivo Mechanical Properties of the Human Trabecular Meshwork,” J. Model. Ophthalmol., 1(3), pp. 100–111. [PubMed]
Lee, C.-H. , Amini, R. , Gorman, R. C. , Gorman, J. H. , and Sacks, M. S. , 2014, “ An Inverse Modeling Approach for Stress Estimation in Mitral Valve Anterior Leaflet Valvuloplasty for In-Vivo Valvular Biomaterial Assessment,” J. Biomech., 47(9), pp. 2055–2063. [CrossRef] [PubMed]
Woo, E. K. , Pavlin, C. J. , Slomovic, A. , Taback, N. , and Buys, Y. M. , 1999, “ Ultrasound Biomicroscopic Quantitative Analysis of Light-Dark Changes Associated With Pupillary Block,” Am. J. Ophthalmol., 127(1), pp. 43–47. [CrossRef] [PubMed]
Amini, R. , Whitcomb, J. E. , Prata, T. S. , Dorairaj, S. , Liebmann, J. M. , Ritch, R. , and Barocas, V. H. , 2010, “ Quantification of Iris Concavity,” J. Ophthalmic Vision Res., 5(3), pp. 211–212.
Salim, S. , and Dorairaj, S. , 2013, “ Anterior Segment Imaging in Glaucoma,” Seminars in Ophthalmology, Vol. 28, Taylor & Francis, Oxford, UK, pp. 113–125. [CrossRef]
Nagel, T. M. , Hadi, M. F. , Claeson, A. A. , Nuckley, D. J. , and Barocas, V. H. , 2014, “ Combining Displacement Field and Grip Force Information to Determine Mechanical Properties of Planar Tissue With Complicated Geometry,” ASME J. Biomech. Eng., 136(11), p. 114501. [CrossRef]
Amini, R. , Whitcomb, J. E. , Al-Qaisi, M. K. , Akkin, T. , Jouzdani, S. , Dorairaj, S. , Prata, T. , Illitchev, E. , Liebmann, J. M. , Ritch, R. , and Barocas, V. H. , 2012, “ The Posterior Location of the Dilator Muscle Induces Anterior Iris Bowing During Dilation, Even in the Absence of Pupillary Block,” Invest. Ophthalmol. Visual Sci., 53(3), pp. 1188–1194. [CrossRef]
Bonet, J. , and Wood, R. D. , 1997, Nonlinear Continuum Mechanics for Finite Element Analysis, Cambridge University Press, Cambridge, UK.
Pence, T. J. , and Gou, K. , 2015, “ On Compressible Versions of the Incompressible Neo-Hookean Material,” Math. Mech. Solids, 20(2), pp. 157–182. [CrossRef]
Beatty, M. F. , and Stalnaker, D. , 1986, “ The Poisson Function of Finite Elasticity,” ASME J. Appl. Mech., 53(4), pp. 807–813. [CrossRef]
Scott, N. , 2007, “ The Incremental Bulk Modulus, Young's Modulus and Poisson's Ratio in Nonlinear Isotropic Elasticity: Physically Reasonable Response,” Math. Mech. Solids, 12(5), pp. 526–542. [CrossRef]
Storn, R. , and Price, K. , 1997, “ Differential Evolution—A Simple and Efficient Heuristic for Global Optimization Over Continuous Spaces,” J. Global Optim., 11(4), pp. 341–359. [CrossRef]
Heys, J. , and Barocas, V. , 1999, “ Mechanical Characterization of the Bovine Iris,” J. Biomech., 32(9), pp. 999–1003. [CrossRef] [PubMed]
OSC, 1987, “ Ohio Supercomputer Center,” Ohio Supercomputer Center, Columbus, OH, accessed Mar. 28, 2018 , http://osc.edu/ark:/19495/f5s1ph73
Zielinski, K. , and Laur, R. , 2008, “ Stopping Criteria for Differential Evolution in Constrained Single-Objective Optimization,” Advances in Differential Evolution, Springer, Berlin, pp. 111–138. [CrossRef]
Amundsen, B. H. , Helle-Valle, T. , Edvardsen, T. , Torp, H. , Crosby, J. , Lyseggen, E. , Støylen, A. , Ihlen, H. , Lima, J. A. , Smiseth, O. A. , and Slørdahl, S. A. , 2006, “ Noninvasive Myocardial Strain Measurement by Speckle Tracking Echocardiography: Validation Against Sonomicrometry and Tagged Magnetic Resonance Imaging,” J. Am. Coll. Cardiol., 47(4), pp. 789–793. [CrossRef] [PubMed]
Zhang, D. , and Arola, D. D. , 2004, “ Applications of Digital Image Correlation to Biological Tissues,” J. Biomed. Opt., 9(4), pp. 691–699. [CrossRef] [PubMed]
Ford, M. R. , Roy, A. S. , Rollins, A. M. , and Dupps, W. J. , 2014, “ Serial Biomechanical Comparison of Edematous, Normal, and Collagen Crosslinked Human Donor Corneas Using Optical Coherence Elastography,” J. Cataract Refractive Surg., 40(6), pp. 1041–1047. [CrossRef]
Pavlatos, E. , Perez, B. C. , Morris, H. J. , Chen, H. , Palko, J. R. , Pan, X. , Weber, P. A. , Hart, R. T. , and Liu, J. , 2016, “ Three-Dimensional Strains in Human Posterior Sclera Using Ultrasound Speckle Tracking,” ASME J. Biomech. Eng., 138(2), p. 021015. [CrossRef]
Midgett, D. E. , Pease, M. E. , Jefferys, J. L. , Patel, M. , Franck, C. , Quigley, H. A. , and Nguyen, T. D. , 2017, “ The Pressure-Induced Deformation Response of the Human Lamina Cribrosa: Analysis of Regional Variations,” Acta Biomater., 53, pp. 123–139. [CrossRef] [PubMed]
Pan, B. , Qian, K. , Xie, H. , and Asundi, A. , 2009, “ Two-Dimensional Digital Image Correlation for In-Plane Displacement and Strain Measurement: A Review,” Meas. Sci. Technol., 20(6), p. 062001. [CrossRef]
Amini, R. , Jouzdani, S. , and Barocas, V. H. , 2012, “ Increased Iris–Lens Contact Following Spontaneous Blinking: Mathematical Modeling,” J. Biomech., 45(13), pp. 2293–2296. [CrossRef] [PubMed]
Huang, E. C. , and Barocas, V. H. , 2004, “ Active Iris Mechanics and Pupillary Block: Steady-State Analysis and Comparison With Anatomical Risk Factors,” Ann. Biomed. Eng., 32(9), pp. 1276–1285. [CrossRef] [PubMed]
Huang, E. C. , and Barocas, V. H. , 2006, “ Accommodative Microfluctuations and Iris Contour,” J. Vision, 6(5), pp. 10–10. [CrossRef]
Liebmann, J. M. , and Ritch, R. , 2002, “ Laser Surgery for Angle Closure Glaucoma,” Seminars in Ophthalmology, Vol. 17, Taylor & Francis, Oxford, UK, pp. 84–91. [CrossRef] [PubMed] [PubMed]
Whitcomb, J. E. , Barnett, V. A. , Olsen, T. W. , and Barocas, V. H. , 2009, “ Ex Vivo Porcine Iris Stiffening Due to Drug Stimulation,” Exp. Eye Res., 89(4), pp. 456–461. [CrossRef] [PubMed]
Whitcomb, J. E. , Amini, R. , Simha, N. K. , and Barocas, V. H. , 2011, “ Anterior–Posterior Asymmetry in Iris Mechanics Measured by Indentation,” Exp. Eye Res., 93(4), pp. 475–481. [CrossRef] [PubMed]
Reeves, P. , 1918, “ Rate of Pupillary Dilation and Contraction,” Psychol. Rev., 25(4), p. 330. [CrossRef]
Crawford, B. , 1936, “ The Dependence of Pupil Size Upon External Light Stimulus Under Static and Variable Conditions,” Proc. R. Soc. London. Ser. B, 121(823), pp. 376–395. [CrossRef]


Grahic Jump Location
Fig. 4

Iris configuration (a) before and (b) after dilation showing changes in the CL and CV

Grahic Jump Location
Fig. 3

Undeformed porcine iris model showing the axis of symmetry

Grahic Jump Location
Fig. 2

Anterior segment optical coherence tomography of a normal eye. Iris CV is defined as the largest distance between the posterior surface of the iris and the iris CL. AS-OCT image is from Dr. Vanita Pathak-Ray, FRCS (Ed), FRCOphth (Lon), Consultant Glaucoma Specialist at VST Center for Glaucoma, LV Prasad Eye Institute, Hyderabad, India with permission.

Grahic Jump Location
Fig. 1

Anteior chamber angle and iris configuration of (a) normal and (b) primary angle closure glaucoma eye (image modified from National Eye institute, NIH, Bethesda, MD). Scale bar is an approximation in the image.

Grahic Jump Location
Fig. 5

(a) Undeformed porcine iris model and (b) deformed porcine iris model. Unique features such as grooves and bumps identified in both configuration are marked by alphabets A to J.

Grahic Jump Location
Fig. 6

Inverse algorithm flowchart

Grahic Jump Location
Fig. 7

Error contour plot shows the uniqueness of the solution for a particular G and ν for noise free nodal displacement based objective function



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In