Stress-free Reference Configuration Research Articles

Stimuli-responsive shell theory

Soft matter mechanics generally involve finite deformations and instabilities of structures in response to a wide range of mechanical and non-mechanical stimuli. Modeling plates and shells is generally a challenge due to their geometrically nonlinear response to loads; however, non-mechanical loads further complicate matters as it is often not clear how they modify the shell’s energy functional. In this work, we demonstrate how to form a mechanical interpretation of these non-mechanical stimuli, in which the standard shell strain energy can be augmented with potentials corresponding to how a non-mechanical stimulus acts to change the shell’s area and curvature via the natural stretch and curvature. As a result, the effect of non-mechanical stimuli to deform shells is transformed into effective external loadings, and this framework allows for the application of analytical and computational tools that are standard within the mechanics community. Furthermore, we generalize the effect of mass change during biological growth to account for its effect on the stress constitution. The theory is formulated based on a standard, stress-free reference configuration which is known a priori, meaning it can be physically observed, and only requires the solution of a single-field equation, the standard mechanical momentum or equilibrium equation, despite capturing the effects of non-mechanical stimuli. We validate the performance of this model by several benchmark problems, and finally, we apply it to complex examples, including the snapping of the Venus flytrap, leaf growth, and the buckling of electrically active polymer plates. Overall, we expect that mechanicians and non-mechanicians alike can use the approach presented here to quickly modify existing computational tools with effective external loadings calculated in this novel theory to study how various types of non-mechanical stimuli impact the mechanics and physics of thin shell structures.

NN-EUCLID: Deep-learning hyperelasticity without stress data

We propose a new approach for unsupervised learning of hyperelastic constitutive laws with physics-consistent deep neural networks. In contrast to supervised learning, which assumes the availability of stress–strain pairs, the approach only uses realistically measurable full-field displacement and global reaction force data, thus it lies within the scope of our recent framework for Efficient Unsupervised Constitutive Law Identification and Discovery (EUCLID) and we denote it as NN-EUCLID. The absence of stress labels is compensated for by leveraging a physics-motivated loss function based on the conservation of linear momentum to guide the learning process. The constitutive model is based on input-convex neural networks, which are capable of learning a function that is convex with respect to its inputs. By employing a specially designed neural network architecture, multiple physical and thermodynamic constraints for hyperelastic constitutive laws, such as material frame indifference, material stability, and stress-free reference configuration are automatically satisfied. We demonstrate the ability of the approach to accurately learn several hidden isotropic and anisotropic hyperelastic constitutive laws – including e.g., Mooney–Rivlin, Arruda–Boyce, Ogden, and Holzapfel models – without using stress data. For anisotropic hyperelasticity, the unknown anisotropic fiber directions are automatically discovered jointly with the constitutive model. The neural network-based constitutive models show good generalization capability beyond the strain states observed during training and are readily deployable in a general finite element framework for simulating complex mechanical boundary value problems with good accuracy.

Benchtop characterization of the tricuspid valve leaflet pre-strains

