A transversely isotropic viscoelastic constitutive equation for brainstem undergoing finite deformation

The objective of this study was to define the constitutive response of brainstem undergoing finite shear deformation. Brainstem was characterized as a transversely isotropic viscoelastic material and the material model was formulated for numerical implementation. Model parameters were fit to shear data obtained in porcine brainstem specimens undergoing finite shear deformation in three directions: parallel, perpendicular, and cross sectional to axonal fiber orientation and determined using a combined approach of finite element analysis (FEA) and a genetic algorithm (GA) optimizing method. The average initial shear modulus of brainstem matrix of 4-week old pigs was 12.7 Pa, and therefore the brainstem offers little resistance to large shear deformations in the parallel or perpendicular directions, due to the dominant contribution of the matrix in these directions. The fiber reinforcement stiffness was 121.2 Pa, indicating that brainstem is anisotropic and that axonal fibers have an important role in the cross-sectional direction. The first two leading relative shear relaxation moduli were 0.8973 and 0.0741, respectively, with corresponding characteristic times of 0.0047 and 1.4538 s, respectively, implying rapid relaxation of shear stresses. The developed material model and parameter estimation technique are likely to find broad applications in neural and orthopaedic tissues.     

Engineering controllable anisotropy in electrospun biodegradable nanofibrous scaffolds for musculoskeletal tissue engineering

Many musculoskeletal tissues exhibit significant anisotropic mechanical properties reflective of a highly oriented underlying extracellular matrix. For tissue engineering, recreating this organization of the native tissue remains a challenge. To address this issue, this study explored the fabrication of biodegradable nanofibrous scaffolds composed of aligned fibers via electrospinning onto a rotating target, and characterized their mechanical anisotropy as a function of the production parameters. The characterization showed that nanofiber organization was dependent on the rotation speed of the target; randomly oriented fibers (33% fiber alignment) were produced on a stationary shaft, whereas highly oriented fibers (94% fiber alignment) were produced when rotation speed was increased to 9.3m/s. Non-aligned scaffolds had an isotropic tensile modulus of 2.1+/-0.4MPa, compared to highly anisotropic scaffolds whose modulus was 11.6+/-3.1MPa in the presumed fiber direction, suggesting that fiber alignment has a profound effect on the mechanical properties of scaffolds. Mechanical anisotropy was most pronounced at higher rotation speeds, with a greater than 33-fold enhancement of the Young's modulus in the fiber direction compared to perpendicular to the fiber direction when the rotation speed reached 8m/s. In cell culture, both the organization of actin filaments of human mesenchymal stem cells and the cellular alignment of meniscal fibroblasts were dictated by the prevailing nanofiber orientation. This study demonstrates that controllable and anisotropic mechanical properties of nanofibrous scaffolds can be achieved by dictating nanofiber organization through intelligent scaffold design.     

Improved prediction of the collagen fiber architecture in the aortic heart valve

Living tissues show an adaptive response to mechanical loading by changing their internal structure and morphology. Understanding this response is essential for successful tissue engineering of load-bearing structures, such as the aortic valve. In this study, mechanically induced remodeling of the collagen architecture in the aortic valve was investigated. It was hypothesized that, in uniaxially loaded regions, the fibers aligned with the tensile principal stretch direction. For biaxial loading conditions, on the other hand, it was assumed that the collagen fibers aligned with directions situated between the principal stretch directions. This hypothesis has already been applied successfully to study collagen remodeling in arteries. The predicted fiber architecture represented a branching network and resembled the macroscopically visible collagen bundles in the native leaflet. In addition, the complex biaxial mechanical behavior of the native valve could be simulated qualitatively with the predicted fiber directions. The results of the present model might be used to gain further insight into the response of tissue engineered constructs during mechanical conditioning.     

Computational model predicts cell orientation in response to a range of mechanical stimuli

