Continuing advances in the fabrication of scaffolds for tissue-engineered vascular grafts (TEVGs) are greatly expanding the scope of potential designs. Increasing recognition of the importance of local biomechanical cues for cell-mediated neotissue formation, neovessel growth, and subsequent remodeling is similarly influencing the design process. This study examines directly the potential effects of different combinations of key geometric and material properties of polymeric scaffolds on the initial mechanical state of an implanted graft into which cells are seeded or migrate. Toward this end, we developed a bilayered computational model that accounts for layer-specific thickness and stiffness as well as the potential to be residually stressed during fabrication or to swell during implantation. We found that, for realistic ranges of parameter values, the circumferential stress that would be presented to seeded or infiltrating cells is typically much lower than ideal, often by an order of magnitude. Indeed, accounting for layer-specific intrinsic swelling resulting from hydrophilicity or residual stresses not relieved via annealing revealed potentially large compressive stresses, which could lead to unintended cell phenotypes and associated maladaptive growth or, in extreme cases, graft failure. Metrics of global hemodynamics were also found to be inversely related to markers of a favorable local mechanobiological environment, suggesting a tradeoff in designs that seek mechanical homeostasis at a single scale. These findings highlight the importance of the initial mechanical state in tissue engineering scaffold design and the utility of computational modeling in reducing the experimental search space for future graft development and testing.

## Introduction

First-in-human reports demonstrate that implantable, biodegradable, polymeric scaffolds can promote neotissue formation and generate living autologous neovessels that have an ability to grow and remodel, thereby establishing tissue-engineered vascular grafts (TEVGs) as the first vascular conduits with growth potential [1]. These grafts have been used to treat congenital heart defects in children and thus have been implanted within the low-pressure caval-pulmonary circulation. Attempts to use similar scaffolds in the higher pressure arterial system have presented a greater challenge, however, with pressure-induced dilatation and/or rupture possible. There is, therefore, a pressing need for improvements in the design of biodegradable scaffolds for use as arterial replacements.

Recent studies using bilayered polymeric designs show promise as a tissue-engineered scaffold for the arterial system [24]. The basic design combines a rapidly degrading, highly porous, compliant core with a slower degrading, stiffer, protective sheath. Although such a bilayered design seems reminiscent of the native arterial structure—characterized by a compliant media surrounded by a stiff adventitia—the differential load-bearing of the two polymeric layers is designed to be very different from native. In a native large artery, the normal media bears much of the load under physiologic conditions, whereas the normal adventitia serves mainly as a sheath to protect the medial smooth muscle cells and elastic fibers from damage in cases of severe acute increases in blood pressure [5]. In contrast, the core of current bilayered TEVGs is meant to serve as a permissive environment for cellular infiltration while the outer sheath limits pressure-induced over-distension or rupture at physiological pressures. Such a design can stress-shield the core, however, where we note that stress-shielding of the media in native arteries can lead to atrophic remodeling, as in the Glagov phenomena [6]. We should thus consider carefully the potential effects of stress-shielding in developing neovessels.

In this paper, we present a parametric study of possible bilayered scaffold designs based on a nonlinear stress analysis. Mechanical stress is a fundamental design parameter in TEVGs for two primary reasons: first, vessels can fail catastrophically if stress exceeds strength at any time during the natural history of neovessel development and, second, stress is a strong stimulus for mechanobiological control of cellular function and thus extracellular matrix integrity [7]. We show that different ratios of layer thickness and material stiffness can affect the predicted stress distributions in nonintuitive ways, particularly when one includes swelling of the core or fabrication-induced residual stresses in the sheath.

In particular, hydrophilic polymers typically swell upon implantation, thus leading to a change in graft geometry and load carrying or even buckling or delamination of the scaffold in extreme cases [8]. Such swelling can also lead to circumferential compression [9,10], which may adversely affect mechanobiological outcomes especially if nascent cells are undergoing phenotypic modulation. Additionally, the fabrication process can impart residual stresses in constructs as solubilized polymers harden or fibers are woven together under tension [11]. These effects can modify the stresses experienced by both layers of the scaffold depending on the layer in which they exist and the bulk material properties. For this reason, we simulated effects of swelling of the core and residually stressing the sheath on overall stress distributions in model grafts under simulated physiologic loading.

While local mechanical stress is of paramount mechanobiological importance to neovessel formation, the implanted scaffold also interacts with native vascular and perivascular tissue, notably the adjacent vessels, flowing blood, and surrounding tissues. The material and geometry of the graft determine its structural compliance, which in turn can directly affect the surrounding tissue. Compliance mismatch has been shown to adversely affect patency rates by changing vessel hemodynamics and remodeling of the adjacent vessels [12]. Moreover, the geometry of the graft can modulate the thrombotic response of the blood by affecting the blood flow field [13]. Thus, we also consider effects of design parameters on graft compliance and geometry to better characterize the broader mechanical effects. By comparing parametric effects on these properties with effects on local metrics (e.g., transmural circumferential stress), we are better able to understand potential tradeoffs for matching homeostatic values locally or globally.

## Methods