The pre-strains of biological soft tissues are important when relating their in vitro and in vivo mechanical behaviors. In this study, we present the first-of-its-kind experimental characterization of the tricuspid valve leaflet pre-strains. We use 3D photogrammetry and the reproducing kernel method to calculate the pre-strains within the central 10×10 mm region of the tricuspid valve leaflets from n=8 porcine hearts. In agreement with previous pre-strain studies for heart valve leaflets, our results show that all the three tricuspid valve leaflets shrink after being explanted from the ex vivo heart. These calculated strains are leaflet-specific and the septal leaflet experiences the most compressive changes. Furthermore, the strains observed after dissection of the central 10×10 mm region of the leaflet are smaller than when the valve is explanted, suggesting that our computed pre-strains are mainly due to the release of in situ annulus and chordae connections. The leaflets are then mounted on a biaxial testing device and preconditioned using force-controlled equibiaxial loading. We show that the employed preconditioning protocol does not 100% restore the leaflet pre-strains as removed during tissue dissection, and future studies are warranted to explore alternative preconditioning methods. Finally, we compare the calculated biomechanically oriented metrics considering five stress-free reference configurations. Interestingly, the radial tissue stretches and material anisotropies are significantly smaller compared to the post-preconditioning configuration. Extensions of this work can further explore the role of this unique leaflet-specific leaflet pre-strains on in vivo valve behavior via high-fidelity in-silico models. Statement of significanceThis study provides a first of its kind benchtop characterization of tricuspid valve leaflet pre-strains. We used 3D photogrammetry to reconstruct the central region of the tricuspid valve leaflets in three configurations. The associated configurational changes revealed compressive, leaflet-specific strains after dissection of the valve from its in situ environment. Interestingly, we found that biaxial preconditioning did not restore the measured pre-strains of the leaflets. Depending on the selection of the stress-free reference configuration, this led to substantial differences in the leaflet mechanics. Our findings and methodology are crucial when it comes to relating in vitro mechanical behaviors to valve function in vivo. Future studies can integrate our quantified pre-strains into in-silico simulations to get more realistic predictions about the valve function.

A cardiac electromechanical model coupled with a lumped-parameter model for closed-loop blood circulation

We propose a novel mathematical and numerical model for cardiac electromechanics, wherein biophysically detailed core models describe the different physical processes concurring to the cardiac function. The core models, once suitably approximated, are coupled by a computationally efficient strategy. Our model is based on: (1) the combination of implicit-explicit (IMEX) schemes to solve the different core cardiac models, (2) an Artificial Neural Network based model, that surrogates a biophysically detailed but computationally demanding microscale model of active force generation and (3) appropriate partitioned schemes to couple the different models in this multiphysics setting. We employ a flexible and scalable intergrid transfer operator, which allows to interpolate Finite Element functions between nested meshes and, possibly, among arbitrary Finite Element spaces for the different core models. Our core 3D electromechanical model of the left ventricle is coupled with a closed-loop 0D model of the vascular network (and the other cardiac chambers) by an approach that is energy preserving. More precisely, we derive a balance law for the mechanical energy of the whole circulatory network. This provides a quantitative insight into the energy utilization, dissipation and transfer among the different compartments of the cardiovascular network and during different stages of the heartbeat. On this ground, a new tool is proposed to validate some energy indicators adopted in the daily clinical practice. A further contribution of this paper is the proposition of a robust algorithm for the reconstruction of the stress-free reference configuration. This feature is fundamental to correctly initialize our electromechanical simulations. As a matter of fact, the geometry acquired from medical imaging typically refers to a configuration affected by residual internal stresses, whereas the elastodynamics equations that govern the mechanics core model are related to a stress-free configuration. To prove the biophysical accuracy of our computational model, we address different scenarios of clinical interest, namely by varying preload, afterload and contractility.

ANCF Multiplicative-Decomposition Thermoelastic Approach for Arbitrary Geometry

The classical approach for the thermal analysis of solids and fluids uses the strain additive decomposition to account for the effect of the temperature. This strain-based approach does not properly capture the effect of complex stress-free reference-configuration geometry, is applicable only to small deformation problems, and leads to simplified expression for the Green–Lagrange strain tensor. In view of the geometric description of the absolute nodal coordinate formulation (ANCF), a new ANCF gradient-based approach is proposed. This approach employs a multiplicative decomposition of the matrix of position-gradient vectors in the stress-free reference configuration into two position-gradient matrices. One matrix is associated with the reference-configuration geometry before the application of the thermal load, and the other accounts for the volumetric change due to the change in temperature. A numerical study demonstrated the implementation of the proposed gradient-based approach.

A geometrically accurate deformable-body approach for the analysis of robotic and parallel-mechanism systems

