A Combination of Constitutive Damage Model and Artificial Neural Networks to Characterize the Mechanical Properties of the Healthy and Atherosclerotic Human Coronary Arteries.

It has been indicated that the content and structure of the elastin and collagen of the arterial wall can subject to a significant alteration due to the atherosclerosis. Consequently, a high tissue stiffness, stress, and even damage/rupture are triggered in the arterial wall. Although many studies so far have been conducted to quantify the mechanical properties of the coronary arteries, none of them consider the role of collagen damage of the healthy and atherosclerotic human coronary arterial walls. Recently, a fiber family-based constitutive equation was proposed to capture the anisotropic mechanical response of the healthy and atherosclerotic human coronary arteries via both the histostructural and uniaxial data. In this study, experimental mechanical measurements along with histological data of the healthy and atherosclerotic arterial walls were employed to determine the constitutive damage parameters and remodeling of the collagen fibers. To do this, the preconditioned arterial tissues were excised from human cadavers within 5-h postmortem, and the mean angle of their collagen fibers was precisely determined. Thereafter, a group of quasistatic axial and circumferential loadings were applied to the arterial walls, and the constrained nonlinear minimization method was employed to identify the arterial parameters according to the axial and circumferential extension data. The remodeling of the collagen fibers during the tensile test was also predicted via Artificial Neural Networks algorithm. Regardless of loading direction, the results presented a noteworthy load-bearing capability and stiffness of the atherosclerotic arteries compared to the healthy ones (P < 0.005). Theoretical fiber angles were found to be consistent with the experimental histological data with less than 2 and 5° difference for the healthy and atherosclerotic arterial walls, respectively. The pseudoelastic damage model data were also compared with that of the experimental data, and interestingly, the arterial mechanical behavior for both the primary loading (up to the elastic region) and the discontinuous softening (up to the ultimate stress) was well addressed. The proposed model predicted well the mechanical response of the arterial tissue considering the damage of collagen fibers for both the healthy and atherosclerotic arterial walls.

Segmental biomechanics of the normal and glaucomatous human aqueous outflow pathway.

The conventional aqueous outflow pathway, encompassing the trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and inner wall endothelium of Schlemm's canal (SC), governs intraocular pressure (IOP) regulation. This study targets the biomechanics of low-flow (LF) and high-flow (HF) regions within the aqueous humor outflow pathway in normal and glaucomatous human donor eyes, using a combined experimental and computational approach. LF and HF TM/JCT/SC complex tissues from normal and glaucomatous eyes underwent uniaxial tensile testing. Dynamic motion of the TM/JCT/SC complex was recorded using customized green-light optical coherence tomography during SC pressurization in cannulated anterior segment wedges. A hyperviscoelastic model quantified TM/JCT/SC complex properties. A fluid-structure interaction model simulated tissue-aqueous humor interaction. FluoSpheres were introduced into the pathway via negative pressure in the SC, with their motion tracked using two-photon excitation microscopy. Tensile test results revealed that the elastic moduli of the LF and HF regions in glaucomatous eyes are 3.5- and 1.5-fold stiffer than the normal eyes, respectively. The FE results also showed significantly larger shear moduli in the TM, JCT, and SC of the glaucomatous eyes compared to the normal subjects. The LF regions in normal eyes demonstrated larger elastic moduli compared to the HF regions in glaucomatous eyes. The resultant strain in the outflow tissues and velocity of the aqueous humor in the FSI models were in good agreement with the digital volume correlation and 3D particle image velocimetry data, respectively. This study uncovers stiffer biomechanical responses in glaucomatous eyes, with LF regions stiffer than HF regions in both normal and glaucomatous eyes. STATEMENT OF SIGNIFICANCE: This study delves into the biomechanics of the conventional aqueous outflow pathway, a crucial regulator of intraocular pressure and ocular health. By analyzing mechanical differences in low-flow and high-flow regions of normal and glaucomatous eyes, this research unveils the stiffer response in glaucomatous eyes. The distinction between regions' properties offers insights into aqueous humor outflow regulation, while the integration of experimental and computational methods enhances credibility. These findings have potential implications for disease management and present a vital step toward innovative ophthalmic interventions. This study advances our understanding of glaucoma's biomechanical basis and its broader impact on ocular health.

