Viscoelastic mechanical measurement of the healthy and atherosclerotic human coronary arteries using DIC technique

Finite Element Analysis of Cervical Spine Kinematic Response during Ejection Utilising a Hill-Type Dynamic Muscle Model

To determine the impact of active muscle on the dynamic response of a pilot's neck during simulated emergency ejection, a detailed three-dimensional (3D) cervical spine (C0-T1) finite element (FE) model integrated with active muscles was constructed. Based on the Hill-type model characterising the muscle force activation mechanics, 13 major neck muscles were modelled. The active force generated by each muscle was simulated as functions of (i) active state (Na), (ii) velocity (Fv(v)), and (iii) length (FL(L)). An acceleration-time profile with an initial acceleration rate of 125 G·s-1 in the 0-80 ms period, reaching peak acceleration of 10 G, then kept constant for a further 70 ms, was applied. The rotational angles of each cervical segment under these ejection conditions were compared with those without muscles and with passive muscles derived from the previous study. Similar trends of segmental rotation were observed with S- and C-curvature of the cervical spine in the 150 ms span analysed. With active muscles, the flexion motion of the C0-C2 segments exhibited higher magnitudes of rotation compared to those without muscle and passive muscle models. The flexion motion increased rapidly and peaked at about 95-105 ms, then decreased rapidly to a lower magnitude. Lower C2-T1 segments exhibited less variation in flexion and extension motions. Overall, during emergency ejections, active muscle activities effectively reduce the variability in rotational angles across cervical segments, except C0-C2 segments in the 60-120 ms period. The role of the active state dynamics of the muscles was crucial to the magnitude of the muscle forces demonstrated. This indicates that it is crucial for pilots to consciously contract their muscles before ejection to prevent cervical spine injuries.

Monitoring the Progression of Early Atherosclerosis Using a Fluorescence Nanoprobe for the Detection and Imaging of Phosphorylation and Glucose Levels

Monitoring the early stage of atherosclerosis (AS) without plaque formation is of great significance. Herein, we developed a metal organic framework (MOF)-based fluorescence nanoprobe to analyze the progression of AS by assessing the levels of protein phosphorylation and glucose in blood and tissue. The probe was prepared by post-modification of the MOF with iodine (I3 - )-rhodamine B (RhB) associate, which realizes the specific recognition of target object through the metal joint ZrIV and I3 - -RhB, respectively. We investigated different stages of target object changes in the early non-plaque stage of AS in blood. It was found that the levels of phosphate and glucose in the blood were higher than those of the normal mice. The results of two-photon images showed that early AS mice had higher levels of protein phosphorylation and glucose than that of the normal mice. The present study provides a suitable fluorescence tool for further revealing the pathogenesis and progression of AS.

A position- and time-dependent pressure profile to model viscoelastic mechanical behavior of the brain tissue due to tumor growth

This study proposed a computational framework to calculate the resultant position- and time-dependent pressure profile on the brain tissue due to tumor growth. A finite element (FE) patch of the brain tissue was constructed and an inverse dynamic FE-optimization algorithm was used to calculate its viscoelastic mechanical properties under compressive uniaxial loading. Two patient-specific post-tumor resection FE models were input to the FE-optimization algorithm to calculate the optimized 3rd-order position-dependent and normal distribution time-dependent pressure profile parameters. The optimized viscoelastic material properties, the most suitable simulation time, and the optimized 3rd-order position- and -time-dependent pressure profiles were calculated.

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 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.

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.

Toward Elucidating the Physiological Impacts of Residual Stresses in the Colorectum

Irritable bowel syndrome afflicts 10-20% of the global population, causing visceral pain with increased sensitivity to colorectal distension and normal bowel movements. Understanding and predicting these biomechanics will further advance our understanding of visceral pain and complement the existing literature on visceral neurophysiology. We recently performed a series of experiments at three longitudinal segments (colonic, intermediate, and rectal) of the distal 30 mm of colorectums of mice. We also established and fitted constitutive models addressing mechanical heterogeneity in both the through-thickness and longitudinal directions of the colorectum. Afferent nerve endings, strategically located within the submucosa, are likely nociceptors that detect concentrations of mechanical stresses to evoke the perception of pain from the viscera. In this study, we aim to: (1) establish and validate a method for incorporating residual stresses into models of colorectums, (2) predict the effects of residual stresses on the intratissue mechanics within the colorectum, and (3) establish intratissue distributions of stretches and stresses within the colorectum in vivo. To these ends we developed two-layered, composite finite element models of the colorectum based on our experimental evidence and validated our approaches against independent experimental data. We included layer- and segment-specific residual stretches/stresses in our simulations via the prestrain algorithm built into the finite element software febio. Our models and modeling approaches allow researchers to predict both organ and intratissue biomechanics of the colorectum and may facilitate better understanding of the underlying mechanical mechanisms of visceral pain.

