Planar biaxial testing of CXL strengthening effects.

The stiffening effect of corneal crosslinking (CXL) treatment, a therapeutic approach for managing the progression of keratoconus, has been primarily investigated using uniaxial tensile experiments. However, this testing technique has several drawbacks and is unable to measure the mechanical response of cornea under a multiaxial loading state. In this work, we used biaxial mechanical testing method to characterize biomechanical properties of porcine cornea before and after CXL treatment. We also investigated the influence of preconditioning on measured properties and used TEM images to determine microstructural characteristics of the extracellular matrix. The conventional method of CXL treatment was used for crosslinking the porcine cornea. The biaxial experiments were done by an ElectroForce TestBench system at a stretch ratio of 1:1 and a displacement rate of 2 mm/min with and without preconditioning. The experimental measurements showed no significant difference in the mechanical properties of porcine cornea along the nasal temporal (NT) and superior inferior (SI) direction. Furthermore, the CXL therapy significantly enhanced the mechanical properties along both directions without creating anisotropic response. The samples tested with preconditioning showed significantly stiffer response than those tested without preconditioning. The TEM images showed that the CXL therapy did not increase the diameter of collagen fibers but significantly decreased their interfibrillar spacing, consistent with the mechanical property improvement of CXL treated samples.

Nonlinear Mechanical Properties of Prestressed Branched Fibrous Networks.

Random fiber networks constitute the solid skeleton of many biological materials such as the cytoskeleton of cells and extracellular matrix of soft tissues. These random networks show unique mechanical properties such as nonlinear shear strain-stiffening and strain softening when subjected to preextension and precompression, respectively. In this study, we perform numerical simulations to characterize the influence of axial prestress on the nonlinear mechanical response of random network structures as a function of their micromechanical and geometrical properties. We build our numerical network models using the microstructure of disordered hexagonal lattices and quantify their nonlinear shear response as a function of uniaxial prestress strain. We consider three different material models for individual fibers and fully characterize their influence on the mechanical response of prestressed networks. Moreover, we investigate both the influence of geometric disorder keeping the network connectivity constant and the influence of the randomness in the stiffness of individual fibers keeping their mean stiffness constant. The effects of network connectivity and bending rigidity of fibers are also determined. Several important conclusions are made, including that the tensile and compressive prestress strains, respectively, increase and decrease the initial network shear stiffness but have no effect on the maximal shear modulus. We discuss the findings in terms of microstructural properties such as the local strain energy distribution.

Constitutive Modeling of Corneal Tissue: Influence of Three-Dimensional Collagen Fiber Microstructure.

The cornea, the transparent tissue in the front of the eye, along with the sclera, plays a vital role in protecting the inner structures of the eyeball. The precise shape and mechanical strength of this tissue are mostly determined by the unique microstructure of its extracellular matrix. A clear picture of the 3D arrangement of collagen fibrils within the corneal extracellular matrix has recently been obtained from the secondary harmonic generation images. However, this important information about the through-thickness distribution of collagen fibrils was seldom taken into account in the constitutive modeling of the corneal behavior. This work creates a generalized structure tensor (GST) model to investigate the mechanical influence of collagen fibril through-thickness distribution. It then uses numerical simulations of the corneal mechanical response in inflation experiments to assess the efficacy of the proposed model. A parametric study is also done to investigate the influence of model parameters on numerical predictions. Finally, a brief comparison between the performance of this new constitutive model and a recent angular integration (AI) model from the literature is given.

Mechanical properties of subisostatic random networks composed of nonlinear fibers.

Fibrous protein networks provide structural integrity to different biological materials such as soft tissues. These networks display an unusual exponential strain-stiffening behavior when subjected to mechanical loads. This nonlinear strain-stiffening behavior has so far been explained in terms of the network microstructure and the flexibility of constituting fibers. Here, we conduct a comprehensive computational study to characterize the importance of material properties of individual fibers in the overall nonlinear mechanical response of random fiber networks. To this end, we consider three nonlinear material models, ranging from an almost linear form to a highly nonlinear one, for the fibers of subisostatic disordered networks. We characterize the amount of strain-stiffening as a function of bending rigidity of the fibers, the amount of nonlinearity of the fibers, and the connectivity of random networks. We find that networks composed of highly nonlinear fibers exhibit much more strain-stiffening than networks made up of linear fibers. Furthermore, the local strain distribution becomes more homogenous as the amount of nonlinearity in the material models increases. Increasing the network connectivity signifies the importance of the nonlinear material response of individual fibers in the overall mechanical behavior of networks. The constitutive behavior of fibers plays an important role in defining the failure response of networks particularly in the damage initiation and evolution. These important findings for how the mechanical response of individual fibers affects the overall mechanical properties of random networks could find applications in designing new biomimetic materials and could help scientists better understand the mechanical properties of biological materials.