Implementing new computational methods for the study of JCT and SC inner wall basement membrane biomechanics and hydrodynamics.

PURPOSE: The conventional aqueous outflow pathway, which includes the trabecular meshwork (TM), juxtacanalicular tissue (JCT), and the inner wall endothelium of Schlemm's canal (SC), regulates intraocular pressure (IOP) by controlling the aqueous humor outflow resistance. Despite its importance, our understanding of the biomechanics and hydrodynamics within this region remains limited. Fluid-structure interaction (FSI) offers a way to estimate the biomechanical properties of the JCT and SC under various loading and boundary conditions, providing valuable insights that are beyond the reach of current imaging techniques. METHODS: In this study, a normal human eye was fixed at a pressure of 7 mm Hg, and two radial wedges of the TM tissues, which included the SC inner wall basement membrane and JCT, were dissected, processed, and imaged using 3D serial block-face scanning electron microscopy (SBF-SEM). Four different sets of images were used to create 3D finite element (FE) models of the JCT and inner wall endothelial cells of SC with their basement membrane. The outer JCT portion was carefully removed as the outflow resistance is not in that region, leaving only the SCE inner wall and a few µm of the tissue, which does contain the resistance. An inverse iterative FE algorithm was then utilized to calculate the unloaded geometry of the JCT/SC complex at an aqueous humor pressure of 0 mm Hg. Then in the model, the intertrabecular spaces, pores, and giant vacuole contents were replaced by aqueous humor, and FSI was employed to pressurize the JCT/SC complex from 0 to 15 mm Hg. RESULTS: In the JCT/SC complex, the shear stress of the aqueous humor is not evenly distributed. Areas proximal to the inner wall of SC experience larger stresses, reaching up to 10 Pa, while those closer to the JCT undergo lower stresses, approximately 4 Pa. Within this complex, giant vacuoles with or without I-pore behave differently. Those without I-pores experience a more significant strain, around 14%, compared to those with I-pores, where the strain is roughly 9%. CONCLUSIONS: The distribution of aqueous humor wall shear stress is not uniform within the JCT/SC complex, which may contribute to our understanding of the underlying selective mechanisms in the pathway.

Developing an experimental-computational workflow to study the biomechanics of the human conventional aqueous outflow pathway.

The aqueous humor actively interacts with the trabecular meshwork (TM), juxtacanalicular tissue (JCT), and Schlemm's canal (SC) through a dynamic fluid-structure interaction (FSI) coupling. Despite the fact that intraocular pressure (IOP) undergoes significant fluctuations, our understanding of the hyperviscoelastic biomechanical properties of the aqueous outflow tissues is limited. In this study, a quadrant of the anterior segment from a normal human donor eye was dynamically pressurized in the SC lumen, and imaged using a customized optical coherence tomography (OCT). The TM/JCT/SC complex finite element (FE) with embedded collagen fibrils was reconstructed based on the segmented boundary nodes in the OCT images. The hyperviscoelastic mechanical properties of the outflow tissues' extracellular matrix with embedded viscoelastic collagen fibrils were calculated using an inverse FE-optimization method. Thereafter, the 3D microstructural FE model of the TM, with adjacent JCT and SC inner wall, from the same donor eye was constructed using optical coherence microscopy and subjected to a flow load-boundary from the SC lumen. The resultant deformation/strain in the outflow tissues was calculated using the FSI method, and compared to the digital volume correlation (DVC) data. TM showed larger shear modulus (0.92 MPa) compared to the JCT (0.47 MPa) and SC inner wall (0.85 MPa). Shear modulus (viscoelastic) was larger in the SC inner wall (97.65 MPa) compared to the TM (84.38 MPa) and JCT (56.30 MPa). The conventional aqueous outflow pathway is subjected to a rate-dependent IOP load-boundary with large fluctuations. This necessitates addressing the biomechanics of the outflow tissues using hyperviscoelastic material-model. STATEMENT OF SIGNIFICANCE: While the human conventional aqueous outflow pathway is subjected to a large-deformation and time-dependent IOP load-boundary, we are not aware of any studies that have calculated the hyperviscoelastic mechanical properties of the outflow tissues with embedded viscoelastic collagen fibrils. A quadrant of the anterior segment of a normal humor donor eye was dynamically pressurized from the SC lumen with relatively large fluctuations. The TM/JCT/SC complex were OCT imaged and the mechanical properties of the tissues with embedded collagen fibrils were calculated using the inverse FE-optimization algorithm. The resultant displacement/strain in the FSI outflow model was validated versus the DVC data. The proposed experimental-computational workflow may significantly contribute to understanding of the effects of different drugs on the biomechanics of the conventional aqueous outflow pathway.

