Research Papers

Collagen Structure and Mechanical Properties of the Human Sclera: Analysis for the Effects of Age

[+] Author and Article Information
Baptiste Coudrillier

Wallace H. Coulter Department of Biomedical
Georgia Institute of Technology and Emory University,
Atlanta, GA 30332
e-mail: baptiste.coudrillier@bme.gatech.edu

Jacek Pijanka, Craig Boote

Structural Biophysics Group,
School of Optometry and Vision Sciences,
Cardiff University,
Cardiff CF24 4HQ, Wales, UK

Joan Jefferys, Harry A. Quigley

Glaucoma Center of Excellence,
Wilmer Ophthalmological Institute,
Johns Hopkins University School of Medicine,
Baltimore, MD 21287

Thomas Sorensen

Diamond Light Source,
Didcot, Oxfordshire, UK

Thao D. Nguyen

Department of Mechanical Engineering,
Johns Hopkins University,
Baltimore, MD 21218
e-mail: vicky.nguyen@jhu.edu

Manuscript received January 20, 2014; final manuscript received December 15, 2014; published online February 11, 2015. Assoc. Editor: Jonathan Vande Geest.

J Biomech Eng 137(4), 041006 (Feb 11, 2015) (14 pages) Paper No: BIO-14-1028; doi: 10.1115/1.4029430 History: Received January 20, 2014

The objective of this study was to measure the collagen fiber structure and estimate the material properties of 7 human donor scleras, from age 53 to 91. The specimens were subjected to inflation testing, and the full-field displacement maps were measured by digital image correlation. After testing, the collagen fiber structure was mapped using wide-angle X-ray scattering. A specimen-specific inverse finite element method was applied to calculate the material properties of the collagen fibers and interfiber matrix by minimizing the difference between the experimental displacements and model predictions. Age effects on the fiber structure and material properties were estimated using multivariate models accounting for spatial autocorrelation. Older age was associated with a larger matrix stiffness (p = 0.001), a lower degree of fiber alignment in the peripapillary sclera (p = 0.01), and a lower mechanical anisotropy in the peripapillary sclera (p = 0.03).

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



Grahic Jump Location
Fig. 1

Angular X-ray scatter intensity profile for a WAXS measurement. The total scatter may be separated into that arising from isotropically arranged collagen fibers (shaded region) and that arising from preferentially oriented fibers (hatched region). The degree of fiber alignment was computed as the ratio of the aligned scatter to the total scatter.

Grahic Jump Location
Fig. 2

(a) Top view of the surface layer of the finite element (FE) mesh used for the IFEA, showing the nodes at which DIC-measured displacements were applied as kinematic boundary conditions. The outer line represents the location of the most peripheral pixels of the sclera detected by the DIC software. The dashed circle is the border between the generic and specimen-specific mesh. (b) Schematic of the transverse view of the finite element mesh showing the thickness profile in the generic model of the peripapillary sclera. The LC and the remant tissues were not included in the FE model but are represented for clarity. In the inverse analysis, we assumed that the displacements on the surface of the generic peripapillary model were identical to the DIC-measured displacements on the surface of the remnant tissues. The geometry was discretized using mixed Q1P0 trilinear hexahedral elements, with three elements spanning the thickness. The mesh had 40 nodes in the meridional direction and 28 nodes in the circumferential direction. The node density was larger toward the LC to capture the stress concentration due to the compliant tissue of the ONH.

Grahic Jump Location
Fig. 3

Flow diagram describing the search algorithm for the value of the matrix stiffness μ and fiber parameters (α, β)

Grahic Jump Location
Fig. 4

(a) Finite element model used to validate that the applied kinematic boundary conditions did not overconstrain the solution of the IFEA. In this model, the LC was modeled as a neo-Hookean material, with shear modulus equal to a tenth of that of the sclera. (b) We used the same mesh for the sclera in the inverse method. Forward displacements were applied to the edges of the mesh (base and scleral canal).

Grahic Jump Location
Fig. 5

(a) Representation of the single element used to calculate the MA. The fiber structure was described using a single WAXS measurement. A biaxial stretch corresponding to the average strain measured in Coudrillier et al. [18] was applied on the ef and ep faces. (b) Representation of the Green–Lagrange strain/Cauchy stress response. The MA was calculated using Eq. (9).

Grahic Jump Location
Fig. 6

(a) A composite polar plot showing the preferred orientations of aligned collagen fibers for the specimen FC53r. The color scale conveys the degree of fiber alignment. FC53r was representative of the other specimens, showing a strong circumferential fiber alignment in the peripapillary sclera, a much more isotropic midposterior sclera, and two symmetrical oblique features emerging from the temporal pole of the peripheral peripapillary sclera, indicated with arrows. (b) Contour map of calculated degree of fiber alignment, defined as the ratio of the aligned scatter to the total scatter as shown in Fig. 1. The degree of fiber alignment was largest in the temporal/superior quadrant and lowest in the superior/nasal quadrant of the peripapillary sclera.

Grahic Jump Location
Fig. 7

Degree of fiber alignment averaged over each quadrant. Age was associated with a significant decrease in mean degree of fiber alignment (p = 0.01).

Grahic Jump Location
Fig. 8

(a) Contour of the cost function in the (μ, γ)-space with wn = 1 for each node (Eq. (2)). (b) Contour of the cost function in the (μ, γ)-space with wn chosen to reflect the degree of anisotropy of the closest WAXS measurement. The weights are represented in the figure on the bottom left. The contribution of the nodes of the peripapillary sclera to the cost function was larger than those of the midposterior sclera. The cost function is represented between its minimum value m and 1.1 m.

Grahic Jump Location
Fig. 9

(a) Step 1, (b) step 3, and (c) step 4 of the sequential algorithm as shown in Fig. 3 and explained in Sec. 2.6. (d) Map of cost function for a dense 3D grid. We evaluated 21,519 combinations of material parameters (α, β, μ), with 0 < α < 200 kPa, 0 < β < 150, and 30 < μ < 1000 kPa. Although hard to see in this figure, the global search algorithm confirmed that the restriction for μ obtained in our sequential algorithm was accurate.

Grahic Jump Location
Fig. 10

Stress–strain curves in the fiber (a) and perpendicular (b) directions obtained from an equibiaxial stress to 60 kPa, corresponding to the hoop stress in a thin-walled sphere of radius 12 mm and thickness 1.2 mm at a pressure of 90 mm Hg. The fiber structure was represented by a von Mises function, which was fit to the specimen-specific averaged fiber distribution in the peripapillary sclera.

Grahic Jump Location
Fig. 11

(a) Matrix shear modulus μ and (b) fiber stiffness 4αβ plotted versus age

Grahic Jump Location
Fig. 12

Map of the MA defined in Eq. (9) for FC77r. The MA was largest in the temporal/superior quadrant of the peripapillary sclera and lowest in the superior/nasal quadrant of the peripapillary sclera.

Grahic Jump Location
Fig. 13

Average MA in each quadrant of the peripapillary sclera plotted versus age




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