A structural constitutive model for smooth muscle contraction in biological tissues

A structural constitutive model that characterizes the active and passive responses of biological tissues with smooth muscle cells (SMCs) is proposed. The model is formulated under the assumption that the contractile units in SMCs and the connected collagen fibers are the active tissue component, while the collagen fibers not connected to the SMCs are the passive tissue component. An evolution law describing the deformation of the active tissue component over time is developed based on the sliding filament theory. In this evolution law the contraction force is the sum of a motor force that initiates contraction, a viscous force that describes the actin–myosin filament sliding, and an elastic force that accounts for the deformation of the cross-bridges. The mechanical response of the collagen fibers is governed by the fiber recruitment process: collagen fibers support load and behave as a linear elastic material only after becoming taut. The proposed structural constitutive model is tested with published active and passive, uniaxial and biaxial experimental data on pig arteries.     

A structural multi-mechanism constitutive equation for cerebral arterial tissue

A structural multi-mechanism constitutive equation is developed to describe the nonlinear, anisotropic, inelastic mechanical behavior of cerebral arterial tissue. Elastin and collagen fibers are treated as separate components (mechanisms) of the artery. Elastin is responsible for load bearing at low strain levels while the collagen mechanism is recruited for load bearing at higher strain levels. This work builds on an earlier model in which both the elastin and collagen mechanisms are represented by isotropic response functions [Wulandana, R., Robertson, A.M., 2005. An inelastic multi-mechanism constitutive equation for cerebral arterial tissue. Biomech. Model. Mechan. 4 (4), 235–248]. Here, the anisotropic material response of the wall is introduced through the collagen mechanism which is composed of helically distributed families of fibers. The orientation of these families is described using either a finite number of fiber orientations or a fiber distribution function. The fiber orientation or dispersion function can be prescribed directly from arterial histology data, or, taking a phenomenological approach, based on data fitting from bi-axial measurements. The activation of the collagen mechanism is specified using a new fiber strain based activation criterion. The multi-mechanism constitutive equation is applied to the simple case of cylindrical inflation and material constants are determined based on available inelastic experimental data for cerebral arteries. While the proposed model captures all features of this inelastic data, there is a pressing need for further experiments to refine the model.     

Elastic–plastic behavior of non-woven fibrous mats

Electrospinning is a novel method for creating non-woven polymer mats that have high surface area and high porosity. These attributes make them ideal candidates for multifunctional composites. Understanding the mechanical properties as a function of fiber properties and mat microstructure can aid in designing these composites. Further, a constitutive model which captures the membrane stress–strain behavior as a function of fiber properties and the geometry of the fibrous network would be a powerful design tool. Here, mats electrospun from amorphous polyamide are used as a model system. The elastic–plastic behavior of single fibers are obtained in tensile tests. Uniaxial monotonic and cyclic tensile tests are conducted on non-woven mats. The mat exhibits elastic–plastic stress–strain behavior. The transverse strain behavior provides important complementary data, showing a negligible initial Poisson's ratio followed by a transverse:axial strain ratio greater than ?1:1 after an axial strain of 0.02. A triangulated framework has been developed to emulate the fibrous network structure of the mat. The micromechanically based model incorporates the elastic–plastic behavior of single fibers into a macroscopic membrane model of the mat. This representative volume element based model is shown to capture the uniaxial elastic–plastic response of the mat under monotonic and cyclic loading. The initial modulus and yield stress of the mat are governed by the fiber properties, the network geometry, and the network density. The transverse strain behavior is linked to discrete deformation mechanisms of the fibrous mat structure including fiber alignment, fiber bending, and network consolidation. The model is further validated in comparison to experiments under different constrained axial loading conditions and found to capture the constraint effect on stiffness, yield, post-yield hardening, and post-yield transverse strain behavior. Due to the direct connection between microstructure and macroscopic behavior, this model should be extendable to other electrospun systems and other two dimensional random fibrous networks.     