Finite element modeling of the complex anisotropic mechanical behavior of the human sclera and pia mater.

BACKGROUND AND OBJECTIVE: Accurate finite element (FE) simulation of the optic nerve head (ONH) depends on accurate mechanical properties of the load-bearing tissues. The peripapillary sclera in the ONH exhibits a depth-dependent, anisotropic, heterogeneous collagen fiber distribution. This study proposes a novel cable-in-solid modeling approach that mimics heterogeneous anisotropic collagen fiber distribution, validates the approach against published experimental biaxial tensile tests of scleral patches, and demonstrates its effectiveness in a complex model of the posterior human eye and ONH. METHODS: A computational pipeline was developed that defines control points in the sclera and pia mater, distributes the depth-dependent circumferential, radial, and isotropic cable elements in the sclera and pia in a pattern that mimics collagen fiber orientation, and couples the cable elements and solid matrix using a mesh-free penalty-based cable-in-solid algorithm. A parameter study was performed on a model of a human scleral patch subjected to biaxial deformation, and computational results were matched to published experimental data. The new approach was incorporated into a previously published eye-specific model to test the method; results were then interpreted in relation to the collagen fibers' (cable elements) role in the resultant ONH deformations, stresses, and strains. RESULTS: Results show that the cable-in-solid approach can mimic the full range of scleral mechanical behavior measured experimentally. Disregarding the collagen fibers/cable elements in the posterior eye model resulted in ∼20-60% greater tensile and shear stresses and strains, and ∼30% larger posterior deformations in the lamina cribrosa and peripapillary sclera. CONCLUSIONS: The cable-in-solid approach can easily be implemented into commercial FE packages to simulate the heterogeneous and anisotropic mechanical properties of collagenous biological tissues.

Relative Contributions of Intraocular and Cerebrospinal Fluid Pressures to the Biomechanics of the Lamina Cribrosa and Laminar Neural Tissues.

Purpose: The laminar region of the optic nerve head (ONH), thought to be the site of damage to the retinal ganglion cell axons in glaucoma, is continuously loaded on its anterior and posterior surfaces by dynamic intraocular pressure (IOP) and orbital cerebrospinal fluid pressure (CSFP), respectively. Thus, translaminar pressure (TLP; TLP = IOP-CSFP) has been proposed as a glaucoma risk factor. Methods: Three eye-specific finite element models of the posterior human eye were constructed, including full 3D microstructures of the load-bearing lamina cribrosa (LC) with interspersed laminar neural tissues (NTs), and heterogeneous, anisotropic, hyperelastic material formulations for the surrounding peripapillary sclera and adjacent pia. ONH biomechanical responses were simulated using three combinations of IOP and CSFP loadings consistent with posture change from sitting to supine. Results: Results show that tensile, compressive, and shear stresses and strains in the ONH were higher in the supine position compared to the sitting position (P < 0.05). In addition, LC beams bear three to five times more TLP-driven stress than interspersed laminar NT, whereas laminar NT exhibit three to five times greater strain than supporting LC (P < 0.05). Compared with CSFP, IOP drove approximately four times greater stress and strain in the LC, NT, and peripapillary sclera, normalized per mm Hg pressure change. In addition, IOP drove approximately three-fold greater scleral canal expansion and anterior-posterior laminar deformation than CSFP per mm Hg (P < 0.05). Conclusions: Whereas TLP has been hypothesized to play a prominent role in ONH biomechanics, the IOP and CSFP effects are not equivalent, as IOP-driven stress, strain, and deformation play a more dominant role than CSFP effects.

Modeling the biomechanics of the conventional aqueous outflow pathway microstructure in the human eye.