Strength and durability investigations of parallel-mechanism (PM) and robotic systems, used as integral part in modern manufacturing, are necessary at the virtual design stage. Nonetheless, because of the limitations of existing finite element (FE) approaches, flexible-link geometries are often simplified or distorted when converting solid models to FE analysis meshes. This article introduces and demonstrates the use of a new unified geometry/analysis approach for the small-deformation analysis of robotic and parallel mechanism (RPM) systems with flexible links. The approach used in this investigation integrates the geometrically accurate absolute nodal coordinates formulation (ANCF) with the computationally efficient floating frame of reference (FFR) formulation that allows for systematic elimination of high-frequency and insignificant deformation modes. The proposed ANCF/FFR approach allows modeling accurately the stress-free reference configuration geometry, captures the dynamic coupling between the rigid-body and elastic displacements, is based on a unified geometry/analysis mesh that eliminates the need for geometry/analysis model conversion or coordinate transformation, and allows for efficient and accurate solution of the nonlinear dynamic equations that govern the RPM motion. The formulation of the spatial ANCF/FFR equations of motion including the elastic (stress) forces using a general continuum-mechanics approach is presented, and use of the strain split method (SSM) as an FE locking-alleviation technique is discussed. A frequency-convergence analysis of flexible links with tapered geometry is performed and the obtained numerical results are compared with solutions obtained using commercial FE software. The application of the new procedure in the nonlinear dynamics simulations of RPM systems is demonstrated using a spatial parallel mechanism that includes flexible links and prismatic, revolute, and spherical joints.

Modeling the non-dissipative response of vulcanized unfilled rubber

Here the possibility of developing a better three dimensional model for rubber when Hencky strain is written as an explicit function of the Cauchy stress is examined. Based on the experimental evidence, vulcanized unfilled rubber is assumed to be isotropic, incompressible, and deforming from a stress-free reference configuration. Vulcanized unfilled rubber undergoes strain-induced crystallization when uniaxial stretch exceeds 2. Therefore, in this study, a non-dissipative model is sought for rubber only till the von Mises stress is below a critical value. Using the biaxial response reported in Kawabata et al. (1981), the material functions are arrived adopting Rivlin and Saunders methodology. The material constants in the determined material functions are found using biaxial tests for four other experimental studies on rubber reported in the literature. The ability of the arrived constitutive model to describe the uniaxial and equal biaxial response of the same rubber is evaluated. The proposed model is benchmarked with established models in the classical framework. Among the models studied the proposed model performed as well as the other models for Kawabata et al. data set and marginally better than the other models for the remaining four experimental data sets. However, while the proposed model performance does not critically depend on the experiments used to deduce the material parameters, for the other models the material parameters need to be found from fixed biaxial experiments for better performance.

On process-step parallel computability and linear superposition of mechanical responses in additive manufacturing process simulation

We study linearity assumptions in the transient macroscale mechanical aspect of additive manufacturing (AM) process simulation. Linearity assumptions are often resorted to in combination with calibrated inelastic deformation components to arrive at computationally tractable yet reasonably accurate AM process models. We point out that linearity assumptions permit the independent computation of the response increment in each step of the AM process, and the total mechanical response is the superposition of all the process-step increments. In effect, process-step increments are computed with respect to the stress-free reference configuration in each step. The implication is that the mechanical response increment in each linearised AM process step may be computed in parallel. Trivial process-step-wise parallel computability breaks down, however, if nonlinearity (i.e. geometric or material) is modelled. In our investigation the influence of geometric nonlinearity on part distortion is small (but this is of course part-geometry specific), and more realistic stresses are obtained by imposing a nonlinear elastoplastic material law after the parallel computation and superposition of the linear AM response increments. It is demonstrated that simulation wall-clock time is reduced by exploiting process-step parallel computability in the linear regime. Moreover, numerical experiments suggest that process-step parallelization scales better (in wall-clock time) than conventional parallelization in the sequential computation of each response increment.