A constitutive model for fibrous tissues considering collagen fiber crimp

A micromechanically based constitutive model for fibrous tissues is presented. The model considers the randomly crimped morphology of individual collagen fibers, a morphology typically seen in photomicrographs of tissue samples. It describes the relationship between the fiber endpoints and its arc-length in terms of a measurable quantity, which can be estimated from image data. The collective mechanical behavior of collagen fibers is presented in terms of an explicit expression for the strain-energy function, where a fiber-specific random variable is approximated by a Beta distribution. The model-related stress and elasticity tensors are provided. Two representative numerical examples are analyzed with the aim of demonstrating the peculiar mechanism of the constitutive model and quantifying the effect of parameter changes on the mechanical response. In particular, a fibrous tissue, assumed to be (nearly) incompressible, is subject to a uniaxial extension along the fiber direction, and, separately, to pure shear. It is shown that the fiber crimp model can reproduce several of the expected characteristics of fibrous tissues.     

On the influence of material non-linearities in geometric modeling of kink band instabilities in unidirectional fiber composites

The axial compressive failure of aligned fiber composites triggered by kink band instabilities is the topic of investigation herein. Particular emphasis is put on the accurate prediction of the post-failure regime, where fiber composites are known to exhibit substantial post-failure strength. In this regard, a previous analytical model, based on geometric relationships and energy principles, is enhanced by consistently taking into account material non-linearities. Therefore, a non-linear constitutive law is introduced herein based on a newly developed exponential formulation. This non-linear constitutive law is subsequently implemented into the stress–strain response in interlaminar shearing as well as the compression response. The model enhancements are validated against published experimental data yielding excellent comparisons. Furthermore, the relevance of modeling non-linear material behavior in interlaminar dilation and bending is assessed using a bilinear constitutive law. However, implementing non-linear material behavior does not yield any improvements and can therefore be neglected.     

A mixed Von Mises distribution for modeling soft biological tissues with two distributed fiber properties

Constitutive modeling of biological tissues plays an important role in the understanding of tissue behavior and the development of synthetic materials for medical and bio-inspired applications. A structural continuum model that incorporates principal structural features of the tissue can potentially provide the link between microstructure and the macroscopic mechanical response of biological tissues. For most soft biological tissues, including arterial walls and skin tissue, the main load-carrying constituent is presumed to be the distributed collagen fibers embedded in a base matrix. It is believed that the organization of the collagen fibers gives rise to the anisotropy of the material. In this paper, a semi-structural constitutive model is proposed to account for planar fiber distributions with more than one distributed planar fiber property. Motivated by histology information of the wing membrane of the bat, a statistical treatment is formulated in this paper to capture the overall effect of the distribution of fiber cross-sectional area and the distribution of the number of fibers. This formulation is suitable for general cases when more than one fiber property varies spatially. Furthermore, this model is a two-dimensional specialization within the framework of a three-dimensional theory, which is different the formulation based on a fundamentally two-dimensional theory.     

Homogenization of long fiber reinforced composites including fiber bending effects

This paper presents a homogenization method, which accounts for intrinsic size effects related to the fiber diameter in long fiber reinforced composite materials with two independent constitutive models for the matrix and fiber materials. A new choice of internal kinematic variables allows to maintain the kinematics of the two material phases independent from the assumed constitutive models, so that stress–deformation relationships, can be expressed in the framework of hyper-elasticity and hyper-elastoplasticity for the fiber and the matrix materials respectively. The bending stiffness of the reinforcing fibers is captured by higher order strain terms, resulting in an accurate representation of the micro-mechanical behavior of the composite. Numerical examples show that the accuracy of the proposed model is very close to a non-homogenized finite-element model with an explicit discretization of the matrix and the fibers.     

Continuous versus discrete (invariant) representations of fibrous structure for modeling non-linear anisotropic soft tissue behavior