BACKGROUND AND OBJECTIVE: Intraocular pressure (IOP) is determined by aqueous humor outflow resistance, which is a function of the combined resistance of Schlemm's canal (SC) endothelium and the trabecular meshwork (TM) and their interactions in the juxtacanalicular connective tissue (JCT) region. Aqueous outflow in the conventional outflow pathway results in pressure gradient across the TM, JCT, and SC inner wall, and induces mechanical stresses and strains that influence the geometry and homeostasis of the outflow system. The outflow resistance is affected by alteration in tissues' geometry, so there is potential for active, two-way, fluid-structure interaction (FSI) coupling between the aqueous humor (fluid) and the TM, JCT, and SC inner wall (structure). However, our understanding of the biomechanical interactions of the aqueous humor with the outflow connective tissues and its contribution to the outflow resistance regulation is incomplete. METHODS: In this study, a microstructural finite element (FE) model of a human eye TM, JCT, and SC inner wall was constructed from a segmented, high-resolution histologic 3D reconstruction of the human outflow system. Three different elastic moduli (0.004, 0.128, and 51.5 MPa based on prior reports) were assigned to the TM/JCT complex while the elastic modulus of the SC inner wall was kept constant at 0.00748 MPa. The hydraulic conductivity was programmed separately for the TM, JCT, and SC inner wall using a custom subroutine. Cable elements were embedded into the TM and JCT extracellular matrix to represent the directional stiffness imparted by anisotropic collagen fibril orientation. The resultant stresses and strains in the outflow system were calculated using fluid-structure interaction method. RESULTS: The higher TM/JCT stiffness resulted in larger stresses, but smaller strains in the outflow connective tissues, and resulted in a 4- and 5-fold larger pressure drop across the SC inner wall, respectively, compared to the most compliant model. Funneling through µm-sized SC endothelial pores was evident in the models at lower tissue stiffness, but aqueous flow was more turbulent in models with higher TM/JCT stiffness. CONCLUSIONS: The mechanical properties of the outflow tissues play a crucial role in the hydrodynamics of the aqueous humor in the conventional outflow system.

Biomechanics of human trabecular meshwork in healthy and glaucoma eyes via dynamic Schlemm's canal pressurization.

BACKGROUND AND OBJECTIVE: The trabecular meshwork (TM) consists of extracellular matrix (ECM) with embedded collagen and elastin fibers providing its mechanical support. TM stiffness is considerably higher in glaucoma eyes. Emerging data indicates that the TM moves dynamically with transient intraocular pressure (IOP) fluctuations, implying the viscoelastic mechanical behavior of the TM. However, little is known about TM viscoelastic behavior. We calculated the viscoelastic mechanical properties of the TM in n = 2 healthy and n = 2 glaucoma eyes. METHODS: A quadrant of the anterior segment was submerged in a saline bath, and a cannula connected to an adjustable saline reservoir was inserted into Schlemm's canal (SC). A spectral domain-OCT (SD-OCT) provided continuous cross-sectional B-scans of the TM/JCT/SC complex during pressure oscillation from 0 to 30 mmHg at two locations. The TM/JCT/SC complex boundaries were delineated to construct a 20-µm-thick volume finite element (FE) mesh. Pre-tensioned collagen and elastin fibrils were embedded in the model using a mesh-free penalty-based cable-in-solid algorithm. SC pressure was represented by a position- and time-dependent pressure boundary; floating boundary conditions were applied to the other cut edges of the model. An FE-optimization algorithm was used to adjust the ECM/fiber mechanical properties such that the TM/JCT/SC model and SD-OCT imaging data best matched over time. RESULTS: Significantly larger short- and long-time ECM shear moduli (p = 0.0032), and collagen (1.82x) and elastin (2.72x) fibril elastic moduli (p = 0.0001), were found in the TM of glaucoma eyes compared to healthy controls. CONCLUSIONS: These findings provide additional clarity on the mechanical property differences in healthy and glaucomatous outflow pathway under dynamic loading. Understanding the viscoelastic properties of the TM may serve as a new biomarker in early diagnosis of glaucoma.

The Effect of Intraocular Pressure Load Boundary on the Biomechanics of the Human Conventional Aqueous Outflow Pathway.