Constrained mixture modeling affects material parameter identification from planar biaxial tests

The constrained mixture theory is an elegant way to incorporate the phenomenon of residual stresses in patient-specific finite element models of arteries. This theory assumes an in vivo reference geometry, obtained from medical imaging, and constituent-specific deposition stretches in the assumed reference state. It allows to model residual stresses and prestretches in arteries without the need for a stress-free reference configuration, most often unknown in patient-specific modeling.A finite element (FE) model requires material parameters, which are classically obtained by fitting the constitutive model to experimental data. The characterization of arterial tissue is often based on planar biaxial test data, to which nonlinear elastic fiber-reinforced material parameters are fitted. However, the introduction of the constrained mixture theory requires an adapted approach to parameter fitting. Therefore, we introduce an iterative fitting method, alternating between nonlinear least squares parameter optimization and an FE prestressing algorithm to obtain the correct constrained mixture material state during the mechanical test.We verify the method based on numerically constructed planar biaxial test data sets, containing ground truth sets of material parameters. The results show that the method converges to the correct parameter sets in just a few iterations. Next, the iterative fitting approach is applied to planar biaxial test data of ovine pulmonary artery tissue. The obtained results demonstrate a convergence towards constrained mixture compatible parameters, which differ significantly from classically obtained parameters.We show that this new modeling approach yields in vivo wall stresses similar to when using classically obtained parameters. However, due to the numerous advantages of constrained mixture modeling, our fitting method is relevant to obtain compatible material parameters, that may not be confused with parameters obtained in a classical way.

A three dimensional constitutive model for plain cement concrete

The validity of some of the assumptions made during systematic development of constitutive relation – homogeneity, isotropy, and non-dissipative response – is examined experimentally by testing plain concrete cylinders in the circumferential displacement controlled uniaxial compression tests. The tested cylinders indicate that the surface strains are not uniform and that the principal direction of the strain varies with the magnitude of the applied load. Hence, these cylinders are either deforming from a stressed state or anisotropic or both. The percentage dissipation computed as the area between the loading and unloading curves normalized using the area under the unloading curve is less than 10 percent up to 90 percent of the peak load. Therefore the mechanical response of concrete can be considered as non-dissipative. Assuming that the magnitude of the residual stresses present in the concrete panels would be small compared to the stresses arising due to the applied load, these residual stresses are ignored. Consequently, using the implicit constitutive theory framework for isotropic and compressible materials undergoing a non-dissipative process from a stress-free reference configuration, a three-dimensional constitutive relation for plain concrete is proposed. The material parameters in the proposed model are estimated from the uniaxial compression test and equal biaxial compression test. These parameters are expressed using initial Young’s modulus, initial Poisson’s ratio, maximum uniaxial compressive stress, the axial and transverse strain corresponding to the maximum uniaxial compressive stress, maximum equal biaxial compressive stress and the equal biaxial strain corresponding to maximum equal biaxial compressive stress. The ability of the proposed model to capture the response of various grade and type of concrete in ten biaxial experiments reported in the literature is examined. It is found that despite the limiting assumptions, the compression-compression, compression-tension and tension-tension stress strain response of the concrete in these ten experiments are adequately captured by the proposed model with mean R2 value of 0.94.

Classical Injective Solutions in the Large in Incompressible Nonlinear Elasticity

We consider a general class of parametrized displacement boundary value problems in incompressible nonlinear elasticity. We prove the existence of an unbounded solution branch of classical injective solutions emanating from the unforced stress-free reference configuration, which constitutes a sharp global implicit function theorem. The ramifications of this are: There exists at least one solution for any given applied loading, and or there exists solutions of arbitrarily large norm for some finite loading. We refer to this as solutions in the large. The nonlinear constraint equation enforcing incompressibility obviates the direct use of either the Leray-Schauder degree or an oriented degree previously developed for compressible problems. Instead we employ the nonlinear Fredholm degree for proper maps developed by Fitzpatrick, et. al. We establish the requisite admissibility properties in the presence of the constraint equation before proving the main results.