The non-linear anisotropic mechanical response of soft tissue is largely dependent on the structure of the underlying collagen network. Collagen structure has been successfully quantified for various tissue types in terms of a locally defined fiber orientation distribution function. The continuous distribution function derived from structural data can be directly incorporated into an integral representation of the strain energy function for modeling tissue behavior. Alternatively, non-integral (often invariant-based) strain energy functions have been developed in which the collagen network structure is approximated using a discrete set of fiber classes. The advantage of such an approach is increased computational efficiency since the values of the strain energy and its derivatives (e.g. stress) can be evaluated without numerical integration. However, because of the structural simplifications such models are presumably unable to predict mechanical data as accurately as the models which incorporate a continuous orientation distribution function. In this work the ability of discrete versus continuous fiber models to capture the non-linear anisotropic response of soft tissue is critically analyzed. Both unimodal and bimodal fiber distributions are considered. A general formulation has been developed in terms of an arbitrary fiber strain energy function, such that the analysis can be performed for any suitable fiber material model. For tissue structures in which a discrete representation is suitable, techniques are presented for establishing the range of loading conditions in which model accuracy is not significantly compromised, thus justifying the use of an invariant-based modeling approach.     

Effects of matrix inelasticity on the overall hysteretic behavior of TiNi-SMA fiber composites

The effects of the inelastic deformation of the matrix on the overall hysteretic behavior of a unidirectional titanium–nickel shape-memory alloy (TiNi-SMA) fiber composite and on the local pseudoelastic response of the embedded SMA fibers are studied under the isothermal loading and unloading condition. The multiaxial phase transformation of the SMA fibers is predicted using the phenomenological constitutive equations which can describe the two-step deformation due to the rhombohedral and martensitic transformations, and the inelastic behavior of the matrix material using the standard nonlinear viscoplastic model. The average behavior of the SMA composite is evaluated with the micromechanical method of cells. It is observed that the inelastic deformation of the matrix due to prior tension results in a compressive stress in the matrix after unloading of the SMA composite and this residual stress impedes the complete recovery of the pseudoelastic strain of the SMA fibers. This explains that a closed hysteresis behavior of the SMA composite is no longer observed in contrast with the case that an elastic behavior of matrix is assumed. The predicted local stress–strain behavior indicates that the cyclic response of matrix is crucial to the design of the hysteretic performance of the SMA composite under the repeated loading conditions.     

Visco-electro-elastic models of fiber-distributed active tissues

We present a constitutive model for stochastically distributed fiber reinforced visco-active tissues, where the behavior of the reinforcement depends on the relative orientation of the electric field. Following our previous works, for the passive behaviors we adopt a second order approximation of the strain energy density associated to the parameters of the fiber distribution. Consistently, we also assume that the active behavior accounts for the stochastic distribution of the fibers. The ensuing mechanical quantities result to be dependent on two average structure tensors. We introduce an extended Helmholtz free energy density characterized by the inclusion of a directional active potential, dependent on a stochastic anisotropic permittivity tensor. The permittivity tensor is expanded in Taylor series up to the second order, allowing to obtain an approximated active potential with the same structure of the passive Helmholtz free energy density. In particular, the explicit expression of active stress and stiffness are dependent on the two average structure tensors that characterize the passive response. Anisotropy follows from the fiber distribution and inherits its stochastic nature through statistics parameters. The active fiber distributed model is extended here to viscous materials by including the contribution of a dual dissipation potential in the variational formulation of the constitutive updates. Additionally, we present a computational example of application of the electro-viscous-mechanical material model by simulating peristaltic contractions on a portion of human intestine.     

Micromechanics-based constitutive modeling for unidirectional laminated composites