UVA/riboflavin collagen crosslinking stiffening effects on anterior and posterior corneal flaps.

The UVA/riboflavin collagen crosslinking (CXL) is one of the treatment procedure for stopping the progression of keratoconus. The inclusion criterion for this procedure is a minimum corneal thickness of 400 µm, which is not often met in patients with advanced keratoconus. Preoperatively swelling thin corneas was shown to stabilize the keratectasia without any postoperative endothelial damage. Recently, we have shown that swelling porcine corneas prior to the CXL treatment had no significant effect on the resulting improvement in their tensile properties. In the present study, we extended this previous study and characterized the stiffening effects of CXL on anterior and posterior flaps as a function of their hydration. A DSAEK system was used to excise 10 mm corneal flaps from 80 porcine corneas. Individual flaps were crosslinked at different initial hydration levels by using riboflavin solutions composed of different dextran concentrations; the thickness was taken as a measure of flap hydration. A DMA machine was used to measure the tensile properties either immediately after the CXL treatment or after the thickness (hydration) of the crosslinked samples was brought down to a specific value. The average thickness of anterior groups was 670 µm, 540 µm, and 410 µm, and the average thickness of posterior groups was 845 µm, 650 µm, and 440 µm. It was found that although CXL significantly increased the tensile properties of all anterior groups, it had an insignificant effect on the stiffness of posterior flaps. Furthermore, except for the posterior flaps in 845 µm and 650 µm thickness groups, decreasing the hydration significantly increased the tensile modulus (p < 0.05). Finally, the anterior flaps that were crosslinked at higher hydration, i.e. swollen before CXL, showed significantly less amount of stiffening in comparison with those crosslinked at lower hydration when the tensile property measurement was done at similar hydration (p < 0.05).

Collagen cross-linking treatment effects on corneal dynamic biomechanical properties.

Cornea is a soft tissue with the principal function of transmitting and refracting light rays. The objective of the current study was to characterize possible effects of the riboflavin/UVA collagen cross-linking on corneal dynamic properties. The original corneal cross-linking protocol was used to induce cross-links in the anterior portion of the bovine cornea. A DMA machine was used to conduct mechanical tensile experiments at different levels of tensile strains. The samples were divided into a control group (n = 5) and a treated group (n = 5). All specimens were first stretched to a strain of 5% and allowed to relax for twenty minutes. After completion of the stress-relaxation experiment, a frequency sweep test with oscillations ranging from 0.01 to 10 Hz was performed. The same procedure was repeated to obtain the stress-relaxation and dynamic properties at 10% strain. It was observed that the collagen cross-linking therapy significantly increased the immediate and equilibrium tensile behavior of the bovine cornea (P < 0.05). Furthermore, for all samples in control and treated groups and throughout the whole range of frequencies, a significantly larger tensile storage modulus was measured at an axial strain of 10% compared to what was obtained at a tensile strain of 5%. Finally, it was noted that although this treatment procedure resulted in a significant increase in the storage and loss modulus at any axial strain and frequency (P < 0.05), it significantly reduced the ratio of the dissipated and stored energy during a single cycle of deformation. Therefore, it was concluded that while the riboflavin/UVA collagen cross-linking increased significantly corneal stiffness, it decreased significantly its damping capability and deformability. This reduced damping ability might adversely interfere with corneal mechanical performance.

Effects of bathing solution on tensile properties of the cornea.