On incorporating osmotic prestretch/prestress in image-driven finite element simulations of cartilage.

Medical imaging performed in vivo captures geometries under Donnan osmotic loading, even when the articulating joint is otherwise mechanically unloaded. Hence patient-specific finite element (FE) models constructed from such medical images of cartilage represent osmotically induced prestretched/prestressed states. When applying classical modeling approaches to patient-specific simulations of cartilage a theoretical inconsistency arises: the in-vivo imaged geometry (used to construct the model) is not an unloaded, stress-free reference configuration. Furthermore when fitting nonlinear constitutive models that include osmotic swelling (to obtain material parameters), if one assumes that experimental data-generated from osmotically loaded cartilage-begin from a stress-free reference configuration the fitted stress-stretch relationship (parameters) obtained will actually describe a different behavior. In this study we: (1) establish a practical computational method to include osmotically induced prestretch/prestress in image-driven simulations of cartilage; and (2) investigate the influence of considering the prestretched/prestressed state both when fitting fiber-reinforced, biphasic constitutive models of cartilage that include osmotic swelling and when simulating cartilage responses. Our results highlight the importance of determining the prestretched/prestressed state within cartilage induced by osmotic loading in the imaged configuration prior to solving boundary value problems of interest. With our new constitutive model and modeling methods, we aim to improve the fidelity of FE-based, patient-specific biomechanical simulations of joints and cartilage. Improved simulations can provide medical researchers with new information often unavailable in a clinical setting, information that may contribute to better insight into the pathophysiology of cartilage diseases.

Patient-specific stress analyses in the ascending thoracic aorta using a finite-element implementation of the constrained mixture theory.

It is now a rather common approach to perform patient-specific stress analyses of arterial walls using finite-element models reconstructed from gated medical images. However, this requires to compute for every Gauss point the deformation gradient between the current configuration and a stress-free reference configuration. It is technically difficult to define such a reference configuration, and there is actually no guarantee that a stress-free configuration is physically attainable due to the presence of internal stresses in unloaded soft tissues. An alternative framework was proposed by Bellini et al. (Ann Biomed Eng 42(3):488-502, 2014). It consists of computing the deformation gradients between the current configuration and a prestressed reference configuration. We present here the first finite-element results based on this concept using the Abaqus software. The reference configuration is set arbitrarily to the in vivo average geometry of the artery, which is obtained from gated medical images and is assumed to be mechanobiologically homeostatic. For every Gauss point, the stress is split additively into the contributions of each individual load-bearing constituent of the tissue, namely elastin, collagen, smooth muscle cells. Each constituent is assigned an independent prestretch in the reference configuration, named the deposition stretch. The outstanding advantage of the present approach is that it simultaneously computes the in situ stresses existing in the reference configuration and predicts the residual stresses that occur after removing the different loadings applied onto the artery (pressure and axial load). As a proof of concept, we applied it on an ideal thick-wall cylinder and showed that the obtained results were consistent with corresponding experimental and analytical results of the well-known literature. In addition, we developed a patient-specific model of a human ascending thoracic aneurysmal aorta and demonstrated the utility in predicting the wall stress distribution in vivo under the effects of physiological pressure. Finally, we simulated the whole process preceding traditional in vitro uniaxial tensile testing of arteries, including excision from the body, radial cutting, flattening and subsequent tensile loading, showing how this process may impact the final mechanical properties derived from these in vitro tests.

An augmented iterative method for identifying a stress-free reference configuration in image-based biomechanical modeling

