On arterial fiber dispersion and auxetic effect

There are two polar contemporary approaches to the constitutive modeling of arterial wall with anisotropy induced by collagen fibers. The first one is based on the angular integration (AI) of the strain energy on a unit sphere for the analytically defined fiber dispersion. The second one is based on the introduction of the generalized structure tensors (GST). AI approach is very involved computationally while GST approach requires somewhat complicated procedure for the exclusion of compressed fibers.We present some middle ground models, which are based on the use of 16 and 8 structure tensors. These models are moderately involved computationally and they allow excluding compressed fibers easily. We use the proposed models to study the role of the fiber dispersion in the constitutive modeling of the arterial wall. Particularly, we study the auxetic effect which can appear in anisotropic materials. The effect means thickening of the tissue in the direction perpendicular to its stretching. Such an effect was not observed in experiments while some simple anisotropic models do predict it. We show that more accurate account of the fiber dispersion suppresses the auxetic effect in a qualitative agreement with experimental observations.     

Anisotropy of fibrous tissues in relation to the distribution of tensed and buckled fibers

Fibrous tissues are characterized by a much higher stiffness in tension than compression. This study uses microstructural modeling to analyze the material symmetry of fibrous tissues undergoing tension and compression, to better understand how material symmetry relates to the distribution of tensed and buckled fibers. The analysis is also used to determine whether the behavior predicted from a microstructural model can be identically described by phenomenological continuum models. The analysis confirms that in the case when all the fibers are in tension in the current configuration, the material symmetry of a fibrous tissue in the corresponding reference configuration is dictated by the symmetry of its fiber angular distribution in that configuration. However, if the strain field exhibits a mix of tensile and compressive principal normal strains, the fibrous tissue is represented by a material body which consists only of those fibers which are in tension; the material symmetry of this body may be deduced from the superposition of the planes of symmetry of the strain and the planes of symmetry of the angular fiber distribution. Thus the material symmetry is dictated by the symmetry of the angular distribution of only those fibers which are in tension. Examples are provided for various fiber angular distribution symmetries. In particular, it is found that a fibrous tissue with isotropic fiber angular distribution exhibits orthotropic symmetry when subjected to a mix of tensile and compressive principal normal strains, with the planes of symmetry normal to the principal directions of the strain. This anisotropy occurs even under infinitesimal strains and is distinct from the anisotropy induced from the finite rotation of fibers. It is also noted that fibrous materials are not stable under all strain states due to the inability of fibers to sustain compression along their axis; this instability can be overcome by the incorporation of a ground matrix. It is concluded that the material response predicted using a microstructural model of the fibers cannot be described exactly by phenomenological continuum models. These results are also applicable to nonbiological fiber-composite materials.     

Role of elastin anisotropy in structural strain energy functions of arterial tissue

The vascular wall exhibits nonlinear anisotropic mechanical properties. The identification of a strain energy function (SEF) is the preferred method to describe its complex nonlinear elastic properties. Earlier constituent-based SEF models, where elastin is modeled as an isotropic material, failed in describing accurately the tissue response to inflation–extension loading. We hypothesized that these shortcomings are partly due to unaccounted anisotropic properties of elastin. We performed inflation–extension tests on common carotid of rabbits before and after enzymatic degradation of elastin and applied constituent-based SEFs, with both an isotropic and an anisotropic elastin part, on the experimental data. We used transmission electron microscopy (TEM) and serial block-face scanning electron microscopy (SBFSEM) to provide direct structural evidence of the assumed anisotropy. In intact arteries, the SEF including anisotropic elastin with one family of fibers in the circumferential direction fitted better the inflation–extension data than the isotropic SEF. This was supported by TEM and SBFSEM imaging, which showed interlamellar elastin fibers in the circumferential direction. In elastin-degraded arteries, both SEFs succeeded equally well in predicting anisotropic wall behavior. In elastase-treated arteries fitted with the anisotropic SEF for elastin, collagen engaged later than in intact arteries. We conclude that constituent-based models with an anisotropic elastin part characterize more accurately the mechanical properties of the arterial wall when compared to models with simply an isotropic elastin. Microstructural imaging based on electron microscopy techniques provided evidence for elastin anisotropy. Finally, the model suggests a later and less abrupt collagen engagement after elastase treatment.     

Collagen fiber alignment does not explain mechanical anisotropy in fibroblast populated collagen gels