To build anisotropic, mechanically functioning tissue, it is essential to understand how cells orient in response to mechanical stimuli. Therefore, a computational model was developed which predicts cell orientation, based on the actin stress fiber distribution inside the cell. In the model, the stress fiber distribution evolves dynamically according to the following: (1) Stress fibers contain polymerized actin. The total amount of depolymerized plus polymerized actin is constant. (2) Stress fibers apply tension to their environment. This active tension is maximal when strain rate and absolute strain are zero and reduces with increasing shortening rate and absolute strain. (3) A high active fiber stress in a direction leads to a large amount of fibers in this direction. (4) The cell is attached to a substrate; all fiber stresses are homogenized into a total cell stress, which is in equilibrium with substrate stress. This model predicts that on a substrate of anisotropic stiffness, fibers align in the stiffest direction. Under cyclic strain when the cellular environment is so stiff that no compaction occurs (1 MPa), the model predicts strain avoidance, which is more pronounced with increasing strain frequency or amplitude. Under cyclic strain when the cellular environment is so soft that cells can compact it (10 kPa), the model predicts a preference for the cyclically strained compared to the compacting direction. These model predictions all agree with experimental evidence. For the first time, a computational model predicts cell orientation in response to this range of mechanical stimuli using a single set of parameters.     

The influence of anisotropy on brain injury prediction

Traumatic Brain Injury (TBI) occurs when a mechanical insult produces damage to the brain and disrupts its normal function. Numerical head models are often used as tools to analyze TBIs and to measure injury based on mechanical parameters. However, the reliability of such models depends on the incorporation of an appropriate level of structural detail and accurate representation of the material behavior. Since recent studies have shown that several brain regions are characterized by a marked anisotropy, constitutive equations should account for the orientation-dependence within the brain. Nevertheless, in most of the current models brain tissue is considered as completely isotropic. To study the influence of the anisotropy on the mechanical response of the brain, a head model that incorporates the orientation of neural fibers is used and compared with a fully isotropic model. A simulation of a concussive impact based on a sport accident illustrates that significantly lowered strains in the axonal direction as well as increased maximum principal strains are detected for anisotropic regions of the brain. Thus, the orientation-dependence strongly affects the response of the brain tissue. When anisotropy of the whole brain is taken into account, deformation spreads out and white matter is particularly affected. The introduction of local axonal orientations and fiber distribution into the material model is crucial to reliably address the strains occurring during an impact and should be considered in numerical head models for potentially more accurate predictions of brain injury.     

Remodelling of the angular collagen fiber distribution in cardiovascular tissues

Understanding collagen fiber remodelling is desired to optimize the mechanical conditioning protocols in tissue-engineering of load-bearing cardiovascular structures. Mathematical models offer strong possibilities to gain insight into the mechanisms and mechanical stimuli involved in these remodelling processes. In this study, a framework is proposed to investigate remodelling of angular collagen fiber distribution in cardiovascular tissues. A structurally based model for collagenous cardiovascular tissues is extended with remodelling laws for the collagen architecture, and the model is subsequently applied to the arterial wall and aortic valve. For the arterial wall, the model predicts the presence of two helically arranged families of collagen fibers. A branching, diverging hammock-type fiber architecture is predicted for the aortic valve. It is expected that the proposed model may be of great potential for the design of improved tissue engineering protocols and may give further insight into the pathophysiology of cardiovascular diseases.     

Multi-scale mechanics of traumatic brain injury: predicting axonal strains from head loads

The length scales involved in the development of diffuse axonal injury typically range from the head level (i.e., mechanical loading) to the cellular level. The parts of the brain that are vulnerable to this type of injury are mainly the brainstem and the corpus callosum, which are regions with highly anisotropically oriented axons. Within these parts, discrete axonal injuries occur mainly where the axons have to deviate from their main course due to the presence of an inclusion. The aim of this study is to predict axonal strains as a result of a mechanical load at the macroscopic head level. For this, a multi-scale finite element approach is adopted, in which a macro-level head model and a micro-level critical volume element are coupled. The results show that the axonal strains cannot be trivially correlated to the tissue strain without taking into account the axonal orientations, which indicates that the heterogeneities at the cellular level play an important role in brain injury and reliable predictions thereof. In addition to the multi-scale approach, it is shown that a novel anisotropic equivalent strain measure can be used to assess these micro-scale effects from head-level simulations only.     