Continued advances in computational power and methods have enabled image-based biomechanical modeling to become an important tool in basic science, diagnostic and therapeutic medicine, and medical device design. One of the many challenges of this approach, however, is identification of a stress-free reference configuration based on in vivo images of loaded and often prestrained or residually stressed soft tissues and organs. Fortunately, iterative methods have been proposed to solve this inverse problem, among them Sellier’s method. This method is particularly appealing because it is easy to implement, convergences reasonably fast, and can be coupled to nearly any finite element package. By means of several practical examples, however, we demonstrate that in its original formulation Sellier’s method is not optimally fast and may not converge for problems with large deformations. Fortunately, we can also show that a simple, inexpensive augmentation of Sellier’s method based on Aitken’s delta-squared process can not only ensure convergence but also significantly accelerate the method.

Identifying hyperelastic and isotropic materials by examining the variation of principal direction of left Cauchy–Green deformation tensor in uniaxial loading

Many bodies of biological and engineering interest have a fibrous and layered structure. Hence, it is believed that these bodies are inhomogeneous and made up of anisotropic material. Classical mechanical experiments used to find the required material symmetry in the constitutive relation cannot distinguish inhomogeneous bodies made of isotropic material and homogeneous bodies made of anisotropic material. Therefore, it is of interest to find an alternative hypothesis so that inhomogeneity and anisotropy can be determined independent of the other. This study finds that the principal (or eigen) direction of the left Cauchy–Green deformation tensor, B does not vary with the magnitude of the applied uniaxial load at a given location whenever the body – homogeneous or inhomogeneous – is made of isotropic and hyperelastic material and the deformations are measured from a stress free reference configuration. In general, the principal direction of the left Cauchy–Green deformation tensor varies with the magnitude of the uniaxial load when the body is made up of anisotropic material. Thus, it is concluded that if the variation in the principal direction of B with the magnitude of the applied uniaxial load is experimentally investigated then one could ascertain whether the body is made up of isotropic or anisotropic material.

A simple model for myocardial changes in a failing heart

In a simplified setting, a multi-network model for remodeling in the left ventricle (LV) is developed that can mimic various pathologies of the heart. The model is an extension of the simple model introduced by Nardinocchi and Teresi [9], Nardinocchi et al. [10,11] that results in an algebraic relation for LV pressure–volume–contraction. We considered two networks, the original tissue and a new tissue, each of which has its own volume fraction, stress-free reference configuration, elastic properties, and contractility. This is used to explore the consequences of microstructural changes in the muscle tissue on LV function in terms of the pressure–volume loop during a single cardiac cycle. Special attention is paid to the stroke volume, which is directly related to cardiac output, and changes in LV wall stress caused by various disease states, including wall thinning (dilated cardiomyopathy), wall thickening (hypertrophic cardiomyopathy), contractility degradation, and stiffness changes (scarring). Various scenarios are considered that are of clinical relevance, and the extent and nature of remodeling that could lead to heart failure are identified.

Mechanical experiments to identify homogeneous bodies

All bodies are inhomogeneous at some scale but experience has shown that some of these bodies can be idealized as a homogeneous body. Here we examine which bodies can be idealized as a homogeneous body when they are subjected to a non-dissipative mechanical process. This is done by studying circumstances in which an inhomogeneous body admits pure stretch homogeneous deformations. Then, we devise experiments wherein these circumstances are prevented. If homogeneous deformation is observed in these devised experiments, the body could be modeled as a homogeneous body. We limit our analysis to a class of isotropic elastic bodies deforming from a stress free reference configuration whose Cauchy stress is explicitly related to left Cauchy–Green deformation tensor. It is further assumed that the constitutive relation is differentiable function of the position vector of material particles in the stress free reference configuration. Then, we find that a cuboid made of compressible and isotropic material could be modeled as a homogeneous body if it deforms homogeneously due to the application of the normal stresses on all of its six faces and the magnitude of the normal stresses on three orthogonal faces are different. A cuboid made of incompressible and isotropic material could be modeled as a homogeneous body, if it deforms homogeneously in two different biaxial experiments, such that the plane in which the forces are applied in the two biaxial experiments is mutually orthogonal.