Many load-bearing soft tissues exhibit mechanical anisotropy. In order to understand the behavior of natural tissues and to create tissue engineered replacements, quantitative relationships must be developed between the tissue structures and their mechanical behavior. We used a novel collagen gel system to test the hypothesis that collagen fiber alignment is the primary mechanism for the mechanical anisotropy we have reported in structurally anisotropic gels. Loading constraints applied during culture were used to control the structural organization of the collagen fibers of fibroblast populated collagen gels. Gels constrained uniaxially during culture developed fiber alignment and a high degree of mechanical anisotropy, while gels constrained biaxially remained isotropic with randomly distributed collagen fibers. We hypothesized that the mechanical anisotropy that developed in these gels was due primarily to collagen fiber orientation. We tested this hypothesis using two mathematical models that incorporated measured collagen fiber orientations: a structural continuum model that assumes affine fiber kinematics and a network model that allows for nonaffine fiber kinematics. Collagen fiber mechanical properties were determined by fitting biaxial mechanical test data from isotropic collagen gels. The fiber properties of each isotropic gel were then used to predict the biaxial mechanical behavior of paired anisotropic gels. Both models accurately described the isotropic collagen gel behavior. However, the structural continuum model dramatically underestimated the level of mechanical anisotropy in aligned collagen gels despite incorporation of measured fiber orientations; when estimated remodeling-induced changes in collagen fiber length were included, the continuum model slightly overestimated mechanical anisotropy. The network model provided the closest match to experimental data from aligned collagen gels, but still did not fully explain the observed mechanics. Two different modeling approaches showed that the level of collagen fiber alignment in our uniaxially constrained gels cannot explain the high degree of mechanical anisotropy observed in these gels. Our modeling results suggest that remodeling-induced redistribution of collagen fiber lengths, nonaffine fiber kinematics, or some combination of these effects must also be considered in order to explain the dramatic mechanical anisotropy observed in this collagen gel model system.     

Microplane constitutive model and computational framework for blood vessel tissue

This paper presents a nonlinearly elastic anisotropic microplane formulation in 3D for computational constitutive modeling of arterial soft tissue in the passive regime. The constitutive modeling of arterial (and other biological) soft tissue is crucial for accurate finite element calculations, which in turn are essential for design of implants, surgical procedures, bioartificial tissue, as well as determination of effect of progressive diseases on tissues and implants. The model presented is defined at a lower scale (mesoscale) than the conventional macroscale and it incorporates the effect of all the (collagen) fibers which are anisotropic structural components distributed in all directions within the tissue material in addition to that of isotropic bulk tissue. It is shown that the proposed model not only reproduces Holzapfel's recent model but also improves on it by accounting for the actual three-dimensional distribution of fiber orientation in the arterial wall, which endows the model with advanced capabilities in simulation of remodeling of soft tissue. The formulation is flexible so that its parameters could be adjusted to represent the arterial wall either as a single material or a material composed of several layers in finite element analyses of arteries. Explicit algorithms for both the material subroutine and the explicit integration with dynamic relaxation of equations of motion using finite element method are given. To circumvent the slow convergence of the standard dynamic relaxation and small time steps dictated by the stability of the explicit integrator, an adaptive dynamic relaxation technique that ensures stability and fastest possible convergence rates is developed. Incompressibility is enforced using penalty method with an updated penalty parameter. The model is used to simulate experimental data from the literature demonstrating that the model response is in excellent agreement with the data. An experimental procedure to determine the distribution of fiber directions in 3D for biological soft tissue is suggested in accordance with the microplane concept. It is also argued that this microplane formulation could be modified or extended to model many other phenomena of interest in biomechanics.     

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.     

Characterizing the mechanical contribution of fiber angular distribution in connective tissue: comparison of two modeling approaches

Modeling of connective tissues often includes collagen fibers explicitly as one of the components. These fibers can be oriented in many directions; therefore, several studies have considered statistical distributions to describe the fiber arrangement. One approach to formulate a constitutive framework for distributed fibers is to express the mechanical parameters, such as strain energy and stresses, in terms of angular integrals. These integrals represent the addition of the contribution of infinitesimal fractions of fibers oriented in a given direction. This approach leads to accurate results; however, it requires lengthy calculations. Recently, the use of generalized structure tensors has been proposed to represent the angular distribution in the constitutive equations of the fibers. Although this formulation is much simpler and fewer calculations are required, such structure tensors can only be used when all the fibers are in tension and the angular distribution is small. However, the amount of error introduced in these cases of non-tensile fiber loading and large angular distributions have not been quantified. Therefore, the objective of this study is to determine the range of values of angular distribution for which acceptable differences (less than 10%) between these two formulations are obtained. It was found, analytically and numerically, that both formulations are equivalent for planar distributions under equal-biaxial stretch. The comparison also showed, for other loading conditions, that the differences decrease when the fiber distribution is very small. Differences of less than 10% were usually obtained when the fiber distribution was very low (κ ≈ 0.03; κ ranges between 0 and 1/3, for aligned and isotropic distributed fibers, respectively). This range of angular distribution greatly limits the types of tissue that can be accurately analyzed using generalized structure tensors. It is expected that the results from this study guide the selection of a proper approach to analyze a particular tissue under a particular loading condition.     