The cornea is a transparent tissue with the major functions of protecting the inner contents of the eye and refracting incoming light. The biomechanical properties of the cornea strongly depend on the microstructure and composition of the stromal layer, a hydrated bio-gel. The uniaxial strip testing is a convenient and well-accepted experimental technique for characterizing corneal material parameters. It is known that the water content of specimens in this method depends on the osmolality of the bathing solution. The present study was designed to investigate the effects of different bathing solutions on uniaxial tensile material properties of the cornea. The tensile behavior of bovine corneal samples was measured in six different bathing solutions, i.e., hypertonic solution (12% NaCl solution), common preserving isotonic solutions (e.g., phosphate buffer saline, ophthalmic balanced salt solution, and 0.9% NaCl solution), hypotonic solution (distilled water), and neutral solution (mineral oil). It was observed that the bathing solution had significant influence on the tensile behavior of the corneal samples. In particular, the specimens tested in bathing solutions causing less swelling had significantly stiffer tensile properties. Furthermore, a simple mathematical model based on Voigt composite material model was developed to represent the measured solution-dependent tensile properties. The present study suggests that extra attention should be paid to corneal thickness (hydration) in uniaxial tensile experiments. It also provides important data on tensile properties of the cornea; such information could significantly contribute to improving the accuracy of numerical predictions of corneal biomechanics.

Regional differences in electroactive response of the sclera.

The sclera exhibits mechanical response when subjected to an external electric stimulation. The scleral electroactive response is a function of its charge density, mechanical properties, thickness, and strength of the applied electric voltage. The primary objective of the present work was to investigate the regional differences in the electroactive response of porcine sclera. To this end, we cut scleral strips in meridional directions from superior-temporal, superior-nasal, inferior-temporal, and inferior-nasal quadrants. In addition, we excised samples circumferentially from the posterior, equatorial, and anterior regions. The electroactive bending response of these samples was measured under 10 and 15 V in 0.15 M NaCl solution. The meridional samples were tested under two different configurations by clamping them either from their anterior or posterior end. It was observed that the scleral electroactive deformation increased with increasing the the electric voltage. Furthermore, regardless of the region from which meridional strips were excised, their electroactive response was considerably larger when they were clamped from their anterior end. Unlike meridional strips, the electroactive response of circumferential samples was significantly dependent on the location, that is, the average maximum bending angle of posterior samples was significantly larger than that of equatorial and anterior strips. The regionally different electroactive bending response of the sclera was discussed in terms of the variation in its biochemical and biomechanical properties throughout the eyeball.

Hydration dependent biomechanical properties of the corneal stroma.

The cornea is responsible for about seventy percent of refractive power of the eye and transmits more than ninety percent of the incident light. The refractive power and transparency of the cornea are largely contingent upon active maintenance of its precise curvature and highly regular microstructure. The biomechanical properties of the cornea are mainly derived from the stromal layer. Over past decades, many computational models have been proposed to predict the corneal behavior in normal and/or diseased state. The predictions of these numerical methods strongly depend on accurate description of corneal mechanical properties. The present study used unconfined compression technique to characterize the dependence of corneal material parameters on thickness variations due to hydration/dehydration. A series of unconfined compression tests was performed on porcine corneal buttons and stress relaxation response of the samples was obtained. A transversely isotropic biphasic model was used to analyze experimental measurements at each ramp-and-hold step. The in-plane Young's modulus, out-of-plane (transverse) Young's modulus, and permeability coefficient of the samples were determined as a function of average thickness. The average thickness variation due to swelling (or dehydration) was between 0.69 mm and 1.27 mm. It was found that corneal material parameters depend strongly on the tissue thickness (hydration). In particular, the in-plane elastic modulus increased with decreasing the thickness (p < 0.05) and the permeability coefficient decreased with decreasing the thickness (p < 0.05). Furthermore, it was observed that although the out-of-plane modulus was almost constant as the average thickness varied between 1.07 mm and 1.27 mm, it increased with tissue thickness outside this range (p < 0.05). The findings of this study provide important information on biomechanical properties of corneal stroma and are useful in computational simulations of the cornea.

On the Mechanical Roles of Glycosaminoglycans in the Tensile Properties of Porcine Corneal Stroma.

Purpose: The cornea is primarily composed of collagen fibrils that are embedded in a ground substance rich in proteoglycans and other glycoproteins. It is known that glycosaminoglycan (GAG) side chains of proteoglycans form anti-parallel duplexes between collagen fibrils. The present work was done in order to investigate the mechanical function of GAGs in defining the tensile properties of porcine corneal stroma. Methods: Porcine corneal stromal strips dissected from the nasal-temporal direction were divided into control, buffer-treated, and enzyme-treated groups. The samples in the control group were used immediately after dissection. However, the buffer-treated and enzyme-treated samples were, respectively, incubated for 18 hours at 37°C in a buffer solution made up of 100-mM sodium acetate at pH 6.0 or in an enzyme solution containing keratanase II. The Blyscan assay was used to quantify the total GAG content and assess GAG depletion in the samples treated with the enzyme and buffer solutions. Uniaxial tensile tests were also performed to determine the effect of GAG removal on mechanical properties of the cornea. Results: The GAG content in enzyme-treated samples was significantly lower than that of the normal and buffer-treated specimens (P < 0.05). Moreover, GAG-depleted strips showed significantly softer mechanical responses in comparison with the control and buffer samples (P < 0.05). Conclusions: Removing GAGs from the corneal extracellular matrix led to significant tensile property reduction; supporting the hypothesis that there exists a strong correlation between the GAG content and mechanical properties of the corneal stroma.