Analysis of mechanical stress within a finitely deformed bilayered cylindrical vessel is a well-posed equilibrium problem [14], thus we simply specialize this general solution for the class of tissue-engineered constructs of interest. In particular, we describe the elastomeric core as a Hart-Smith material that can swell when hydrated, as seen for previously used core materials [8,15,16]. Swelling of an otherwise incompressible material can be modeled by imposing a constraint on the volume change once swollen [9]. Conversely, we describe the stiffer sheath as a neo-Hookean material that can be residually stressed during fabrication, as noted previously [11,17]. We model the residual stress using the concept of an opening-angle associated deformation [18]. In all the cases, we satisfy overall radial and axial equilibrium under the assumption of an axisymmetric geometry and loading, namely,
$P=∫rirI(tθθ−trr)rdr+∫rIro(tθθ−trr)rdr$
(1)
where $P$ is the distending pressure, and
$L=π∫rirI(2tzz−trr−tθθ)rdr+∫rIro(2tzz−trr−tθθ)rdr+Pπri2$
(2)
where $L$ is the applied axial load, with the core delimited by the inner and interfacial radii $r∈[ri,rI]$ and the sheath by the interfacial and outer radii $r∈(rI,ro]$. The general constitutive relation for Cauchy stress is, as usual
$t=−pI+2detFF∂W∂CFT$
(3)

where $C=FTF$ is the right Cauchy-Green tensor and $F$ is the deformation gradient tensor; $p(r)$ is a scalar Lagrange multiplier that enforces incompressibility during transient motions, which can be calculated using the radial equilibrium equation, by integrating from $ri$ to any value $r$ within the wall, and boundary conditions. Once swollen, nearly incompressible behavior is expected given the bulk stiffness of the polymers.

Specifically, let the layer-specific stored energy function $W$ be given by
$Wc=μc2∫ exp(c2(I1−3)2)dI1Ws=μs2(I1−3)$
(4)
where $Wc$ is a Hart-Smith model for the core, with $μc$ as the shear modulus and $c2$ modulating the material nonlinearity, and $I1=trC$ is the first invariant of the right Cauchy-Green tensor. $Ws$ is the neo-Hookean model for the sheath, with $μs$ as the shear modulus. For an isotropic material that undergoes swelling, the stored energy function may be modified as [9]
$W*=m(υ*)W(υ*−1/3F)$
(5)
where $υ*=detF$ is the volume change of the material due to swelling. If the material behavior further depends on $I1$ alone, we have
$W*=m(υ*)W(υ*−2/3I1)$
(6)
A convenient form for $m(υ*)$ is $υ*q$ with $q∈[0,1]$. Such an assumption (with q = 2/3) recovers Treloar's postulated form for a swelling neo-Hookean material based on statistical arguments and fits well with experimental data [19]. The Hart-Smith material can thus be modified to include swelling as
$Wc*=μcυ*2/32∫ exp(c2(υ*−2/3I1−3)2)dI1$
(7)

The deformation gradient $F$ is visualized in Fig. 1 for each of the assumed motions, which remain continuous from layer to layer (i.e., no delamination considered) during the sequence of motions: opening angle, pressurization, and swelling related.

Specifically, the deformation of the core was calculated from the prescribed swelling-induced volumetric deformation and a pressure-driven deformation according to
$Fc=[υ*Rλr*000r*R000λ]$
(8)
This swelling in the core was assumed to be homogeneous with a negligible “opening-angle induced deformation,” thus coordinates in the stress-free configuration match those in the unloaded configuration: $R=ϱ$, $Θ=φ$, and $Z=ζ$. Radial locations in the undeformed configuration are mapped to the deformed configuration via
(9)
The deformation for the sheath was calculated using the incompressibility constraint ($detFs=1)$ and a multiplicative combination of an opening angle ($Φo$) induced and pressure-induced deformation, namely,
$Fs=[(π−Φo)RπrΛλ000πr(π−Φo)R000Λλ]$
(10)
The fabrication-induced opening angle in the sheath was assumed to be homogeneous; it was also assumed that the sheath did not swell. Axial stretch $Λ$, caused by the opening-angle related deformation, was calculated from the axial equilibrium equation in the traction-free configuration. Radial locations in the sheath were found via integration to give
(11)
The equations of interest were nondimensionalized using the Buckingham-Pi method. Scales for length, mass, and time were taken as
$Ls=RD,Ms=ρRD3,Ms=ρRD3$
(12)

where $RD$ is the inner radius of the native (host) vessel at diastole (D), $ρ$ is the mass density of native tissue, and $PD$ is the native blood pressure at diastole. Table 1 lists values for these normalizing scales.

Parametric studies examined five parameters in the equilibrium solution: shear modulus of the core $μc$, undeformed thickness of the core and sheath $Hc=RI−Ri$ and $Hs=Ro−RI$, respectively, swelling ratio for the core $υ*$, and opening angle for the sheath $Φo$. Values for $μc$, $υ*$, and $Φo$ encompassed normal ranges for polymers and tissues, as reported in the literature, and allowed the examination of potentially coupled nonlinear behaviors. Ranges for thickness for each layer reflected those used in past bilayered graft experiments in mice [4]. The shear modulus of the sheath $μs$ was held at a high constant value to simulate the stiffness required to resist arterial pressure and limit strains, as in prior designs. Due to previous experimental results in rubberlike materials, $q$ was set to a single value (=2/3). Because strains in the core were limited by the stiffness of the sheath, a fixed value for the core nonlinearity parameter $c2$ was prescribed for all the simulations based on prior bilayer scaffold data. Similarly, the high stiffness of the sheath typically leads to minimal axial stretch of the graft upon implantation [20]. Hence, axial stretch was assumed to be unity for all the simulations, and effects of parametric variation on axial force are not shown. Table 1 also lists the range of values for each simulated parameter. Finally, simulations included variations in both $υ*$ and $Φo$ to illustrate their combined effects.