New interpretation of arterial stiffening due to cigarette smoking using a structurally motivated constitutive model

Cigarette smoking is the leading self-inflicted risk factor for cardiovascular diseases; it causes arterial stiffening with serious sequelea including atherosclerosis and abdominal aortic aneurysms. This work presents a new interpretation of arterial stiffening caused by smoking based on data published for rat pulmonary arteries. A structurally motivated "four fiber family" constitutive relation was used to fit the available biaxial data and associated best-fit values of material parameters were estimated using multivariate nonlinear regression. Results suggested that arterial stiffening caused by smoking was reflected by consistent increase in an elastin-associated parameter and moreover by marked increase in the collagen-associated parameters. That is, we suggest that arterial stiffening due to cigarette smoking appears to be isotropic, which may allow simpler phenomenological models to capture these effects using a single stiffening parameter similar to the approach in isotropic continuum damage mechanics. There is a pressing need, however, for more detailed histological information coupled with more complete biaxial mechanical data for a broader range of systemic arteries.     

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.     

Passive material properties of intact ventricular myocardium determined from a cylindrical model

The equatorial region of the canine left ventricle was modeled as a thick-walled cylinder consisting of an incompressible hyperelastic material with homogeneous exponential properties. The anisotropic properties of the passive myocardium were assumed to be locally transversely isotropic with respect to a fiber axis whose orientation varied linearly across the wall. Simultaneous inflation, extension, and torsion were applied to the cylinder to produce epicardial strains that were measured previously in the potassium-arrested dog heart. Residual stress in the unloaded state was included by considering the stress-free configuration to be a warped cylindrical arc. In the special case of isotropic material properties, torsion and residual stress both significantly reduced the high circumferential stress peaks predicted at the endocardium by previous models. However, a resultant axial force and moment were necessary to cause the observed epicardial deformations. Therefore, the anisotropic material parameters were found that minimized these resultants and allowed the prescribed displacements to occur subject to the known ventricular pressure loads. The global minimum solution of this parameter optimization problem indicated that the stiffness of passive myocardium (defined for a 20 percent equibiaxial extension) would be 2.4 to 6.6 times greater in the fiber direction than in the transverse plane for a broad range of assumed fiber angle distributions and residual stresses. This agrees with the results of biaxial tissue testing. The predicted transmural distributions of fiber stress were relatively flat with slight peaks in the subepicardium, and the fiber strain profiles agreed closely with experimentally observed sarcomere length distributions. The results indicate that torsion, residual stress and material anisotropy associated with the fiber architecture all can act to reduce endocardial stress gradients in the passive left ventricle.     

Mechanical stresses in abdominal aortic aneurysms: influence of diameter, asymmetry, and material anisotropy

Biomechanical studies suggest that one determinant of abdominal aortic aneurysm (AAA) rupture is related to the stress in the wall. In this regard, a reliable and accurate stress analysis of an in vivo AAA requires a suitable 3D constitutive model. To date, stress analysis conducted on AAA is mainly driven by isotropic tissue models. However, recent biaxial tensile tests performed on AAA tissue samples demonstrate the anisotropic nature of this tissue. The purpose of this work is to study the influence of geometry and material anisotropy on the magnitude and distribution of the peak wall stress in AAAs. Three-dimensional computer models of symmetric and asymmetric AAAs were generated in which the maximum diameter and length of the aneurysm were individually controlled. A five parameter exponential type structural strain-energy function was used to model the anisotropic behavior of the AAA tissue. The anisotropy is determined by the orientation of the collagen fibers (one parameter of the model). The results suggest that shorter aneurysms are more critical when asymmetries are present. They show a strong influence of the material anisotropy on the magnitude and distribution of the peak stress. Results confirm that the relative aneurysm length and the degree of aneurysmal asymmetry should be considered in a rupture risk decision criterion for AAAs.     