A micromechanics-based constitutive model is developed to predict the effective mechanical behavior of unidirectional laminated composites. A newly developed Eshelby’s tensor for an infinite circular cylindrical inclusion [Cheng, Z.Q., Batra, R.C., 1999. Exact Eshelby tensor for a dynamic circular cylindrical inclusion. J. Appl. Mech. 66, 563–565] is adopted to model the unidirectional fibers and is incorporated into the micromechanical framework. The progressive loss of strength resulting from the partial fiber debonding and the nucleation of microcracks is incorporated into the constitutive model. To validate the proposed model, the predicted effective stiffness of transversely isotropic composites under far field loading conditions is compared with analytical solutions. The constitutive model incorporating the damage models is then implemented into a finite element code to numerically characterize the elastic behavior of laminated composites. Finally, the present predictions on the stress–strain behavior of laminated composite plate containing an open hole is compared with experimental data to verify the predictive capability of the model.     

Failure criterion of collagen fiber: Viscoelastic behavior simulated by using load control data

A nonlinear Zener model is developed to model the viscoelastic behavior of collagen fibers, a building block of the biological soft tissues in the skeletal system. The effects of the strain rate dependency, the loading history, rest, and recovery on the stress-strain relationship of collagen fibers were investigated using the Zener model. The following loading conditions were simulated: (1) the stress relaxation after cyclic loading, (2) the constant strain rate loading before and after cyclic loading (stabilization) and post recovery, and (3) the constant strain rate loading over a wide range of loading rates. In addition, we explored the critical values of stress and strain using different failure criteria at different strain rates. Four major findings were derived from these simulations. First of all, the stress relaxation is larger with a smaller number cycles of preloading. Second, the strain rate sensitivity diminishes after the stabilization and recovery from resting. Third, the stress-strain curve is dependent on the strain rate except for extreme loading conditions (very fast or slow rates of loading). Finally, the strain energy density (SED) criteria may be a more practical failure criterion than the ultimate stress or strain criterion for collagen fiber. These results provide the basis for interpretation of the viscoelastic and failure behaviors of complex structures such as spinal functional units with more economical CPU than full finite element modeling of the whole structure would have required.     

A continuum mechanics framework and a constitutive model for remodelling of collagen gels and collagenous tissues

Collagen is a very important protein of the human body and is responsible for the structural stability of many body components. Furthermore, collagen fibre networks are able to grow and remodel themselves, which enables them to adjust to varying physiological conditions. This remodelling is accomplished by fibre-producing cells, such as fibroblasts. The ability to adjust to new physiological conditions is very important, for example in wound healing. In the present paper, a theoretical framework for modelling collagenous tissues and collagen gels is proposed. Continuum mechanics is employed to describe the kinematics of the collagen, and affine deformations of fibres are assumed. Biological soft tissues can be approximated as being hyperelastic, and the constitutive model for the collagen fabric is therefore formulated in terms of a strain energy function. This strain energy function includes a density function that describes the distribution of the collagen fibre orientation. The density function evolves according to an evolution law, where fibres tend to reorient towards the direction of maximum Cauchy stress. The remodelling of the collagen network is also assumed to include a pre-stretching of collagen fibres, accomplished by fibroblasts. The theoretical framework is applied to experiments performed on collagen gels, where gels were exposed to remodelling under both biaxial and uniaxial constraints. The proposed model was able to predict both the resulting collagen distribution and the resulting stress-strain relationships obtained for the remodelled collagen gels. The influence of the most important model parameters is demonstrated, and it appears that there is a fairly unique set of model parameters that gives an optimal fit to the experimental data.     

Anisotropic plasticity model coupled with Lode angle dependent strain-induced transformation kinetics law

A phenomenological macroscopic plasticity model is developed for steels that exhibit strain-induced austenite-to-martensite transformation. The model makes use of a stress-state dependent transformation kinetics law that accounts for both the effects of the stress triaxiality and the Lode angle on the rate of transformation. The macroscopic strain hardening is due to nonlinear kinematic hardening as well as isotropic hardening. The latter contribution is assumed to depend on the dislocation density as well as the current martensite volume fraction. The constitutive equations are embedded in the framework of finite strain isothermal rate-independent anisotropic plasticity. Experimental data for an anisotropic austenitic stainless steel 301LN is presented for uniaxial tension, uniaxial compression, transverse plane strain tension and pure shear. The model parameters are identified using a combined analytical–numerical approach. Numerical simulations are performed of all calibration experiments and excellent agreement is observed. Moreover, we make use of experimental data from ten combined tension and shear experiments to validate the proposed constitutive model. In addition, punch and notched tension tests are performed to evaluate the model performance in structural applications with heterogeneous stress and strain fields.     