Finite Deformation of Scleral Tissue under Electrical Stimulation: An Arbitrary Lagrangian-Eulerian Finite Element Method.

The sclera is considered as the principal load-bearing tissue within the eye. The sclera is negatively charged; thus, it exhibits mechanical response to electrical stimulation. We recently demonstrated the electroactive behavior of sclera by performing experimental measurements that captured the deformation of the tip of scleral strips subjected to electric voltage. We also numerically analyzed the electromechanical response of the tissue using a chemo-electro-mechanical model. In the pre-sent study, we extended our previous work by experimentally characterizing the deformation profile of scleral strips along their length under electrical stimulation. In addition, we improved our previous mathematical model such that it could numerically capture the large deformation of samples. For this purpose, we considered the transient variability of the fixed charge density and the coupling between mechanical and chemo-electrical phenomena. These improvements in-creased the accuracy of the computational model, resulting in a better numerical representation of experimentally measured bending angles.

Modeling and experimental investigation of electromechanical properties of scleral tissue; a CEM model using an anisotropic hyperelastic constitutive relation.

The sclera is a soft tissue primarily consisting of collagen fibers, elastin, and proteoglycans. The proteoglycans are composed of a core protein and negatively charged glycosaminoglycan side chains. The fixed electric charges inside the scleral extracellular matrix play a key role in its swelling and are expected to cause the tissue to deform in response to an electric field. However, the electroactive response of the sclera has not yet been investigated. The present work experimentally demonstrates that sclera behaves similar to an anionic electrosensitive hydrogel and develops a chemo-electro-mechanical (CEM) mathematical framework for its electromechanical response. In the numerical model, a hyperelastic constitutive law with distributed collagen fibers is used to capture the nonlinear mechanical properties of the sclera, and the coupled Poisson-Nernst-Planck equations represent the distribution of mobile ions throughout the domain. After calibrating the proposed numerical CEM model against the experimental measurements, we employ it to investigate the effects of different parameters on the scleral electromechanical response including the voltage and fixed charge density. The experimental and numerical findings of the present study confirm that sclera behaves as an electroactive hydrogel and provide new insight into the mechanical response of this ocular tissue.

Effect of corneal collagen crosslinking on viscoelastic shear properties of the cornea.

The cornea is responsible for most of the refractive power in the eye and acts as a protective layer for internal contents of the eye. The cornea requires mechanical strength for maintaining its precise shape and for withstanding external and internal forces. Corneal collagen crosslinking (CXL) is a treatment option to improve corneal mechanical properties. The primary objective of this study was to characterize CXL effects on viscoelastic shear properties of the porcine cornea as a function of compressive strain. For this purpose, corneal buttons were prepared and divided into three groups: control group (n = 5), pseudo-crosslinked group (n = 5), and crosslinked group (n = 5). A rheometer was used to perform dynamics torsional shear experiments on corneal disks at different levels of compressive strain (0%-40%). Specifically, strain sweep experiments and frequency sweep tests were done in order to determine the range of linear viscoelasticity and frequency dependent shear properties, respectively. It was found that the shear properties of all samples were dependent on the shear strain magnitude, loading frequency, and compressive strain. With increasing the applied shear strain, all samples showed a nonlinear viscoelastic response. Furthermore, the shear modulus of samples increased with increasing the frequency of the applied shear strain and/or increasing the compressive strain. Finally, the CXL treatment significantly increased the shear storage and loss moduli when the compressive strain was varied from 0% to 30% (p < 0.05); larger shear moduli were observed at compressive 40% strain but the difference was not significant (P = 0.12).

Experimental and numerical analysis of electroactive characteristics of scleral tissue.