Reference configurations of growing bodies

The growing bodies are considered as bodies with induced inhomogeneity caused by junction of inconsistently deformed parts. The body is formalized as an abstract smooth manifold, and all possible affine connections on it are classified. A method is shown for introducing a special connection—material connection—for which neighborhoods of all material points of a growing body pass into a stress-free state. A method for constructing a global stress-free reference configuration of a growing body as an embedding in a space with absolute parallelism is proposed. It is shown that in the case of layered accretion, the material connection corresponding to the stress-free embedding is determined by three independent functions and is in general non-Euclidean. The property of being non-Euclidean is determined by the fact that the torsion of the material connection is nonzero. We suggest to formalize the growing body as a fibration of a three-dimensional smooth manifold over a one-dimensional base, and this formalization characterizes the structure of the material connection.

Extra-fibrillar matrix mechanics of annulus fibrosus in tension and compression

The annulus fibrosus (AF) of the disk is a highly nonlinear and anisotropic material that undergoes a complex combination of loads in multiple orientations. The tensile mechanical behavior of AF in the lamellar plane is dominated by collagen fibers and has been accurately modeled using exponential functions. On the other hand, AF mechanics perpendicular to the lamella, in the radial direction, depend on the properties of the ground matrix with little to no fiber contribution. The ground matrix is mainly composed of proteoglycans (PG), which are negatively charged macromolecules that maintain the tissue hydration via osmotic pressure. The mechanical response of the ground matrix can be divided in the contribution of osmotic pressure and an elastic solid part known as extra-fibrillar matrix (EFM). Mechanical properties of the ground matrix have been measured using tensile and confined compression tests. However, EFM mechanics have not been measured directly. The objective of this study was to measure AF nonlinear mechanics of the EFM in tension and compression. To accomplish this, a combination of osmotic swelling and confined compression in disk radial direction, perpendicular to the lamella, was used. For this type of analysis, it was necessary to define a stress-free reference configuration. Thus, a brief analysis on residual stress in the disk and a procedure to estimate the reference configuration are presented. The proposed method was able to predict similar swelling deformations when using different loading protocols and models for the EFM, demonstrating its robustness. The stress-stretch curve of the EFM was linear in the range 0.9 < λ₃ < 1.3 with an aggregate modulus of 10.18±3.32 kPa; however, a significant nonlinearity was observed for compression below 0.8. The contribution of the EFM to the total aggregate modulus of the AF decreased from 70 to 30% for an applied compression of 50% of the initial thickness. The properties obtained in this study are essential for constitutive and finite element models of the AF and disk and can be applied to differentiate between functional degeneration effects such as PG loss and stiffening due to cross-linking.

Biaxial biomechanical properties of self-assembly tissue-engineered blood vessels

Along with insights into the potential for graft success, knowledge of biomechanical properties of small diameter tissue-engineered blood vessel (TEBV) will enable designers to tailor the vessels' mechanical response to closer resemble that of native tissue. Composed of two layers that closely mimic the native media and adventitia, a tissue-engineered vascular adventitia (TEVA) is wrapped around a tissue-engineered vascular media (TEVM) to produce a self-assembled tissue-engineered media/adventia (TEVMA). The current study was undertaken to characterize the biaxial biomechanical properties of TEVM, TEVA and TEVMA under physiological pressures as well as characterize the stress-free reference configuration. It was shown that the TEVA had the greatest compliance over the physiological loading range while the TEVM had the lowest compliance. As expected, compliance of the SA-TEBV fell in between with an average compliance of 2.73 MPa(-1). Data were used to identify material parameters for a microstructurally motivated constitutive model. Identified material parameters for the TEVA and TEVM provided a good fit to experimental data with an average coefficient of determination of 0.918 and 0.868, respectively. These material parameters were used to develop a two-layer predictive model for the response of a TEVMA which fit well with experimental data.