A stochastic-structurally based three dimensional finite-strain damage model for fibrous soft tissue

A fully three dimensional finite-strain damage model for fibrous soft tissue is developed. The model assumes uncoupled contributions for the matrix and collagen fibers, and uncoupled bulk and deviatoric response over any range of deformations. A simple isotropic damage mechanism within the framework of continuum damage mechanics has been used to describe the softening behavior under deformation for the matrix. On the other hand, statistical aspects related to the length distribution of the reinforcing fibers lead to a damage model for the reinforcing material. As a result, a general theoretical framework for constitutive modeling of biological soft tissue with continuum damage is obtained. A theoretical example consisting of a biaxial test of a soft tissue reinforced with two families of collagen fibers has been considered to demonstrate the capabilities of the proposed model and to study the sensitivity to changes in the statistical parameters associated with the reinforcing material. Also, a preliminary numerical example is included to demonstrate the model on a inhomogeneous boundary value problem. Results show that the model is able to capture the typical stress-strain behavior observed in fibrous soft tissue and seems to confirm the soundness of the proposed formulation.     

Micro–macro analysis of shape memory alloy composites

The aim of the paper is to develop a micro–macro approach for the analysis of the mechanical behavior of composites obtained embedding long fibers of Shape Memory Alloys (SMA) into an elastic matrix. In order to determine the overall constitutive response of the SMA composites, two homogenization techniques are proposed: one is based on the self-consistent method while the other on the analysis of a periodic composite. The overall response of the SMA composites is strongly influenced by the pseudo-elastic and shape memory effects occurring in the SMA material. In particular, it is assumed that the phase transformations in the SMA are governed by the wire temperature and by the average stress tensor acting in the fiber. A possible prestrain of the fibers is taken into account in the model.Numerical applications are developed in order to analyze the thermo-mechanical behavior of the SMA composite. The results obtained by the proposed procedures are compared with the ones determined through a micromechanical analysis of a periodic composite performed using suitable finite elements.Then, in order to study the macromechanical response of structural elements made of SMA composites, a three-dimensional finite element is developed implementing at each Gauss point the overall constitutive laws of the SMA composite obtained by the proposed homogenization procedures. Some numerical applications are developed in order to assess the efficiency of the proposed micro–macro model.     

A computational micromechanics constitutive model for the unloading behavior of paper

In this study, a computational micromechanics material model for the unloading behavior of paper and other nonwoven materials is presented. The asymptotic fiber and bond (AFB) model for paper elastic–plastic behavior [Sinha, S.K., Perkins, R.W., 1995. Micromechanics constitutive model for use in finite element analysis, In: Proceedings of the 1995, Joint ASME Applied Mechanics and Materials Summer Meeting, Los Angeles, CA, USA, Jun 28–30, 1995] has been extended to model the unloading process through a computational algorithm and implemented using the UMAT subroutine in ABAQUS finite element code. For every unloading increment, the material model assumes elastic unloading with a slope equal to the initial elastic modulus. The Jacobian matrix of the constitutive model is updated at every unloading increment by applying the incremental form of AFB model for a planar element with an elastic fiber and bond condition. A uniaxial tensile and a biaxial Mullen burst loading–unloading experiments were carried out for a paperboard sample and simulated using the model. The stress–strain curve and residual strain for the uniaxial loading were in good agreement with experimental results. The finite element model of the burst test with the AFB unloading material model predicted the general shape of the pressure versus deflection curve. However, the model over predicted the residual deflection by more than 50%. The loading portion of the pressure–deflection curve had a significant offset from experimental curves, and the nonlinearity in the unloading curve towards the end was not predicted. The discrepancies with experimental results are attributed to the burst test itself, model parameter estimation inadequacies, boundary conditions used in the FEA, and neglecting time-dependant effects. Nevertheless, the model can be useful in parametric studies relating microstructure to unloading behavior in structural problems.     