BACKGROUND: Aqueous humor outflow resistance in the trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and Schlemm's canal (SC) endothelium of the conventional outflow pathway actively contribute to intraocular pressure (IOP) regulation. Outflow resistance is actively affected by the dynamic outflow pressure gradient across the TM, JCT, and SC inner wall tissues. The resistance effect implies the presence of a fluid-structure interaction (FSI) coupling between the outflow tissues and the aqueous humor. However, the biomechanical interactions between viscoelastic outflow tissues and aqueous humor dynamics are largely unknown. METHODS: A 3D microstructural finite element (FE) model of a healthy human eye TM/JCT/SC complex was constructed with elastic and viscoelastic material properties for the bulk extracellular matrix and embedded elastic cable elements. The FE models were subjected to both idealized and a physiologic IOP load boundary using the FSI method. RESULTS: The elastic material model for both the idealized and physiologic IOP load boundary at equal IOPs showed similar stresses and strains in the outflow tissues as well as pressure in the aqueous humor. However, outflow tissues with viscoelastic material properties were sensitive to the IOP load rate, resulting in different mechanical and hydrodynamic responses in the tissues and aqueous humor. CONCLUSIONS: Transient IOP fluctuations may cause a relatively large IOP difference of ~20 mmHg in a very short time frame of ~0.1 s, resulting in a rate stiffening in the outflow tissues. Rate stiffening reduces strains and causes a rate-dependent pressure gradient across the outflow tissues. Thus, the results suggest it is necessary to use a viscoelastic material model in outflow tissues that includes the important role of IOP load rate.

Viscoelastic Biomechanical Properties of the Conventional Aqueous Outflow Pathway Tissues in Healthy and Glaucoma Human Eyes.

BACKGROUND: Although the tissues comprising the ocular conventional outflow pathway have shown strong viscoelastic mechanical response to aqueous humor pressure dynamics, the viscoelastic mechanical properties of the trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and Schlemm's canal (SC) inner wall are largely unknown. METHODS: A quadrant of the anterior segment from two human donor eyes at low- and high-flow (LF and HF) outflow regions was pressurized and imaged using optical coherence tomography (OCT). A finite element (FE) model of the TM, the adjacent JCT, and the SC inner wall was constructed and viscoelastic beam elements were distributed in the extracellular matrix (ECM) of the TM and JCT to represent anisotropic collagen. An inverse FE-optimization algorithm was used to calculate the viscoelastic properties of the ECM/beam elements such that the TM/JCT/SC model and OCT imaging data best matched over time. RESULTS: The ECM of the glaucoma tissues showed significantly larger time-dependent shear moduli compared to the heathy tissues. Significantly larger shear moduli were also observed in the LF regions of both the healthy and glaucoma eyes compared to the HF regions. CONCLUSIONS: The outflow tissues in both glaucoma eyes and HF regions are stiffer and less able to respond to dynamic IOP.

Biomechanics of the keratoconic cornea: Theory, segmentation, pressure distribution, and coupled FE-optimization algorithm.

Understanding of the corneal biomechanical properties is of high interest due to its potential application in the early diagnosis of keratoconus (KC). KC by itself is a non-inflammatory eye disorder causes corneal structural and/or compositional anomalies. The biomechanically weakened cornea is no longer able to preserve the normal shape of the cornea against the intraocular pressure (IOP) and gradually starts to bulge outward, invoking a conical shape and subsequent distorted vision. The most popular way to measure the in vivo corneal biomechanical properties is the CorVis-ST, which enables to analyze the dynamic response of the cornea under a temporary air puff pressure. However, the complications, such as the lack of knowledge on the accurate air-puff pressure distribution on the cornea's surface as a function of the distance from the apex of the cornea as well as the time, hinder us to have a reliable estimation of the cornea's mechanical parameters. This study aims to establish patient-specific geometries of the healthy and KC corneas and calculate the pressure distribution on the cornea's surface as a function of both the distance from the apex of the cornea and time, and thereafter, the viscoelastic mechanical properties of both the healthy and KC corneas using a coupled finite element (FE)-optimization algorithm. To do that, the dynamic deformation response of six healthy and six KC corneas were measured via CorVis-ST. The videos of the in vivo deformation of the corneas under the applied air puff pressure were segmented using our segmentation algorithm to determine the anterior and posterior curvatures of the corneas during the dynamic movement of the cornea. The FE model of the corneas were established using the segmented data and subjected to a negative (pre-stress), positive IOP, and air puff pressure while the floating boundary conditions were applied to the two ends of the corneas' FE models. The simulation results were imported into a loop of FE-optimization algorithm and analyzed until the deformation amplitude at the apex of the cornea reaches to its minimum difference compared to the clinical data by CorVis-ST. The results revealed that the pressure distributions found in the literature as a function of the distance from the apex of the cornea and time unable to provide satisfactory results. Therefore, the pressure distributions both as a function of the distance and time were optimized using our coupled FE-optimization algorithm and employed to estimate the viscoelastic properties of the healthy and KC corneas. The mean percentage error (MPE) of 8.45% and 10.79% were found for the healthy and KC corneas compared to the clinical data of CorVis-ST, respectively. The results also revealed a significantly higher short-time shear modulus for the KC (62.33 MPa) compared to the healthy (37.45 MPa) corneas while the long-time shear modulus of both the healthy and KC corneas were almost the same (4.01 vs. 3.91 MPa). The proposed algorithm is a noninvasive technique capable of accurately estimating the viscoelastic mechanical properties of the cornea, which can contribute to understand the mechanism of KC development and improve diagnosis and intervention in KC.