To examine effects of these parameters on model outputs of clinical interest, six key metrics were selected as representative of graft performance. At the functional biomechanical scale, we calculated the average transmural stress in the core, $tcavg$, the pointwise stress in the core, $tθθc$, and the pointwise stress in the sheath, $tθθs$. On the gross hemodynamics scale, we considered the diastolic luminal radius, $riD*$, and the radial compliance of the graft from diastole to systole, $Δri*/ΔP$. Also, to account for structural effects of core swelling and sheath residual stresses, the ratio of overall deformed to undeformed thickness, $h/H$, was calculated for these motions. Native values for these metrics are included or referenced in Table 1. Figure 2 schematically illustrates the input parameters and calculated metrics.

## Results

### Effects of Core Shear Modulus.

Generally, for a single layered neo-Hookean cylinder, a simple stress analysis reveals that increasing the stiffness of the material results in less distension at a fixed value of pressure and axial stretch and, thus, lower wall stress. Intuition may similarly suggest that increasing the shear modulus of the core $μc$ in a bilayered cylinder would lead to a lower mean stress within the core, but one must consider relative effects of stiffness and thickness in the core and sheath. Indeed, given a high $μs$ (fixed at a normalized value of 100), increases in $μc$ from a normalized value of 1–50 increase the share of the load borne by the core and, therefore, increase $tcavg$ (Fig. 3(a)). This behavior holds over a wide range of relative thicknesses and moduli, only reversing for very thin sheaths ($Hs/Hc<0.025$ ; not shown) and values of $μc$ approaching those of the sheath ($μc/μs>0.95$ ; not shown). Although the pressurized graft remains stretched circumferentially for all the values of modulus and thickness considered, the average circumferential stresses in the core can become negative when the relative structural stiffness of the sheath is significantly higher than that of the core. These negative values result from what is essentially a confined compression of an incompressible core since axial length was assumed to remain constant while the stiff sheath distended minimally.

Increases in sheath thickness $Hs$ decrease $tcavg$, as would be expected when increasing the overall structural stiffness of the graft. Therefore, increases in core thickness $Hc$, which would typically be thought to increase the structural stiffness of the graft, could lower the stress in both layers. Yet, the effect of $Hc$ depends on both $μc$ and $Hs$, whereby increases in $Hc$ with low $μc$ actually increase $tcavg$. This finding was supported further by the increasing value of diastolic radius $riD*$ experienced with increasing $Hc$ (Fig. 3(c)—i and ii). The opposite effect was seen for the highest value of $μc$, such that an increase in $Hc$ decreases $tcavg$ at any $Hs$. At intermediate values of $μc$, $tcavg$ increases and decreases depending on the value of $Hs$, which created transition points in Fig. 3(a)—ii and iii. While these transitions were also seen in graft compliance $Δri*/ΔP$ and diastolic radius $riD*$ (Figs. 3(b) and 3(c)—iii), the major effects of increasing $μc$ on these metrics were a marked decrease in magnitude and attenuation of thickness dependence.

### Effects of Core Swelling.

Previous work on swelling in cylindrical models focused on single layered constructs [9,10]. For a swelling core and nonswelling stiff sheath, simulations revealed lower values of $tcavg$ with increases in swelling (Fig. 4(a)). Furthermore, the compressive stresses in the core were amplified by increases in $Hs$. Depending on the degree of amplification, changes in $Hc$ raised or lowered the value of $tcavg$. At low values of $υ*$, increasing $Hc$ led to increased values of $tcavg$ if $Hs$ was low enough to allow the sheath to dilate as the core swelled. With higher values of $υ*$, increasing $Hc$ actually decreased $tcavg$ for a wide range of $Hs$. These transitions seen with $υ*$ were similar to those identified with changes in core shear modulus $μc$.

The radial compliance of the graft decreased slightly with swelling of the core (Fig. 4(b)). This suggests that distension upon pressurization is reduced by increasing the hydration of the core. Swelling of the core in the presence of a stiff or thick sheath led to a decrease in $riD*$ (Fig. 4(c)), that is, the swelling was inward, which in a clinical setting could be interpreted as stenosis. The effects of both thickness parameters on diastolic radius were qualitatively similar to effects on average core stress. Another important geometric effect arising with swelling was a progressive increase in thickness ratio $h/H$ with increasing $υ*$ (Figs. 5(a) and 5(b)). Despite this intuitive thickening, changes in $Hs$ at high fixed values of $υ*$ actually decreased $h/H$ in sufficiently thin cores (Fig. 5(b)—ii and iii). Such behavior can be ascribed to a balance between circumferential compression and thickening, where less stress is induced in the graft by distending and thinning rather than protruding further into the lumen.