How to asses, visualize and compare the anisotropy of linear structures reconstructed from optical sections--a study based on histopathological quantification of human brain microvessels

Three-dimensional analyses of the spatial arrangement, spatial orientation and preferential directions of systems of fibers are frequent tasks in many scientific fields, including the textile industry, plant biology and tissue modeling. In biology, systems of oriented and branching lines are often used to represent the three-dimensional directionality and topology of microscopic blood vessels supplying various organs. In our study, we present a novel p(χ2) (chi-square) method for evaluating the anisotropy of line systems that involves comparing the observed length densities of lines with the discrete uniform distribution of an isotropic line system with the χ2-test. Using this method in our open source software, we determined the rose of directions, preferential directions and level of anisotropy of linear systems representing the microscopic blood vessels in samples of various regions from human brains (cortex, subcortical gray matter and white matter). The novel method was compared with two other methods used for anisotropy quantification (ellipsoidal and fractional anisotropy). All three methods detected different levels of anisotropy of blood microvessels in human brain. The microvascular bed in the cortex was closer to an isotropic network, while the microvessels supplying the white matter appeared to be an anisotropic and direction-sensitive system. All three methods were able to determine the differences between various brain regions. The advantage of our p(χ2) method is its high correlation with the number of preferential directions of the line system. However, the software, named esofspy, is able to calculate all three of the measures of anisotropy compared and documented in this paper, thus making the methods freely available to the scientific community.     

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.     

Biaxial mechanical properties of the native and glutaraldehyde-treated aortic valve cusp: Part II--A structural constitutive model

We have formulated the first constitutive model to describe the complete measured planar biaxial stress-strain relationship of the native and glutaraldehyde-treated aortic valve cusp using a structurally guided approach. When applied to native, zero-pressure fixed, and low-pressure fixed cusps, only three parameters were needed to simulate fully the highly anisotropic, and nonlinear in-plane biaxial mechanical behavior. Differences in the behavior of the native and zero- and low-pressure fixed cusps were found to be primarily due to changes in the effective fiber stress-strain behavior. Further, the model was able to account for the effects of small (< 10 deg) misalignments in the cuspal specimens with respect to the biaxial test axes that increased the accuracy of the model material parameters. Although based upon a simplified cuspal structure, the model underscored the role of the angular orientation of the fibers that completely accounted for extreme mechanical anisotropy and pronounced axial coupling. Knowledge of the mechanics of the aortic cusp derived from this model may aid in the understanding of fatigue damage in bioprosthetic heart valves and, potentially, lay the groundwork for the design of tissue-engineered scaffolds for replacement heart valves.     

Reciprocal influence of ethylene and gibberellins on response-gene expression in <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis>

A cortical band of fiber cells originate de novo in tendrils of redvine [Brunnichia ovata (Walt.) Shiners] when these convert from straight, supple young filaments to stiffened coiled structures in response to touch stimulation. We have analyzed the cell walls of these fibers by in situ localization techniques to determine their composition and possible role(s) in the coiling process. The fiber cell wall consists of a primary cell wall and two lignified secondary wall layers (S₁ and S₂) and a less lignified gelatinous (G) layer proximal to the plasmalemma. Compositionally, the fibers are sharply distinct from surrounding parenchyma as determined by antibody and affinity probes. The fiber cell walls are highly enriched in cellulose, callose and xylan but contain no homogalacturonan, either esterified or de-esterified. Rhamnogalacturonan-I (RG-I) epitopes are not detected in the S layers, although they are in both the gelatinous layer and primary wall, indicating a further restriction of RG-I in the fiber cells. Lignin is concentrated in the secondary wall layers of the fiber and the compound middle lamellae/primary cell wall but is absent from the gelatinous layer. Our observations indicate that these fibers play a central role in tendril function, not only in stabilizing its final shape after coiling but also generating the tensile strength responsible for the coiling. This theory is further substantiated by the absence of gelatinous layers in the fibers of the rare tendrils that fail to coil. These data indicate that gelatinous-type fibers are responsible for the coiling of redvine tendrils and a number of other tendrils and vines.     

A Gauss-Kronrod-Trapezoidal integration scheme for modeling biological tissues with continuous fiber distributions