A multi-modular tensegrity model of an actin stress fiber

Stress fibers are contractile bundles in the cytoskeleton that stabilize cell structure by exerting traction forces on the extracellular matrix. Individual stress fibers are molecular bundles composed of parallel actin and myosin filaments linked by various actin-binding proteins, which are organized end-on-end in a sarcomere-like pattern within an elongated three-dimensional network. While measurements of single stress fibers in living cells show that they behave like tensed viscoelastic fibers, precisely how this mechanical behavior arises from this complex supramolecular arrangement of protein components remains unclear. Here we show that computationally modeling a stress fiber as a multi-modular tensegrity network can predict several key behaviors of stress fibers measured in living cells, including viscoelastic retraction, fiber splaying after severing, non-uniform contraction, and elliptical strain of a puncture wound within the fiber. The tensegrity model can also explain how they simultaneously experience passive tension and generate active contraction forces; in contrast, a tensed cable net model predicts some, but not all, of these properties. Thus, tensegrity models may provide a useful link between molecular and cellular scale mechanical behaviors and represent a new handle on multi-scale modeling of living materials.     

Scleral anisotropy and its effects on the mechanical response of the optic nerve head

This paper presents a computational modeling study of the effects of the collagen fiber structure on the mechanical response of the sclera and the adjacent optic nerve head (ONH). A specimen-specific inverse finite element method was developed to determine the material properties of two human sclera subjected to full-field inflation experiments. A distributed fiber model was applied to describe the anisotropic elastic behavior of the sclera. The model directly incorporated wide-angle X-ray scattering measurements of the anisotropic collagen structure. The converged solution of the inverse method was used in micromechanical studies of the mechanical anisotropy of the sclera at different scales. The effects of the scleral collagen fiber structure on the ONH deformation were evaluated by progressively filtering out local anisotropic features. It was found that the majority of the midposterior sclera could be described as isotropic without significantly affecting the mechanical response of the tissues of the ONH. In contrast, removing local anisotropic features in the peripapillary sclera produced significant changes in scleral canal expansion and lamina cribrosa deformation. Local variations in the collagen structure of the peripapillary sclera significantly influenced the mechanical response of the ONH.     

A recruitment-based rheological model for mechanical behavior of soft tissues

When lung tissue is subjected to finite deformations, phenomena appear that can only be described using nonlinear models. This paper considers the lung as a material composed of two elements, a continuous phase that acts uninterruptedly and a second phase composed of fiber elements that are recruited progressively into the mechanical process. Each individual fiber participates in the mechanical response of the set only when the deformation is above a certain value. A nine-parameter model was designed adopting standard viscoelastic elements both for the matrix and for each of the fibers. The mechanical behavior of the lung can be reproduced by a fitting process with standard numerical procedures in both dynamic-mechanical measurements and stress relaxation processes. Mechanical stress relaxation tests and dynamic-mechanical measurements have been carried out on subpleural parenchymal strips from rat lung. The model permits the reproduction of lung behavior in both types of measurements. The results show a recruitment ratio that decreases with deformation and the nonparticipation of the parallel matrix fraction in the lung's mechanical response so that a uniaxial transmission of force in the lung occurs via the recruited elements and the matrix series.     

Dynamics of mechanical signal transmission through prestressed stress fibers