The sclera provides mechanical support to retina and protects internal contents of the eye against external injuries. The scleral extracellular matrix is mainly composed of collagen fibers and proteoglycans (PGs). At physiological pH, collagen molecules are neutral but PGs contain negatively charged glycosaminoglycan chains. Thus, the sclera can be considered as a polyelectrolyte hydrogel and is expected to exhibit mechanical response when subjected to electrical stimulations. In this study, we mounted scleral strips, dissected from the posterior part of porcine eyes, at the center of a custom-designed container between two electrodes. The container was filled with NaCl solution and the bending deformation of scleral strips as a function of the applied electric voltage was measured experimentally. It was found that scleral strips reached to an average bending angle of 3°, 10° and 23° when subjected to 5V, 10V, and 15V, respectively. We also created a chemo-electro-mechanical finite element model for simulating the experimental measurements by solving coupled Poisson-Nernst-Plank and equilibrium mechanical field equations. The scleral fixed charge density and modulus of elasticity were found by fitting the experimental data. The ion concentration distribution inside the domain was found numerically and was used to explain the underlying mechanisms for the scleral electroactive response. The numerical simulations were also used to investigate the effects of various parameters such as the electric voltage and fixed charge density on the scleral deformation under an electric field. STATEMENT OF SIGNIFICANCE: This manuscript investigates the electroactive response of scleral tissue. It demonstrates that the sclera deforms mechanically when subjected to electrical stimulations. A chemo-electro-mechanical model is also presented in order to numerically capture the electromechanical response of the sclera. This numerical model is used to explain the experimental observations by finding the ion distribution inside the tissue under an electric field. This work is significant because it shows that the sclera is an electroactive polyanionic hydrogel and it provides new information about the underlying mechanisms governing its mechanical and electrical properties.

Regional Differences in the Glycosaminoglycan Role in Porcine Scleral Hydration and Mechanical Behavior.

Purpose: This study characterized the role of glycosaminoglycans (GAGs) in the hydration, thickness, and biomechanical properties of posterior and anterior porcine sclera. Methods: The scleral discs and strips were obtained from the anterior and posterior parts of porcine eyes, and their initial hydration and thickness were measured. The anterior and posterior scleral discs were used to show the efficacy of the GAG removal protocol by quantifying their GAG content. The strips were divided into three groups of PBS treatment, buffer treatment, and enzyme treatment in order to assess the effects of different treatment procedures on the thickness, hydration, and viscoelastic properties of the samples. The mechanical properties of the strips were determined by performing uniaxial tensile stress relaxation experiments. Results: It was found that the control and buffer groups had insignificant differences in all measured quantities. The samples from the posterior region had a significantly larger GAG content and thickness in comparison with those from anterior region; however, there was an insignificant difference in their hydration. The GAG depletion process decreased the hydration of both anterior and posterior samples significantly (P < 0.05). Furthermore, the mechanical tests showed that the removal of GAGs resulted in stiffer mechanical behavior in both anterior and posterior samples (P < 0.05). In particular, the peak stress and equilibrium stress were significantly larger for the strips in the enzyme treatment group. Conclusions: GAGs and their interaction with the collagen network are important in defining the hydration and mechanical properties of both posterior and anterior sclera.

Tensile Viscoelastic Properties of the Sclera after Glycosaminoglycan Depletion.

PURPOSE: Fibrillar collagen network and glycosaminoglycans (GAGs) are the primary components of extracellular matrix (ECM) of the sclera. The main goal of this study was to investigate the possible structural roles of GAGs in the scleral tensile properties as a function of preconditioning and displacement rate. METHODS: Four-step uniaxial stress-relaxation tests were used for characterizing the viscoelastic tensile response of the posterior porcine sclera with and without enzymatic GAG removal. The scleral strips were divided into different groups based on the displacement rate and the presence or absence of a preconditioning step in the loading protocol. The groups were (1) displacement rate of 0.2 mm/min without preconditioning, (2) displacement rate of 1 mm/min without preconditioning, (3) displacement rate of 0.2 mm/min with preconditioning, and (4) displacement rate of 1 mm/min with preconditioning. The peak stress, equilibrium stress, and the equilibrium elastic modulus were calculated for all specimens and compared against each other. RESULTS: Increasing the displacement rate from 0.2 mm/min to 1.0 mm/min was found to cause an insignificant change in the equilibrium stress and equilibrium elastic modulus of porcine scleral strips. Removal of GAGs resulted in an overall stiffer tensile behavior independent of the displacement rate in samples that were not preconditioned (P < .05). The behavior of preconditioned samples with and without GAG removal was not significantly different from each other. CONCLUSIONS: The experimental measurements of the present study showed that GAGs play an important role in the mechanical properties of the posterior porcine sclera. Furthermore, using a preconditioning step in the uniaxial testing protocol resulted in not being able to identify any significant difference in the tensile behavior of GAG depleted and normal scleral strips.