### Effects of Sheath Residual Stress.

Similar to results in Sec. 3.2, increasing $Φo$ (i.e., sheath residual stress) led to dramatic decreases in $tcavg$, particularly with high values of $Hs$ (Fig. 6(a)). While residual stress in the sheath produced similar effects on $tcavg$ as did swelling in the core, a few key differences emerged. Radial compliance decreased, relative to values for both the swelling and the control cases, with increasing $Φo$ (Fig. 6(b)). Furthermore, at low values of $Hs$, Fig. 6(c) shows that $riD*$ decreased with increasing values of opening angle, whereas $riD*$ increased for this range of $Hs$ for increasing $υ*$. Opening angles also had less of an effect on the thickening of the scaffold when compared to swelling of the core, as seen by the lower magnitude of the thickness ratio for increasing opening angles (Fig. 5(c)). The pronounced decreases in $h/H$ seen with high swelling ratios and thin cores were also absent from opening-angle simulations. These effects can likely be attributed to the volume conserving nature of the assumed opening-angle induced deformation, while swelling fundamentally assumes an increase in volume.

### Coupled Effects of Core Swelling and Sheath Residual Stress.

Together, core swelling and sheath residual stresses synergistically reduce $tcavg$ to a greater degree than either effect individually. Decreases in $Δri*/ΔP$ were seen at low values of $Hs$ due to the effects of $Φo$, but little change was evident in the higher range of sheath thickness. At low values of $Hs$, the reduction in diastolic radius $riD*$ due to increasing $Φo$ was averaged with the increase of $riD*$ seen with increasing $υ*$ such that the diastolic radius changed minimally when both motions were included. For higher sheath thicknesses, $riD*$ behaved similarly to $tcavg$ and was lower than for either motion alone.

### Pointwise Stress Comparisons.

In all the cases considered thus far, it is immediately apparent that values of sheath stress $tθθs$ were much higher than values of core stress $tθθc$ (Fig. 7). These results stem from the relatively high value of sheath modulus (1 MPa) and confirm that stress-shielding is largely a function of material property imbalance, though relative values of thickness of the layers also plays a role. Increases in $μc$ raise $tθθc$ while decreasing $tθθs$. Increases in $Hs$ reduce $tθθ$ for both the core and the sheath. Yet, an increase in $Hc$ decreases $tθθc$ slightly while having a differential effect on $tθθs$ depending on the value of $μc$. With low $μc$, an increase in $Hc$ increases $tθθs$ (Fig. 7(a)—i and iii). This trend holds for an increase in $Hc$ at a higher value of $Hs$ (Fig. 7(a)—iii and iv). As typically calculated for nonresidually stressed cylinders experiencing a pressure load, stresses were slightly higher at the innermost radial location for each layer and decreased with increasing radius.

With the addition of core swelling, increasing values of $υ*$ lower $tθθc$ and raise $tθθs$ (Fig. 7(b)). The increasingly negative values of $tθθc$ highlighted a change in stress distribution when circumferential stretch becomes compressive. Swelling causes circumferential stretch to fall below 1, such that $tθθc$ values at the luminal radial location become the most negative and $tθθc$ increases toward the sheath. This effect does not change the distribution of stress in the sheath, which follows the trend from the nonswelling cases. Increasing $Hs$ lowered $tθθc$ and $tθθs$ in all the swelling cases. Increasing $Hc$ had differential effects on $tθθc$ as predicted by average stress comparisons in Fig. 4, though $tθθs$ was increased at each simulated $υ*$ value regardless of core thickness.

Changes in residual stress in the sheath resulted in similar effects as those for core swelling, that is, increases in $Φo$ generally reduce $tθθc$ and raise $tθθs$ (Fig. 7(c)). Also, changes in thickness led to similar trends in both cases. A major difference between the two motions was the change in radial dependence of $tθθs$ with increasing $Φo$ and the thickness dependence of this change. As shown in Fig. 7(c)—i, $tθθs$ tended to be highest at the layer interface in the absence of sheath residual stress, but increasing $Φo$ shifted the distribution such that the highest value of $tθθs$ was at the outer surface. For a thicker sheath, the radial stress distributions for each corresponding $Φo$ had more positive slopes in comparison to the thinner sheath cases (Fig. 7(c)—ii). Thickening of the core had the opposite effect, where the slope of the stress distribution decreased with increasing $Hc$.

## Discussion

A basic tenet of functional tissue engineering has been that one should design scaffolds so as to mimic the native tissue structure and material properties. For arteries, this implies that one should design multilayered constructs to mimic the native biomechanical properties of the medial and adventitial layers while ensuring a nonthrombogenic inner surface [21]. In this way, replacement of large (elastic) arteries should attempt to create vessels having a stress-bearing, energy-storing core (i.e., medialike structure) that is protected from acute pressure-induced over-distension by an initially less-loaded but nonlinearly stiffening, strong outer sheath (i.e., adventitialike structure). The primary energy-storing constituent within a large artery is the elastic fiber, consisting primarily of elastin (90%) and elastin-associated glycoproteins (10%) that are organized in a repeating lamellar fashion. Functional elastin is produced primarily during the perinatal period, however, its production and organization in tissue engineering remain a significant challenge [22]. The protective adventitia consists primarily of undulated type I collagen, oriented in multiple directions to endow the wall with necessary multiaxial strength while allowing normal distensibility.