Transmission of mechanical stimuli through the actin cytoskeleton has been proposed as a mechanism for rapid long-distance mechanotransduction in cells; however, a quantitative understanding of the dynamics of this transmission and the physical factors governing it remains lacking. Two key features of the actin cytoskeleton are its viscoelastic nature and the presence of prestress due to actomyosin motor activity. We develop a model of mechanical signal transmission through prestressed viscoelastic actin stress fibers that directly connect the cell surface to the nucleus. The analysis considers both temporally stationary and oscillatory mechanical signals and accounts for cytosolic drag on the stress fibers. To elucidate the physical parameters that govern mechanical signal transmission, we initially focus on the highly simplified case of a single stress fiber. The results demonstrate that the dynamics of mechanical signal transmission depend on whether the applied force leads to transverse or axial motion of the stress fiber. For transverse motion, mechanical signal transmission is dominated by prestress while fiber elasticity has a negligible effect. Conversely, signal transmission for axial motion is mediated uniquely by elasticity due to the absence of a prestress restoring force. Mechanical signal transmission is significantly delayed by stress fiber material viscosity, while cytosolic damping becomes important only for longer stress fibers. Only transverse motion yields the rapid and long-distance mechanical signal transmission dynamics observed experimentally. For simple networks of stress fibers, mechanical signals are transmitted rapidly to the nucleus when the fibers are oriented largely orthogonal to the applied force, whereas the presence of fibers parallel to the applied force slows down mechanical signal transmission significantly. The present results suggest that cytoskeletal prestress mediates rapid mechanical signal transmission and allows temporally oscillatory signals in the physiological frequency range to travel a long distance without significant decay due to material viscosity and/or cytosolic drag.     

A micromechanical procedure for modelling the anisotropic mechanical properties of brain white matter

This paper proposes a micromechanics algorithm utilising the finite element method (FEM) for the analysis of heterogeneous matter. The characterisation procedure takes the material properties of the constituents, axons and extracellular matrix (ECM) as input data. The material properties of both the axons and the matrix are assumed to have linear viscoelastic behaviour with a perfect bonding between them. The results of the modelling have been validated with experimental data with material white input from brainstem by considering the morphology of brainstem in which most axons are oriented in longitudinal direction in the form of a uniaxial fibrous composite material. The method is then employed to examine the undulations of axons within different subregions of white matter and to study the impact due to axon/matrix volume fractions. For such purposes, different unit cells composed of wavy geometries and with various volume factions have been exposed to the six possible loading scenarios. The results will clearly demonstrate the undulation and axon volume fraction impacts. In this respect, undulation affects the material stiffness heavily in the axon longitudinal direction, whereas the axons' volume fraction has a much greater impact on the mechanical properties of the white matter in general. Also the results show that the created stresses and strains in the axons and matrix under loading will be impacted by undulation change. With increase in undulation the matrix suffers higher stresses when subjected to tension, whereas axons suffer higher stresses in shear. The axons always exhibit higher stresses whereas the matrix exhibits higher strains. The evaluated time-dependent local stress and strain concentrations within a repeating unit cell of the material model are indicative of the mechanical behaviour of the white tissue under different loading scenarios.     

A microstructurally based continuum model of cartilage viscoelasticity and permeability incorporating measured statistical fiber orientations

The remarkable mechanical properties of cartilage derive from an interplay of isotropically distributed, densely packed and negatively charged proteoglycans; a highly anisotropic and inhomogeneously oriented fiber network of collagens; and an interstitial electrolytic fluid. We propose a new 3D finite strain constitutive model capable of simultaneously addressing both solid (reinforcement) and fluid (permeability) dependence of the tissue’s mechanical response on the patient-specific collagen fiber network. To represent fiber reinforcement, we integrate the strain energies of single collagen fibers—weighted by an orientation distribution function (ODF) defined over a unit sphere—over the distributed fiber orientations in 3D. We define the anisotropic intrinsic permeability of the tissue with a structure tensor based again on the integration of the local ODF over all spatial fiber orientations. By design, our modeling formulation accepts structural data on patient-specific collagen fiber networks as determined via diffusion tensor MRI. We implement our new model in 3D large strain finite elements and study the distributions of interstitial fluid pressure, fluid pressure load support and shear stress within a cartilage sample under indentation. Results show that the fiber network dramatically increases interstitial fluid pressure and focuses it near the surface. Inhomogeneity in the tissue’s composition also increases fluid pressure and reduces shear stress in the solid. Finally, a biphasic neo-Hookean material model, as is available in commercial finite element codes, does not capture important features of the intra-tissue response, e.g., distributions of interstitial fluid pressure and principal shear stress.     