Fibrous biological tissues may be modeled using a continuous fiber distribution (CFD) to capture tension–compression nonlinearity, anisotropic fiber distributions, and load-induced anisotropy. The CFD framework requires spherical integration of weighted individual fiber responses, with fibers contributing to the stress response only when they are in tension. The common method for performing this integration employs the discretization of the unit sphere into a polyhedron with nearly uniform triangular faces (finite element integration or FEI scheme). Although FEI has proven to be more accurate and efficient than integration using spherical coordinates, it presents three major drawbacks: First, the number of elements on the unit sphere needed to achieve satisfactory accuracy becomes a significant computational cost in a finite element (FE) analysis. Second, fibers may not be in tension in some regions on the unit sphere, where the integration becomes a waste. Third, if tensed fiber bundles span a small region compared to the area of the elements on the sphere, a significant discretization error arises. This study presents an integration scheme specialized to the CFD framework, which significantly mitigates the first drawback of the FEI scheme, while eliminating the second and third completely. Here, integration is performed only over the regions of the unit sphere where fibers are in tension. Gauss–Kronrod quadrature is used across latitudes and the trapezoidal scheme across longitudes. Over a wide range of strain states, fiber material properties, and fiber angular distributions, results demonstrate that this new scheme always outperforms FEI, sometimes by orders of magnitude in the number of computational steps and relative accuracy of the stress calculation.     

A computational study of invariant I5 in a nearly incompressible transversely isotropic model for white matter

The aligned axonal fiber bundles in white matter make it suitable to be modeled as a transversely isotropic material. Recent experimental studies have shown that a minimal form, nearly incompressible transversely isotropic (MITI) material model, is capable of describing mechanical anisotropy of white matter. Here, we used a finite element (FE) computational approach to demonstrate the significance of the fifth invariant (I₅) when modeling the anisotropic behavior of white matter in the large-strain regime. We first implemented and validated the MITI model in an FE simulation framework for large deformations. Next, we applied the model to a plate-hole structural problem to highlight the significance of the invariant I₅ by comparing with the standard fiber reinforcement (SFR) model. We also compared the two models by fitting the experiment data of asymmetric indentation, shear test, and uniaxial stretch of white matter. Our results demonstrated the significance of I₅ in describing shear deformation/anisotropy, and illustrated the potential of the MITI model to characterize transversely isotropic white matter tissues in the large-strain regime.     

Cellulose Content in Fibers of Cottons which Differ in their Lint Lengths and Extent of Fuzz

Three cottons differing in their extent of fuzz fibers (linters) and final length of lint fibers were analyzed for amount of fiber cell walls and fiber cellulose at various times postanthesis. Cellulose determinations were performed directly on whole fibers and on fiber cell wall preparations. The data suggest that the presence of fuzz fibers does not account for a rise, followed by a plateau, followed by a rise, in cellulose content expressed as a percentage of cell wall material. It is concluded that: (1) under our greenhouse conditions, all fuzz fibers are initiated by day eight after anthesis; (2) weight per mm length of all fibers increases up to the point of secondary wall deposition and increases even more rapidly after that; (3) deposition of secondary wall cellulose in fuzz fibers probably does not begin until after similar deposition begins in lint fibers; (4) the actual amount of cellulose in primary walls of all elongating fibers (fuzz and lint) is a constant value, about 1 × 10^?16 mg/mm; and (5) secondary wall cellulose deposition in lint fibers begins very sharply, in advance of cessation of elongation, at a time closely related to final lint fiber length. It is speculated that: (1) cell wall preparation procedures may remove significant amounts of noncellulosic wall material, thus making it difficult to define all functional constituents on the basis of what is left in a cell wall residue; and/or (2) primary walls may lose to the cytoplasm some of their constituents in advance of secondary wall deposition, the extent of loss varying due to developmental age of the elongating fibers.     

A new constitutive formulation for characterizing the mechanical behavior of soft tissues.

We present a new constitutive formulation that combines certain desirable features of two previously used approaches (phenomenological and microstructural). Specifically, we assume that certain soft tissues can be idealized as composed of various families of noninteracting fibers and a homogeneous matrix. Both the fibers and the matrix are assumed to follow the gross deformation. Within the usual framework of pseudoelasticity, incompressibility, homogeneity, and the continuum hypothesis, a pseudostrain-energy function (W) is proposed wherein W is expressed in terms of matrix and fibrous contributions. Unlike phenomenological approaches where a W is usually chosen in an ad hoc manner, the present approach can be used to postulate reasonable forms of W based on limited structural information and multiaxial stress-strain data. Illustrative applications of the theory are discussed for visceral pleura and myocardium. Concise structurally motivated constitutive relations result, wherein load-dependent anisotropy, nonlinear material behavior, finite deformations, and incompressibility are accounted for.