The design of a tissue-engineered polymeric construct for implantation has necessarily been very different, however. The primary desire is to create a rapidly degrading porous environment that allows cells to infiltrate and synthesize native matrix as the scaffold degrades. Ideally, such synthesis and degradation should be balanced to enable the material and/or structural stiffness of the TEVG to continually mimic values of the distal and proximal host so as to not adversely disturb the (local) mechanobiology or the (global) pressure pulse wave [2325]. Toward this end, the core material is often designed to serve as the permissive porous environment while the sheath is designed to ensure structural integrity until the neotissue has formed and is capable of bearing the in vivo loads. Issues of swelling and residual stress arise from the choice of materials and fabrication methods rather than a desire to mimic biological phenomenon. Yet, the presence of residual stresses in an implanted scaffold has been proposed to offer benefits (though without detailed mechanical analysis), such as smooth muscle cell alignment in the circumferential direction and additional strength [26,27]. Given the many different design parameters—materials, porosities, symmetries, geometry, and fabrication methods—and possible nonlinear couplings, it is difficult to intuit potentially beneficial effects.

In this paper, we showed that even if a compliant core material (e.g., poly(glycerol sebacate)) is modeled by an elastomeric Hart-Smith strain energy function, its ability to distend and act as a “media” is severely limited by the stiffness of the sheath, which must currently be nonelastomeric (e.g., poly(caprolactone)) to bear the majority of the magnitude of blood pressure loading without dilatation. Preliminary simulations showed minimal sensitivity to the degree of nonlinearity in the core, as represented by the parameter $c2$. The lack of load bearing by the core contributes greatly to the low average stress values therein seen over the range of simulations. In comparison to native medial stresses, on the order of 100–200 kPa, core stresses were frequently an order of magnitude lower. Since the core is the primary site of infiltrating cells in most experimentally tested grafts, these cells may exhibit attenuated mechanostimulated matrix production while the sheath remains [2]. Given that the sheathing materials often have extended degradation times, this could both stress-shield the infiltrating cells for long periods and promote a long-term foreign body (inflammatory) response [28]. Indeed, the presence of stress-shielding suggests that much of the matrix production in the core would likely be immunomediated. Finally, when the sheathing materials eventually degrade, they often do so with a rapid loss of mechanical integrity [29]. If the core neotissue is not “preconditioned” with a significant chronic load, its capacity to suddenly bear the majority of the blood pressure induced load may be inadequate. Hence, potential stress-shielding of the core by a stiff sheath could lead to multiple problems.

Sudden rupture can arise whenever the stress experienced by the graft exceeds its strength. In the case of long-term failure noted earlier, unstressed neotissue can lack sufficient organization and crosslinking to bear normal arterial loads. Acute rupture can also occur due to a failure of the graft materials themselves shortly after implantation. The stiff polymers typically used as sheathing materials have moderate strain-at-break values, often below 20% [30]. Using a thin sheath may allow higher stresses in the core, but at an increased risk for rupture. Considering effects of either swelling in the core or fabrication-induced residual stresses in the sheath increases the stress experienced by the sheath, which could contribute to its failure. Residual stresses in the sheath were of particular importance due to the change in the stress distribution seen with thicker sheaths (Fig. 7(c)—ii and iv). In this case, the stress was highest at the outer edge of the graft. Experience in fabrication has shown that defects in the sheath, such as tears or holes, typically concentrate on the outer surface of the graft. With the highest value of stress at the outer surface of the graft, the effect of these stress concentrating defects could be amplified.

The large negative values of circumferential stress experienced by swollen and residually stressed grafts can similarly play a role in adverse outcomes. Compressive stresses have been associated with chondrocyte/osteocyte differentiation and proliferation [31]. If infiltrating cells experience compressive loads, the graft could become calcified. Local calcifications could serve as stress concentrations and promote failure; diffuse calcification could increase graft stiffness and adversely affect the mechanobiology and hemodynamics. Indeed, previous bilayered grafts have exhibited calcified regions in the core [4]. Furthermore, compressive circumferential stresses could lead to buckling of the graft. Extensive analysis of growing bilayered cylinders has shown that exceeding critical values of stress leads to collapse of the open lumen into a wrinkled geometry [32]. Such wrinkling would likely increase the thrombotic potential of the graft. In addition to thrombosis from collapse, the decrease in lumen diameter that comes with increased swelling or residual stress (cf. Figs. 4(c) and 6(c)) could contribute to stenosis or thrombosis at the proximal or distal anastomosis [13]. Finally, axial buckling from circumferentially constrained swelling can cause delamination or tearing of the scaffold, resulting in graft failure due to dissection or focal aneurysmal dilation and rupture.

Another important mechanobiological consideration is swelling alone, which increasingly thickens the graft (Fig. 5(b)). Thickening provides more wall volume for infiltrating cells to fill with inflammatory matrix [33], which may be difficult to resolve since the ability of over-thick vessels to reverse remodel appears to be minimal [34]. Implanted grafts are already typically much thicker than the native vessel to provide enough structural stiffness to bear the pressure load. A thicker graft at long time points could also affect the surrounding adjacent vessels by reducing the compliance of the construct.