Toward "strong" green nanocomposites: polyvinyl alcohol reinforced with extremely oriented cellulose whiskers

To exploit the maximum potential of cellulose whiskers (CWs), we report here for the first time the successful fabrication of nanocomposites reinforced with highly oriented CWs in a polymer matrix. The nanocomposites were prepared using polyvinyl alcohol (PVA) and a colloidal suspension of cotton-derived CWs. The macroscopically homogeneous PVA-CW suspensions were extruded into cold methanol to form gel fibers followed by a hot drawing. Compared to the neat PVA fiber, the as-spun fiber containing a small amount of CWs (5 wt % of solid PVA) showed higher drawability, leading to an extremely high orientation of CWs with the matrix PVA. The stress-transfer mechanism, a prime determining factor for high mechanical properties of nanocomposites, was studied by X-ray diffraction. The stress on the incorporated CWs was monitored by applying an in situ nondestructive load to the composite fibers. The applied stress to the whole sample was found to be effectively transferred to the CWs inside the composites, suggesting strong interfacial bonding between the filler and the matrix. Effective stress transfer to the oriented whiskers resulted in outstanding enhancement in mechanical properties of the nanocomposites.     

A methodology for developing anisotropic AAA phantoms via additive manufacturing

An Abdominal Aortic Aneurysm (AAA) is a permanent focal dilatation of the abdominal aorta at least 1.5 times its normal diameter. The criterion of maximum diameter is still used in clinical practice, although numerical studies have demonstrated the importance of biomechanical factors for rupture risk assessment. AAA phantoms could be used for experimental validation of the numerical studies and for pre-intervention testing of endovascular grafts. We have applied multi-material 3D printing technology to manufacture idealized AAA phantoms with anisotropic mechanical behavior. Different composites were fabricated and the phantom specimens were characterized by biaxial tensile tests while using a constitutive model to fit the experimental data. One composite was chosen to manufacture the phantom based on having the same mechanical properties as those reported in the literature for human AAA tissue; the strain energy and anisotropic index were compared to make this choice. The materials for the matrix and fibers of the selected composite are, respectively, the digital materials FLX9940 and FLX9960 developed by Stratasys. The fiber proportion for the composite is equal to 0.15. The differences between the composite behavior and the AAA tissue are small, with a small difference in the strain energy (0.4%) and a maximum difference of 12.4% in the peak Green strain ratio. This work represents a step forward in the application of 3D printing technology for the manufacturing of AAA phantoms with anisotropic mechanical behavior.     

A linear material model for fiber-induced anisotropy of the anulus fibrosus

The anulus fibrosus (AF) is a lamellar, fibrocartilaginous component of the intervertebral disc, which exhibits highly anisotropic behaviors in tension. These behaviors arise from the material's unique collagen structure. We have investigated the use of a linear, fiber-induced anisotropic model for the AF using a quadratic strain energy density formulation with an explicit representation of the collagen fiber populations. We have proposed a representative set of intrinsic material properties using independent datasets of the AF from the literature and appropriate thermodynamic constraints. The model was validated by comparing predictions with previous experimental data for AF behavior and its dependence on fiber angle. The model predicts that compressible effects may exist for the AF, and suggests that physical effects of the equivalent "matrix," "fiber," "fiber-matrix," and "fiber-fiber," interactions may be important contributors to the mechanical behavior of the AF.     