Modeling the biomechanics of the lamina cribrosa microstructure in the human eye.

Glaucoma is among the leading causes of blindness worldwide that is characterized by irreversible damage to the retinal ganglion cell axons in the lamina cribrosa (LC) region of the optic nerve head (ONH), most often associated with elevated intraocular pressure (IOP). The LC is a porous, connective tissue structure that provides mechanical support to the axons as they exit the eye and the biomechanics of the LC microstructure likely play a crucial role in protecting the axons passing through it. There is a limited knowledge of the IOP-driven biomechanics of the LC microstructure, primarily due to its small size and the difficulty with imaging the LC both in vitro and in vivo. We present finite element (FE) models of three human eye posterior poles that include the LC microstructure and interspersed neural tissues (NT) composed of retinal axons that are constructed directly from segmented, binary images of the LC. These models were used to estimate the stresses and strains in the LC and NT for an acute IOP elevation from 0 to 45 mmHg and compared with identical models except that the LC was represented as a homogenized continuum material with either homogeneous isotropic neo-Hookean properties or heterogeneous properties derived from local connective tissue volume fraction (CTVF) and predominant LC beam orientation. Stresses and strains in the LC and NT microstructure were investigated, and results were compared against those from the models wherein the LC was represented as a homogenized continuum. The regionalized volumetric average stresses and strains showed that the microstructural model yielded similar patterns to our prior approach using an LC continuum representation with mapped LC CTVF/anisotropy, but the microstructural modeling approach allows analysis of the stresses and strains in the LC and NT separately. As expected, the LC beams carried most of the IOP load in the microstructural models but exhibited less strain, while the encapsulated NT exhibited lower stresses and much higher strains. Results also revealed that the continuum models underestimate the maximum strains in the LC beams and NT by a factor of 2-3. Microstructural modeling should provide greater insight into the biomechanical factors driving damage to the axons (NT) and LC connective tissue remodeling that occur in glaucoma. The methods presented are ideal for modeling any structure with a complex microstructure composed of different materials, such as trabecular bone, lung, and tissue engineering scaffolds such as decellularized LC. Matlab code for mesh generation from a segmented image stack of the microstructure is included as Supplemental Material. STATEMENT OF SIGNIFICANCE: Glaucoma is among the leading causes of blindness worldwide that is characterized by axon damage in the lamina cribrosa (LC) region of the eye. We present a new approach for finite element modeling the entire eye-specific 3D LC microstructure and the interspersed neural tissues, incorporated into an eye-specific posterior eye model that provides appropriate boundary and loading conditions. Results are presented for three human donor eyes, showing that prior modeling approaches underestimate the stresses and strains in the laminar microstructure. We constructed models from image stacks of the segmented microstructure (Matlab code included) using an approach that is ideal for modeling any structure with a complex microstructure composed of different materials, such as trabecular bone, lung, and tissue engineering scaffolds.

Analysis of the effects of finite element type within a 3D biomechanical model of a human optic nerve head and posterior pole.