Adjacent vessels anastomosed to the implanted graft are constrained to move with the ends of the construct. While area compliance of the graft can approach that of the native artery for thin sheaths and compliant cores, the strain associated with such deformations could exceed the safe limit for many stiff polymeric materials, including poly(caprolactone). Increasing the stiffness of the core seems to be a reasonable solution to the stress mismatch between the core and the sheath, but this can greatly reduce the compliance of the graft, even when the core is an order of magnitude less stiff than the sheath. Designing a successful implant will likely require a balancing act to prevent instabilities or acute rupture, to engender mechanobiological viability, and to promote normal hemodynamics. We submit that the present simulations represent another step toward the ultimate goal of identifying an optimal design.

The challenge of balance will likely need to be met by the increasing ability of material scientists to tailor microscale properties to a particular application. For example, aligning fibers circumferentially could allow increased core stresses while reducing the effect of radial compression [35]. Furthermore, residual stresses introduced into polymeric materials during fabrication can be relaxed by annealing. Such processing can also affect the crystallinity of the polymer chains, however, leading to changes in modulus and strain-at-break [36]. Understanding the microstructure and how fabrication techniques can affect the implant is critical. Due to the many possible scaffold designs, including different combinations of porosity, fiber diameter, anisotropy, and so forth, parametric studies involving more complex microstructural models may be essential in determining which designs hold promise without undue experimental burden. Herein, we used computational modeling to consider a subset of parameters and we were able to contrast over 1000 different design scenarios in a time- and cost-efficient manner.

Finally, the present study focused solely on the initial structure and properties of the scaffold, not any consideration of how these parameters evolve in vivo. Such studies can and should be pursued as well [24]. Nevertheless, we submit that “initial conditions” greatly influence subsequent growth and remodeling of an implanted TEVG and the present study provides important new insight into how a bilayered composition can be designed to interact initially with the native vasculature and infiltrating cells.

## Acknowledgment

The authors would like to thank Ramak Khosravi and Yen-Lin Wu for discussions related to this work.

## Funding Data

• National Institutes of Health (Grant Nos. R01 HL098228 and R01 HL128602).

• National Science Foundation Graduate Research Fellowship (Grant No. NSF DGE1122492).

## References