On the anisotropy of the canine diaphragmatic central tendon

We studied the mechanical and anatomical anisotropy of the canine diaphragmatic central tendon (CT). Dumb-bell-shaped strips with effective dimensions of 10 x 2 mm (length x width) were cut from different regions of the canine diaphragmatic CT in two different orientations relative to the direction of neighboring muscle fibers. Specimens sampled with their long axial dimension oriented parallel to the neighboring muscle fibers were named Group-1 and those sampled with an orientation perpendicular to the neighboring muscle fibers were named Group-2. Results from one-dimensional stress-strain and tensile failure strength tests revealed that the CT is a nonlinear, inelastic, and anisotropic material. Group-1 specimens were found to have a higher stiffness, higher failure strength and higher strain energy density at failure than Group-2 specimens. Polarized microscopy showed that multiple sheets of collagen fiber bundles formed an orthogonal network in the tendon. Collagen fiber bundles along Group-1 direction formed parallel trajectory lines connecting the neighboring costal and crural muscles; bundles along Group-2 direction were observed to orient 90 degrees away. At the central apex region of the CT, collagen bundles of Group-1 formed a fan-like trajectory pattern. This collagen network architecture was compared favorably to the trajectories of an approximated principal stress field in the CT due to simulated contractile forces from its adjacent costal and crural muscles. These combined results suggest a structure-function relationship for the anatomical and mechanical anisotropy in the canine diaphragmatic CT.     

Effect of the anisotropic permeability in the frequency dependent properties of the superficial layer of articular cartilage

Abstract

Articular cartilage is a tissue of fundamental importance for the mechanics of joints, since it provides a smooth and lubricated surface for the proper transfer of loads. From a mechanical point of view, this tissue is an anisotropic poroviscoelastic material: its characteristics at the macroscopic level depend on the complex microscopic architecture. With the ability to probe the local microscopic features, dynamic nanoindentation test is a powerful tool to investigate cartilage mechanics. In this work we focus on a length scale where the time dependent behaviour is regulated by poroelasticity more than viscoelasticity and we aim to understand the effect of the anisotropic permeability on the mechanics of the superficial layer of the articular cartilage. In a previous work, a finite element model for the dynamic nanoindentation test has been presented. In this work, we improve the model by considering the presence of an anisotropic permeability tensor that depends on the collagen fibers distribution. Our sensitivity analysis highlights that the permeability decreases with increasing indentation, thus making the tissue stiffer than the case of isotropic permeability, when solicited at the same frequency. With this improved model, a revised identification of the mechanical and physical parameters for articular cartilage is provided. To this purpose the model was used to simulate experimental data from tests performed on bovine tissue, giving a better estimation of the anisotropy in the elastic properties. A relation between the identified macroscopic anisotropic permeability properties and the microscopic rearrangement of the fiber/matrix structure during indentation is also provided.     

Search for critical loading condition of the spine--a meta analysis of a nonlinear viscoelastic finite element model

The relative vulnerability of spinal motion segments to different loading combinations remains unknown. The meta-analysis described here using the results of a validated L2-L3 nonlinear viscoelastic finite element model was designed to investigate the critical loading and its effect on the internal mechanics of the human lumbar spine. A Box-Behnken experimental design was used to design the magnitude of seven independent variables associated with loads, rotations and velocity of motion. Subsequently, an optimization method was used to find the primary and secondary variables that influence spine mechanical output related to facet forces, disc pressure, ligament forces, annulus matrix compressive/shear stresses and anulus fibers strain. The mechanical responses with respect to the two most-relevant variables were then regressed linearly using the response surface quadratic model. Axial force and sagittal rotation were identified as the most-relevant variables for mechanical responses. The procedure developed can be used to find the critical loading for finite element models with multi input variables. The derived meta-models can be used to predict the risk associated with various loading parameters and in setting safer load limits.