BACKGROUND AND OBJECTIVE: Biomechanical stresses and strains can be simulated in the optic nerve head (ONH) using the finite element (FE) method, and various element types have been used. This study aims to investigate the effects of element type on the resulting ONH stresses and strains. METHODS: A single eye-specific model was constructed using 3D delineations of anatomic surfaces in a high-resolution, fluorescent, 3D reconstruction of a human posterior eye, then meshed using our simple meshing algorithm at various densities using 4- and 10-noded tetrahedral elements, as well as 8- and 20-noded hexahedral elements. A mesh-free approach was used to assign heterogeneous, anisotropic, hyperelastic material properties to the lamina cribrosa, sclera and pia. The models were subjected to elevated IOP of 45 mmHg after pre-stressing from 0 to 10 mmHg, and solved in the open-source FE package Calculix; results were then interpreted in relation to computational time and simulation accuracy, using the quadratic hexahedral model as the reference standard. RESULTS: The 10-noded tetrahedral and 20R-noded hexahedral elements exhibited similar scleral canal and laminar deformations, as well as laminar and scleral stress and strain distributions; the quadratic tetrahedral models ran significantly faster than the quadratic hexahedral models. The linear tetrahedral and hexahedral elements were stiffer compared to the quadratic element types, yielding much lower stresses and strains in the lamina cribrosa. CONCLUSIONS: Prior studies have shown that 20-noded hexahedral elements yield the most accurate results in complex models. Results show that 10-noded tetrahedral elements yield very similar results to 20-noded hexahedral elements and so they can be used interchangeably, with significantly lower computational time. Linear element types did not yield acceptable results.

Finite element modeling of the periodontal ligament under a realistic kinetic loading of the jaw system.

PURPOSE: The stresses and deformations in the periodontal ligament (PDL) under the realistic kinetic loading of the jaw system, i.e., chewing, are difficult to be determined numerically as the mechanical properties of the PDL is variably present in different finite element (FE) models. This study was aimed to conduct a dynamic finite element (FE) simulation to investigate the role of the PDL (PDL) material models in the induced stresses and deformations using a simplified patient-specific FE model of a human jaw system. METHODS: To do that, a realistic kinetic loading of chewing was applied to the incisor point, contralateral, and ipsilateral condyles, through the experimentally proven trajectory approach. Three different material models, including the elasto-plastic, hyperelastic, and viscoelastic, were assigned to the PDL, and the resulted stresses of the tooth FE model were computed and compared. RESULTS: The results revealed the highest von Mises stress of 620.14 kPa and the lowest deformation of 0.16 mm in the PDL when using the hyperelastic model. The concentration of the stress in the elastoplastic and viscoelastic models was in the mid-root and apex of the PDL, while for the hyperelastic model, it was concentrated in the cervical margin. The highest deformation in the PDL regardless of the employed material model was located in the caudal direction of the tooth. The viscoelastic PDL absorbed the transmitted energy from the dentine and led to lower stress in the cancellous bone compared to the elastoplastic and hyperelastic material models. CONCLUSION: These results have implications not only for understanding the stresses and deformations in the PDL under chewing but also for providing comprehensive information for the medical and biomechanical experts in regard of the role of the material models being used to address the mechanical behavior of the PDL in other components of the tooth.

Mechanical measurement of the human cerebellum under compressive loading.

The cerebellum is responsible for controlling the posture and walking stability of the body. The cerebellum can subject to the traumatic injuries following by complicated clinical problems, i.e., the cerebellar pathologies. Application of the computational models can be helpful to figure out the injury mechanisms of the cerebellum, however, there is a lack of knowledge on the mechanical properties of the cerebellum under compressive loading. Therefore, this study aimed to perform an experimental study to measure the mechanical properties of 17 male individuals' cerebellum under the series of compressive loadings. The resulted stress-strain data of the cerebellum revealed the elastic modulus and maximum/failure stress of 13.48 ± 2.65 (Mean ± SD) and 19.65 ± 3.89 kPa, respectively. The findings of this study have implications not only for understanding the mechanical properties of the human cerebellum tissue under the compressive loading, but also for providing a raw data for the doctors and biomechanical experts as the mechanical threshold of the cerebellum as well as computational modelling of the traumatic brain injuries.

Prediction of human gait trajectories during the SSP using a neuromusculoskeletal modeling: A challenge for parametric optimization.