References
1.
Hibino
,
N.
,
McGillicuddy
,
E.
,
Matsumura
,
G.
,
Ichihara
,
Y.
,
Naito
,
Y.
,
Breuer
,
C.
, and
Shinoka
,
T.
,
2010
, “
Late-Term Results of Tissue-Engineered Vascular Grafts in Humans
,”
J. Thorac. Cardiovasc. Surg.
,
139
(
2
), pp.
431
436
.
2.
Wu
,
W.
,
Allen
,
R. A.
, and
Wang
,
Y.
,
2012
, “
Fast-Degrading Elastomer Enables Rapid Remodeling of a Cell-Free Synthetic Graft Into a Neoartery
,”
Nat. Med.
,
18
(
7
), pp.
1148
1153
.
3.
Yang
,
X.
,
Wei
,
J.
,
Lei
,
D.
,
Liu
,
Y.
, and
Wu
,
W.
,
2016
, “
Appropriate Density of PCL Nano-Fiber Sheath Promoted Muscular Remodeling of PGS/PCL Grafts in Arterial Circulation
,”
Biomaterials
,
88
, pp.
34
47
.
4.
Khosravi
,
R.
,
Best
,
C. A.
,
Allen
,
R. A.
,
Stowell
,
C. E.
,
Onwuka
,
E.
,
Zhuang
,
J. J.
,
Lee
,
Y. U.
,
Yi
,
T.
,
Bersi
,
M. R.
,
Shinoka
,
T.
,
Humphrey
,
J. D.
,
Wang
,
Y.
, and
Breuer
,
C. K.
,
2016
, “
Long-Term Functional Efficacy of a Novel Electrospun Poly(Glycerol Sebacate)-Based Arterial Graft in Mice
,”
Ann. Biomed. Eng.
,
44
(
8
), pp.
2402
2416
.
5.
Bellini
,
C.
,
Ferruzzi
,
J.
,
Roccabianca
,
S.
,
Martino
,
E. S. D.
, and
Humphrey
,
J. D.
,
2013
, “
A Microstructurally Motivated Model of Arterial Wall Mechanics With Mechanobiological Implications
,”
Ann. Biomed. Eng.
,
42
(
3
), pp.
488
502
.
6.
Glagov
,
S.
,
Weisenberg
,
E.
, and
Zarins
,
C. K.
,
1987
, “
Compensatory Enlargement of Human Atherosclerotic Coronary Arteries
,”
N. Engl. J. Med.
,
316
(
22
), pp.
1371
1375
.
7.
Humphrey
,
J. D.
,
2007
, “
Vascular Adaptation and Mechanical Homeostasis at Tissue, Cellular, and Sub-Cellular Levels
,”
Cell Biochem. Biophys.
,
50
(
2
), pp.
53
78
.
8.
Liu
,
Q.
,
Tian
,
M.
,
Shi
,
R.
,
Zhang
,
L.
,
Chen
,
D.
, and
Tian
,
W.
,
2007
, “
Structure and Properties of Thermoplastic Poly(Glycerol Sebacate) Elastomers Originating From Prepolymers With Different Molecular Weights
,”
J. Appl. Polym. Sci.
,
104
(
2
), pp.
1131
1137
.
9.
Demirkoparan
,
H.
, and
Pence
,
T. J.
,
2007
, “
Swelling of an Internally Pressurized Nonlinearly Elastic Tube With Fiber Reinforcing
,”
Int. J. Solids Struct.
,
44
(
11
), pp.
4009
4029
.
10.
Demirkoparan
,
H.
, and
Pence
,
T. J.
,
2007
, “
The Effect of Fiber Recruitment on the Swelling of a Pressurized Anisotropic Non-Linearly Elastic Tube
,”
Int. J. Non-Linear Mech.
,
42
(
2
), pp.
258
270
.
11.
Surrao
,
D. C.
,
Hayami
,
J. W. S.
,
Waldman
,
S. D.
, and
Amsden
,
B. G.
,
2010
, “
,”
Biomacromolecules
,
11
(
12
), pp.
3624
3629
.
12.
Abbott
,
W. M.
,
Megerman
,
J.
,
Hasson
,
J. E.
,
L'Italien
,
G.
, and
Warnock
,
D. F.
,
1987
, “
Effect of Compliance Mismatch on Vascular Graft Patency
,”
J. Vasc. Surg.
,
5
(
2
), pp.
376
382
.
13.
Binns
,
R. L.
,
Ku
,
D. N.
,
Stewart
,
M. T.
,
Ansley
,
J. P.
, and
Coyle
,
K. A.
,
1989
, “
Optimal Graft Diameter: Effect of Wall Shear Stress on Vascular Healing
,”
J. Vasc. Surg.
,
10
(
3
), pp.
326
337
.
14.
Humphrey
,
J. D.
,
2013
,
Cardiovascular Solid Mechanics: Cells, Tissues, and Organs
,
, New York.
15.
Fumio
,
U.
,
Hiroshi
,
Y.
,
Kumiko
,
N.
,
Sachihiko
,
N.
,
Kenji
,
S.
, and
Yasunori
,
M.
,
1990
, “
Swelling and Mechanical Properties of Poly(Vinyl Alcohol) Hydrogels
,”
Int. J. Pharm.
,
58
(
2
), pp.
135
142
.
16.
Chagnon
,
G.
,
Marckmann
,
G.
, and
Verron
,
E.
,
2004
, “
A Comparison of the Hart-Smith Model With Arruda-Boyce and Gent Formulations for Rubber Elasticity
,”
Rubber Chem. Technol.
,
77
(
4
), pp.
724
735
.
17.
Liu
,
W.
,
Lipner
,
J.
,
Moran
,
C. H.
,
Feng
,
L.
,
Li
,
X.
,
Thomopoulos
,
S.
, and
Xia
,
Y.
,
1986
, “
Generation of Electrospun Nanofibers With Controllable Degrees of Crimping Through a Simple, Plasticizer-Based Treatment
,”
,
27
(
16
), pp.
2583
2588
.
18.
Chuong
,
C. J.
, and
Fung
,
Y. C.
,
1986
, “
Residual Stress in Arteries
,”
Frontiers in Biomechanics
,
Springer
, New York, pp.
117
129
.
19.
Treloar
,
L. R. G.
,
1975
,
The Physics of Rubber Elasticity
,
Oxford University Press
, New York.
20.
Udelsman
,
B. V.
,
Khosravi
,
R.
,
Miller
,
K. S.
,
Dean
,
E. W.
,
Bersi
,
M. R.
,
Rocco
,
K.
,
Yi
,
T.
,
Humphrey
,
J. D.
, and
Breuer
,
C. K.
,
2014
, “
Characterization of Evolving Biomechanical Properties of Tissue Engineered Vascular Grafts in the Arterial Circulation
,”
J. Biomech.
,
47
(
9
), pp.
2070
2079
.
21.
Humphrey
,
J. D.
,
2003
,
Functional Tissue Engineering
,
Springer
, New York, pp.
35
45
.
22.
Huang
,
A. H.
,
Balestrini
,
J. L.
,
Udelsman
,
B. V.
,
Zhou
,
K. C.
,
Zhao
,
L.
,
Ferruzzi
,
J.
,
Starcher
,
B.
,
Levene
,
M. J.
,
Humphrey
,
J. D.
, and
Niklason
,
L. E.
,
2016
, “
Biaxial Stretch Improves Elastic Fiber Maturation, Collagen Arrangement, and Mechanical Properties in Engineered Arteries
,”
Tissue Eng. Part C
,
22
(
6
), pp.
524
533
.
23.
Miller
,
K. S.
,
Lee
,
Y. U.
,
Naito
,
Y.
,
Breuer
,
C. K.
, and
Humphrey
,
J. D.
,
2014
, “
Computational Model of the In Vivo Development of a Tissue Engineered Vein From an Implanted Polymeric Construct
,”
J. Biomech.
,
47
(
9
), pp.
2080
2087
.
24.
Miller
,
K. S.
,
Khosravi
,
R.
,
Breuer
,
C. K.
, and
Humphrey
,
J. D.
,
2015
, “
A Hypothesis-Driven Parametric Study of Effects of Polymeric Scaffold Properties on Tissue Engineered Neovessel Formation
,”
Acta Biomater.
,
11
, pp.
283
294
.
25.
Swillens
,
A. S.
,
Lanoye
,
L.
,
De Backer
,
J.
,
Stergiopulos
,
N.
,
Verdonck
,
P. R.
,
Vermassen
,
F.
, and
Segers
,
P.
,
2008
, “
Effect of an Abdominal Aortic Aneurysm on Wave Reflection in the Aorta
,”
IEEE Trans. Biomed. Eng.
,
55
(
5
), pp.
1602
1611
.
26.
Stitzel
,
J. D.
,
Pawlowski
,
K. J.
,
Wnek
,
G. E.
,
Simpson
,
D. G.
, and
Bowlin
,
G. L.
,
2001
, “
Arterial Smooth Muscle Cell Proliferation on a Novel Biomimicking, Biodegradable Vascular Graft Scaffold
,”
J. Biomater. Appl.
,
16
(
1
), pp.
22
33
.
27.
Teo
,
W.-E.
,
He
,
W.
, and
Ramakrishna
,
S.
,
2006
, “
Electrospun Scaffold Tailored for Tissue-Specific Extracellular Matrix
,”
Biotechnol. J.
,
1
(
9
), pp.
918
929
.
28.
Lam
,
C. X. F.
,
Hutmacher
,
D. W.
,
Schantz
,
J.-T.
,
Woodruff
,
M. A.
, and
Teoh
,
S. H.
,
2009
, “
Evaluation of Polycaprolactone Scaffold Degradation for 6 Months In Vitro and In Vivo
,”
J. Biomed. Mater. Res. A
,
90A
(
3
), pp.
906
919
.
29.
Sun
,
H.
,
Mei
,
L.
,
Song
,
C.
,
Cui
,
X.
, and
Wang
,
P.
,
2006
, “
The In Vivo Degradation, Absorption and Excretion of PCL-based Implant
,”
Biomaterials
,
27
(
9
), pp.
1735
1740
.
30.
Neppalli
,
R.
,
Maregaa
,
C.
,
Marigoa
,
A.
,
Bajgai
,
M. P.
,
Kim
,
H. Y.
, and
Causina
,
V.
,
2011
, “
Improvement of Tensile Properties and Tuning of the Biodegradation Behavior of Polycaprolactone by Addition of Electrospun Fibers
,”
Polymer
,
52
(
18
), pp.
4054
4060
.
31.
Yanagisawa
,
M.
,
Suzuki
,
N.
,
Mitsui
,
N.
,
Koyama
,
Y.
,
Otsuka
,
K.
, and
Shimizu
,
N.
,
2007
, “
Effects of Compressive Force on the Differentiation of Pluripotent Mesenchymal Cells
,”
Life Sci.
,
81
(
5
), pp.
405
412
.
32.
Moulton
,
D. E.
, and
Goriely
,
A.
,
2011
, “
Circumferential Buckling Instability of a Growing Cylindrical Tube
,”
J. Mech. Phys. Solids
,
59
(
3
), pp.
525
537
.
33.
Khosravi
,
R.
,
Miller
,
K. S.
,
Best
,
C. A.
,
Shih
,
Y. C.
,
Lee
,
Y. U.
,
Yi
,
T.
,
Shinoka
,
T.
,
Breuer
,
C. K.
, and
Humphrey
,
J. D.
,
2015
, “
Biomechanical Diversity Despite Mechanobiological Stability in Tissue Engineered Vascular Grafts Two Years Post-Implantation
,”
Tissue Eng. Part A
,
21
(
9
), pp.
1529
1538
.
34.
Sakao
,
S.
,
Tatsumi
,
K.
, and
Voelkel
,
N. F.
,
2010
, “
Reversible or Irreversible Remodeling in Pulmonary Arterial Hypertension
,”
Am. J. Respir. Cell Mol. Biol.
,
43
(
6
), pp.
629
634
.
35.
Hu
,
J.-J.
,
2015
, “
Constitutive Modeling of an Electrospun Tubular Scaffold Used for Vascular Tissue Engineering
,”
Biomech. Model. Mechanobiol.
,
14
(
4
), pp.
897
913
.
36.
Thomas
,
V.
,
Jose
,
M. V.
,
Chowdhury
,
S.
,
Sullivan
,
J. F.
,
Dean
,
D. R.
, and
Vohra
,
Y. K.
,
2006
, “
Mechano-Morphological Studies of Aligned Nanofibrous Scaffolds of Polycaprolactone Fabricated by Electrospinning
,”
J. Biomater. Sci. Polym. Ed.
,
17
(
9
), pp.
969
984
.
37.
Harrison
,
S.
,
Tamimi
,
E.
,
Uhlorn
,
J.
,
Leach
,
T.
, and
Geest
,
J. P. V.
,
2016
, “
Computationally Optimizing the Compliance of a Biopolymer Based Tissue Engineered Vascular Graft
,”
ASME J. Biomech. Eng.
,
138
(
1
), p.
014505
.