In Vivo Based Fluid-Structure Interaction Biomechanics of the Left Anterior Descending Coronary Artery

A fluid-structure interaction-based biomechanical model of the entire left anterior descending coronary artery is developed from in vivo imaging via the finite element method in this paper. Included in this investigation is ventricle contraction, three-dimensional motion, all angiographically visible side branches, hyper/viscoelastic artery layers, non-Newtonian and pulsatile blood flow, and the out-of-phase nature of blood velocity and pressure. The fluid-structure interaction model is based on in vivo angiography of an elite athlete's entire left anterior descending coronary artery where the influence of including all alternating side branches and the dynamical contraction of the ventricle is investigated for the first time. Results show the omission of side branches result in a 350% increase in peak wall shear stress and a 54% decrease in von Mises stress. Peak von Mises stress is underestimated by up to 80% when excluding ventricle contraction and further alterations in oscillatory shear indices are seen, which provide an indication of flow reversal and has been linked to atherosclerosis localization. Animations of key results are also provided within a video abstract. We anticipate that this model and results can be used as a basis for our understanding of the interaction between coronary and myocardium biomechanics. It is hoped that further investigations could include the passive and active components of the myocardium to further replicate in vivo mechanics and lead to an understanding of the influence of cardiac abnormalities, such as arrythmia, on coronary biomechanical responses.

Measurement of the mechanical properties of the human gallbladder

Gallbladder is a small organ of the body which is located in the right side of the liver. It is responsible of storing the bile and releasing it to the intestine. The gallbladder can subject to the mechanical deformation/loading as a result of the cholecystitis, cholesterolosis of the gallbladder, etc. However, so far the mechanical properties of the human gallbladder have not been measured. This study was aimed at conducting an experimental study to measure the mechanical properties of the human gallbladder under the axial and transversal tensile loadings. To do that, the gallbladder tissue of 16 male individuals was excised during the autopsy and subjected to a series of axial and transversal loadings under the strain rate of 5 mm/min. The amount of elastic modulus as well as the maximum/failure stress of the tissues were calculated via the resulted stress-strain diagrams and reported. The results revealed that the axial and transversal elastic modulus were 641.20 ± 28.12 (mean ± SD) and 255 ± 24.55 kPa, respectively. The amount of maximum stresses was also 1240 ± 99.94 and 348 ± 66.75 kPa under the axial and transversal loadings, respectively. The results revealed a significantly higher axial stiffness (p < .05, post hoc Scheffe method) compared to the transversal one. These findings have implications not only for understanding the axial and transversal mechanical properties of the human gallbladder tissue, but also for providing a diagnosis tool for the doctors to have a suitable threshold value of the healthy gallbladder tissue.

Optimizing through computational modeling to reduce dogboning of functionally graded coronary stent material

Coronary artery disease is the leading cause of death among the men and women. One of the most suitable treatments for this problem is balloon angioplasty with stenting. Functionally graded material (FGM) stents have shown suitable mechanical behavior in simulations. While their deformation was superior to uniform materials, the study was aimed at finding the most suitable configuration to reach the optimum performance. A combination of finite element method (FEM) and optimization algorithm have been used to fulfil this objective. To do that, three different conditions have been investigated in a Palmaz-Schatz geometry, where in the first and second ones the stent was a combination of steel and CoCr alloy (L605), and the third condition was a combination of CoCr alloy (L605) and CoCr alloy (F562). In the first and third conditions, dogboning was the objective function, but in the second condition a non-uniform deformation indicator was chosen as the objective function. In all three conditions the heterogeneous index was the control variable. The stent in the third condition showed a poor performance. While in the steel/CoCr alloy (L605) stents the heterogeneous index of 0.374 showed the lowest maximum dogboning, the heterogeneous index of 5 had more uniform deformation. Overall due to the lower dogboning of the steel/CoCr alloy (L605) stent with heterogeneous index of 0.374, this stent is recommended as the optimum stent in this geometrical configuration.