The parametric optimization techniques have been widely employed to predict human gait trajectories; however, their applications to reveal the other aspects of gait are questionable. The aim of this study is to investigate whether or not the gait prediction model is able to justify the movement trajectories for the higher average velocities. A planar, seven-segment model with sixteen muscle groups was used to represent human neuro-musculoskeletal dynamics. At first, the joint angles, ground reaction forces (GRFs) and muscle activations were predicted and validated for normal average velocity (1.55 m/s) in the single support phase (SSP) by minimizing energy expenditure, which is subject to the non-linear constraints of the gait. The unconstrained system dynamics of extended inverse dynamics (USDEID) approach was used to estimate muscle activations. Then by scaling time and applying the same procedure, the movement trajectories were predicted for higher average velocities (from 2.07 m/s to 4.07 m/s) and compared to the pattern of movement with fast walking speed. The comparison indicated a high level of compatibility between the experimental and predicted results, except for the vertical position of the center of gravity (COG). It was concluded that the gait prediction model can be effectively used to predict gait trajectories for higher average velocities.

Biomechanics of the Healthy and Keratoconic Corneas: A Combination of the Clinical Data, Finite Element Analysis, and Artificial Neural Network.

BACKGROUND: Keratoconus is recognized by asymmetrical thinning and bulging of the cornea, resulting in distortion in the surface of the cornea. Keratoconus also alters the biomechanical properties of the cornea, which can be an indicator of the healthy and keratoconus eyes. This study was aimed at employing a combination of clinical data, finite element method (FEM), and artificial neural network (ANN) to establish a novel biomechanical- based diagnostic method for the keratoconus eyes. METHODS: To do that, the clinical-biomechanical parameters of 40 healthy and 40 keratoconus eyes were obtained via the Pentacam and non-contact tonometer (Corvis ST, Oculus Optikgeräte, Wetzlar, Germany) devices. Intraocular pressure (IOP) was measured using a Goldmann applanation tonometer as well as Corvis. According to the geometry of the cornea, the FE model of each cornea was made and the same boundary and loading conditions were applied not only to confirm the FE model in terms of the biomechanical parameters but also to calculate the amount of von Mises stress in the apex of the cornea. The clinical-biomechanical data of the Corvis along with the von Mises stresses were then incorporated into the ANN algorithm to distinguish the healthy and keratoconus corneas on a basis of the resulted von Mises stresses. The proposed programming code, according to the input data from the Corvis, enabled to predict whether the cornea is keratoconus or not. Finally, to verify the results of the proposed method, 155 individuals were examined. RESULTS: The clinical and biomechanical results of the Corvis revealed that the healthy corneas have a higher thickness compared to the keratoconus ones. No significant differences were observed among the IOPs, 1st applanation length, and pick distance in the highest concavity. The 2nd applanation length and radius in the highest concavity of the healthy cornea were higher than the keratoconus ones. Conversely, the 1st and 2nd applanation velocities and deformation amplitudes of the keratoconus corneas were higher than the healthy ones. The FE results also showed higher stresses for the healthy corneas compared to the keratoconus ones. The ANN was also well verified since it demonstrated more than 95.5% accuracy on diagnosing the keratoconus eyes. CONCLUSION: These findings have implications not only for identifying the keratoconus corneas as an important clinical and surgical tool for eye care professionals but also for providing both a quantitative and an accurate approach to the problem of understanding the biomechanical nature of keratoconus.

A combination of experimental measurement, constitutive damage model, and diffusion tensor imaging to characterize the mechanical properties of the human brain.

Understanding the mechanical properties of the human brain is deemed important as it may subject to various types of complex loadings during the Traumatic Brain Injury (TBI). Although many studies so far have been conducted to quantify the mechanical properties of the brain, there is a paucity of knowledge on the mechanical properties of the human brain tissue and the damage of its axon fibers under the various types of complex loadings during the Traumatic Brain Injury (TBI). Although many studies so far have been conducted to quantify the mechanical properties of the brain, there is a paucity of knowledge on the mechanical properties of the human brain tissue and the damage of its axon fibers under the frontal lobe of the human brain. The constrained nonlinear minimization method was employed to identify the brain coefficients according to the axial and transversal compressive data. The pseudo-elastic damage model data was also well compared with that of the experimental data and it not only up to the primary loading but also the discontinuous softening could well address the mechanical behavior of the brain tissue.