Cyclic behavior of unidirectional and cross-ply titanium matrix composites

Relatively simple and efficient micromechanical models are used to obtain the uniaxial response of SCS-6/Timetal 21S with [0]₄ and [0/90]_s laminates when subjected to isothermal and thermomechanical fatigue (TMF) loadings. Features of the modeling that are required to obtain the accurate deformation behavior for this class of materials under these loadings are highlighted. To this end, a comparison is made between the concentric cylinder model and the uniaxial stress model for representing the [0] laminate. The axial stresses from the two models are very similar under mechanical loading. The greatest differences appear under thermal loading alone. The differences on the composite response between a time-independent elastic-plastic and a viscoplastic matrix constitutive model are also examined. The latter is based on the Bodner-Partom unified constitutive model. The [0/90] laminate is treated by adding a parallel element with smeared [90] ply properties to the [0] model and invoking axial strain compatibility as well as stress equilibrium. The proposed constitutive law for the [90] ply includes both matrix viscoplasticity and fiber/matrix separation damage and is based on damage mechanics concepts. The effect of cyclic frequency on TMF behavior is examined. The in-phase TMF life is shown to be very sensitive to frequency due to the relaxation of matrix stress and the attendant increase in fiber stress.     

Influence of random uncertainties of anisotropic fibrous model parameters on arterial pressure estimation

This paper deals with a stochastic approach based on the principle of the maximum entropy to investigate the effect of the parameter random uncertainties on the arterial pressure. Motivated by a hyperelastic, anisotropic, and incompressible constitutive law with fiber families, the uncertain parameters describing the mechanical behavior are considered. Based on the available information, the probability density functions are attributed to every random variable to describe the dispersion of the model parameters. Numerous realizations are carried out, and the corresponding arterial pressure results are compared with the human non-invasive clinical data recorded over a mean cardiac cycle. Furthermore, the Monte Carlo simulations are performed, the convergence of the probabilistic model is proven. The different realizations are useful to define a reliable confidence region, in which the probability to have a realization is equal to 95%. It is shown through the obtained results that the error in the estimation of the arterial pressure can reach 35% when the estimation of the model parameters is subjected to an uncertainty ratio of 5%. Finally, a sensitivity analysis is performed to identify the constitutive law relevant parameters for better understanding and characterization of the arterial wall mechanical behaviors.     

Viscoelasticity of brain corpus callosum in biaxial tension

Computational models of the brain rely on accurate constitutive relationships to model the viscoelastic behavior of brain tissue. Current viscoelastic models have been derived from experiments conducted in a single direction at a time and therefore lack information on the effects of multiaxial loading. It is also unclear if the time-dependent behavior of brain tissue is dependent on either strain magnitude or the direction of loading when subjected to tensile stresses. Therefore, biaxial stress relaxation and cyclic experiments were conducted on corpus callosum tissue isolated from fresh ovine brains. Results demonstrated the relaxation behavior to be independent of strain magnitude, and a quasi-linear viscoelastic (QLV) model was able to accurately fit the experimental data. Also, an isotropic reduced relaxation tensor was sufficient to model the stress-relaxation in both the axonal and transverse directions. The QLV model was fitted to the averaged stress relaxation tests at five strain magnitudes while using the measured strain history from the experiments. The resulting model was able to accurately predict the stresses from cyclic tests at two strain magnitudes. In addition to deriving a constitutive model from the averaged experimental data, each specimen was fitted separately and the resulting distributions of the model parameters were reported and used in a probabilistic analysis to determine the probability distribution of model predictions and the sensitivity of the model to the variance of the parameters. These results can be used to improve the viscoelastic constitutive models used in computational studies of the brain.