A computational study of the behavior of colloidal gel networks at low volume fraction.

Colloidal gel networks appear in different scientific and industrial applications because of their unique properties. Molecular dynamics simulations could reveal the relation between molecular level and macroscopic properties of these systems. Nevertheless, the predictions of numerical simulations might depend on the specific form and parameters of interaction potentials. In this paper, a new effective interaction potential is used for characterizing the mechanical behavior of low volume fraction colloidal gels under large shear deformation. The findings are compared with those obtained from other available forms of interaction potentials in order to determine gel characteristics that are interaction potential independent. Furthermore, the macroscopic stress-strain behavior is discussed in terms of the behavior of different terms of the proposed interaction potential. The correlation between the stretch of interparticle bonds and their alignment in the direction of the maximum principal stress is also computed in order to provide microscopic explanations for the initial strain softening behavior. It is concluded that, in addition to topology, local mechanical interactions between colloidal particles are important in defining the mechanical response of soft gels.

Hydration related changes in tensile response of posterior porcine sclera.

It has been shown that there exists significance dependence between hydration and biomechanical properties of hydrated tissues such as cornea. The primary purpose of this study was to determine hydration effects on mechanical properties of sclera. Scleral strips, dissected from the posterior part of pig eyes along the superior-inferior direction, were divided into four hydration groups by first drying them and then soaking them in PBS until their hydration reached to 75%, 100%, 150%, and 200%. The strips were subjected to ten consecutive cycles of loading and unloading up to 1 MPa. The response of samples at the tenth cycle was used to compute the tangent modulus, maximum strain, and hysteresis as a function of hydration. The experiments were done in oil in order to prevent hydration changes during the mechanical tests. The mechanical response of strips right after dissection, control group, was also measured. In general, significant softening of scleral strips was found with increasing hydration (p < 0.05). The stress-strain response of control group was between those of samples with hydration 150% and 200%. The experimental stress-strain data were successfully represented numerically with an exponential mathematical relation with R2 > 0.99. The present study showed that hydration would significantly alter the tensile response of scleral tissue. Thus, the hydration of scleral specimens during mechanical experimental measurements should be carefully controlled.

The contribution of sGAGs to stress-controlled tensile response of posterior porcine sclera.

Despite the significant progress in characterizing mechanical functions of individual scleral extracellular matrix (ECM) components, the biomechanical contribution of sulfated glycosaminoglycans (sGAGs) is still poorly understood. The primary purpose of this study was to determine the possible function of sGAGs in scleral mechanical response by characterizing the tensile behavior of normal and sGAG-depleted samples. We used chondroitinase ABC solution to remove sGAGs from scleral samples that were dissected from posterior porcine eyes. We performed biochemical analyses for assessing the efficacy of sGAG removal protocol. Furthermore, we conducted stress-controlled uniaxial tensile tests to characterize the influence of sGAG removal on mechanical properties of sclera. The tensile behavior of scleral strips right after dissection and after being soaked in buffer was also determined. Biochemical analyses confirmed that 18 hour incubation in 0.125 U/ml Chondroitinase ABC solution removed over 90% of chondroitin and dermatan sGAGs. No significant difference was observed in the thickness/hydration of samples because of enzyme- and buffer-treated samples. Furthermore, it was found that sGAG depletion did not significantly alter the tangent modulus, energy dissipation, and peak strain of posterior scleral strips. It was concluded that sGAGs did not influence the stress-controlled viscoelastic tensile response of sclera.

Colloidal particle gel models using many-body potential interactions.

Many-body effective interactions are commonly used in a molecular dynamics simulation study of gel networks formed by colloidal particles. Here we report an interaction potential that can be used to investigate the mechanical response of colloidal gel networks under shear deformation. We then investigate the dependence of the numerical simulation results on the form of mathematical expression used to define the interparticle interactions. This work reveals insight into particle gel models by discussing the physical origins of their mechanical response.