Morphing of soft tubes by anisotropic growth

We present a study of smart growth in layered cylindrical structures. We start from the characterization of a compatible growth field in an anisotropic growing tube with the aim to show a small perturbation in the compatible growth field that may produce a controlled deprivation of compatibility and localization of elastic energy storage in a composite structure made up of anisotropic growing tubes.


Introduction
Compatible growth in soft structures is challenging as it may drive large changes of shape, maintaining the structure stress-free; nevertheless, growth-induced morphings which allow to store elastic energy are noteworthy, too. Nematic elastomers are a nice example of how to design the hybrid alignment of the nematic units in an elastomeric strip so to produce, under thermal stimuli, stress-free bending without any storage of elastic energy [1,2]. On the other hand, heart mechanics delivers an amazing example of a pumping device which stores elastic energy during contraction and releases it during relaxation. Therein, anisotropic growth produces a spring-like effect which extends, dilates, and rotates the ventricle while storing elastic energy; release of elastic energy determines recovery of the initial configuration [3][4][5].
In both the situations, key issue is the control of the anisotropic growth process, of the induced morphing and of the elastic energy storage. This latter is strictly related to growth: whereas compatible growth determines stress-free morphings, a defect of compatibility produces deformations which are accompanied with elastic energy storage.
The research question we set here is about the control of the defect of compatibility of an anisotropic growth process and its effects on the storage of elastic energy occurring in the body along the growth-driven morphing. In our opinion, this is an interesting question when a layer-by-layer manufacture of a soft device is produced with each layer equipped by active fibers which can remodel under appropriate stimuli to produce a change of shape which depends on the fiber structure. On the other hand, remodeling of fibers and energy storage are two mechanisms which are strongly coupled in biological structures like human heart and chameleon tongue, just to cite some more biological structures where the active mechanism of elastic energy build-up is driven by muscular contraction corresponding to spiral fibers [6][7][8]. Actually, heart mechanics has driven our work. Indeed, it initiates from [9]. Therein, to highlight the relationships between compatibility issues and elastic energy storage, we investigated the general structure of compatible anisotropic growth fields in Fig. 1 (Left) A volume element of the body (Right). The three-states map evidencing the reference volume element (light violet), the corresponding ground volume element under growth G (violet), and the current volume element (pink). Note that the reference and the ground volume elements are both attached at the same point p; for the sake of clarity, they are depicted in two different places (color figure online) cylindrical tubes [9]. Successively, during our studies on heart mechanics [4,10,11], we discussed the role of compatibility of the growth processes which determine muscle contraction in the epicardium and endocardium layers of the human left ventricle. We conjectured that the defect of compatibility between the two layers, due to the different fiber orientations, is determinant in the storing of the elastic energy occurring during the contraction phase of the cardiac cycle when the left ventricle largely twists.
On the other end, incompatibility is likewise important during the successive release of the stored energy occurring along the relaxation phase which produces the back torsion of the left ventricle. In the present paper, we start studying layered tubes with differently fibered tube-like layers, as it holds for epicardium and endocardium layers in the heart and for the chameleon tongue. Firstly, we discuss the form of the compatible growth fields which would compete to each layer if it were free from the surroundings and of the corresponding morphings. We show how, fixed the fibers orientation, compatible growth fields can be fully explicitly characterized by two scalar parameters which are the values of the surface components of the growth fields at the inner radius of the tube. We also show how these parameters can be chosen to produce morphings with specific characteristics.
Secondly, we introduce a controlled defect of compatibility in the layered tube by gluing together tube-like layers, possibly with different fibers orientations, each of them suffering an anisotropic growth having the same structure of the compatible growth fields above explicitly characterized. The solution of the elastic problem evidences how the storage of elastic energy in the tube is localized at the interface of the tubes where also the defect of compatibility is concentrated. We present and discuss a few numerical experiments based on differently fibered tubes to highlight our results.

Background
Within the theory of elasticity, growth is introduced to model a change of the zero-stress state, or ground state, of any volume elements of a body B, identified with its reference configuration in the 3D Euclidean ambient space E. The mechanical state of that body is described by two smooth fields, the placement χ and the growth tensor G: The placement χ associates to any material point p ∈ B its position in space x = χ( p) ∈ E and determines the actual configuration χ(B) of the body. Key geometrical functions are performed by its gradient F, its adjugate F * , and its Jacobian determinant J : Let a line element a, a face a × b and a volume element a × b · c, built out of the vectors a, b, c ∈ V = VE , be attached to a place p ∈ B , as Fig. 1 (left panel) shows, then their images under the action of χ , attached to x = χ( p) , are gauged, respectively, by: On the other hand, the growth tensor G( p) associates to the line, face, and volume elements at p ∈ B their ground states (see Fig. 1): It is worth noting that these geometric elements belong to the tangent space of B at p. These two instances of 1D, 2D, and 3D body elements are in turn related by a tensor field F e such that: 3) represents the multiplicative decomposition of the deformation gradient F into the elastic deformation F e and the growth tensor G. 1 It holds: with J e = det F e . In general, a relaxed configuration cannot be realized, not even locally, and thus, it cannot be described by a placement: precisely, the symmetric tensor C o = G T G may not be the metric tensor of any possible configuration. Each neighborhood of a material point p is allowed to grow (i.e., to change its relaxed state) as the growth tensor G at p prescribes; hence, the integrity of the body may not be preserved, unless an elastic deformation F e takes place to accommodate the growth G, and realizes the final visible deformation F. Thus, the (right Cauchy-Green) strain tensor C = F T F is the metric tensor of the realized configuration χ(B).

Compatible growth fields
A special class of growth fields are the compatible growth fields which realize relaxed configurations which are different from the trivial reference configuration identified by G = I, with I the identity tensor. They are distinguished solutions of the elastic problem characterized by a minimum of energy and zero stress and are the gradient of a placement χ c .
Precisely, let us assume that the elastic strain energy density ϕ o per unit ground volume associated with the mechanical state (χ , G) be a convex functionφ o of the elastic strain measure C e = F T e F e (right Cauchy-Green strain): (2.5) The corresponding Cauchy stress tensor T is Based on previous results [9,15], we define a smooth growth field G as compatible, if there exists a placement χ c such that: If this is the case, C o is an Euclidean metric tensor, and the relaxed configuration is realizable, provided χ c satisfies any possible boundary conditions. Straightforward consequences of (2.7) are that: (i) compatible growth fields realize zero elastic strains: A compatible growth field G yields a global relaxed configuration χ c (B) of the body, and we say that the body suffers a stress-free morphing, which can be controlled by an accurate choice of the growth field G which must be compatible.

Anisotropic growth fields
We study compatible growth of incompressible cylindrical fibered tubes which twist, extend (shorten), and inflate without storing any elastic energy. Introducing a Cartesian frame {o; i 1 , i 2 , i 3 } with i i ∈ V and origin o ∈ E, the map with e R (θ ) = cos θ i 1 + sin θ i 2 the unit normal to the cylindrical surface, and e Z = i 3 the longitudinal axis of the tube, describes a tube T ⊂ E of inner radius R i , external radius R e , and axis i 3 in terms of cylindrical coordinates (R, θ, Z ) and the corresponding orthonormal polar basis (e R (θ ), e R (θ ), e Z ). To each R, it corresponds a tube-like surface S R : the inner M i and outer M e boundaries of the tube correspond to the surfaces S R i and S R e , respectively. A helicoidal fiber field in T is identified by its tangent vector field with α(R) the possibly radially varying angle measured with respect to the circumferential direction. We assume that the tube is fibered with radially uniform helicoidal fibers (α(R) = α) and introduce the anisotropic growth tensor field G as: with E = e ⊗ e, E = e * ⊗ e * , and E R = e R ⊗ e R the orientation tensors corresponding to the fields e, e * = e R × e, and e R , respectively. The anisotropic growth tensor G is completely parameterized by the eigenvectors (e(θ ), e * (θ ), and e R (θ )) and the corresponding eigenvalues (γ (R), γ ⊥ (R), and γ (R)). It is worth noting that for α = π/2 (as well as for α = 0), G describes the growth occurring along the longitudinal, circumferential, and radial direction by the fields γ (γ ⊥ ), γ ⊥ (γ ), and γ , respectively.

Anisotropic compatible growth fields
We look for stress-free morphing of the tube which admits a description in terms of extension, torsion, and inflation of the tube. In [9], we solved and discussed the general representation problem of some classes of anisotropic growth fields and characterized the placement χ c which realizes the anisotropic growth G. Before studying growth and remodeling paths, we shortly sum up the solution.
We considered placements χ c describing inflation, torsion, and extension of the tubes, with radially symmetric and longitudinally uniform gradients, keeping fixed the center of the end Z = 0 of the tube as well as hampering longitudinal displacements and rotation of that end around the vertical axis. We represent these placements as with r (R) the radial-symmetric inflation, τ the torsion, ε the axial strain of the tube, and the current cylindrical basis (a r (ϑ), a r (ϑ), a z ) defined by a r (ϑ) = cos ϑ i 1 + sin ϑ i 2 and a z = e Z = i 3 . With this, in order that an anisotropic growth G be realized by a placement χ c , the following equation must hold: It corresponds to a system of six nonlinear differential-trigonometric equations in the six unknowns: the field r (R), and the constants τ and ε characterizing χ c , and the fields γ (R), γ ⊥ (R), and γ R (R) characterizing G. From Eq. (3.14), with some manipulations, we firstly obtained [9]: and ε =γ (a 2 cos 2 α + sin 2 α) 1/2 , (3.15.1,2) withγ = γ (R i ),γ ⊥ = γ ⊥ (R i ) and the anisotropic growth factor a =γ /γ ⊥ . With this, we obtained the representation forms of the surface growth fields γ and γ ⊥ as and of the inflation r as in terms of three parameters:γ ,γ ⊥ , and α. From (3.14), we also obtained the equation γ 2 (R) = (∂r (R)/∂ R) 2 which delivers the growth along the thickness from the inflation r (R). The requests that γ 2 (R) and γ 2 ⊥ (R) be positive determine an admissibility range for a 2 : for any As both γ − (R) and γ + (R) are monotone, the limit conditions are satisfied for any R if they are and only if are satisfied for R = R e . With this, we can identify an admissible range R a of the parametersγ andγ ⊥ defined by the following inequalities Fort → 0, lim γ − = 0 and lim γ + = ∞, that is, R a covers the full parameters plane. 2 Moreover, for thin tubes, once introduced the radial dimensionless coordinate ζ such that R = R i (1 + ζ ), we can identify the fields γ and γ ⊥ with their values at ζ = 0, that is, R = R i , and write: and growth-induced thickeningγ as In thin tubes, once chosen a pair (γ ,γ ⊥ ) and the fiber orientation α, we get the thickening growthγ as well as the inflation¯ of the tube by Eqs. (3.23)-(3.24).

Growth and remodeling paths
Given (γ ,γ ⊥ ) ∈ R a , we can evaluate the extension ε and the torsion τ of the tube, the growths γ (R) and γ ⊥ (R) along the directions E and E spanning the surface S R of the tube, the inflation r (R), and the growthinduced thickening γ (R) of the tube such that the deformation process χ c is realizable. So, the pair (γ ,γ ⊥ ) is the control parameter of the morphing process. We'll follow the stress-free morphing under a continuous change of the pair (γ ,γ ⊥ ) for fixed α; and we say that we are following a growth path. Alternatively, we'll do the same under a continuous change of the fiber orientation for a fixed pair (γ ,γ ⊥ ), and we say that we are following a remodeling path [16]. So, the bisectrixγ =γ ⊥ divides the parameters plane into two portions: for anyγ <γ ⊥ (bottom portion) we get a negative (clockwise) torsion; for anyγ >γ ⊥ (top portion) we get a positive (counterclockwise) torsion, if α = 0, π/2. Following [16], we say that the bisectrix is a discriminant of perversion as torsion τ changes sign across it. Figure 2 shows tubes corresponding to different pairs (γ ⊥ , γ ) along the τ isolines, for α = π/4 and R i = 1. To each dot on the isotorsion lines, there corresponds a different morphing, depicted in the inset. In particular, the reference tube, gray dot on the dashed line, is stress-free morphed onto different configurations exhibiting thickening, inflation, and extension. Finally, the light pink-colored part of the parameters plane identifies the admissibility range R a corresponding tot = 0.1.
On the other side, the isolines corresponding to the axial strain ε = cost satisfy the equation arranged from (3.15.2). It is easy to verify that: for α = 0, isolines correspond to straight linesγ ⊥ =cost; for α = π/2, isolines correspond to straight linesγ =cost; and for α = π/4, isolines correspond to a family of hyperboles parameterized by ε (see Fig. 3). We can identify both the straight lines and the hyperboles corresponding to ε =const= 1; they always divide the plane in two portions corresponding to axial shortening and extension, and, due to it, their points are perversion points. In Fig. 3, the black dashed line corresponds to ε = 1 for α = π/4. The black dots on the isoextension lines identify the values of the parametersγ ⊥ and γ which determine the shapes represented as insets in the Figure. Each shape is fully determined by a pair (γ ⊥ ,γ ) which delivers the extension ε corresponding to the isoextension line other than torsion, inflation, surface growth, and thickening. It is worth noting that to all the shapes corresponding to dots on the bisectrix it does not correspond to any torsion as forγ ⊥ =γ , τ = 0. Comparing the shapes corresponding to the isoextension lines above and below the dashed black line allows to appreciate the lengthening opposed to the shortening of the tube. Finally, the light pink-colored part of the parameters plane identifies the admissibility range R a corresponding tot = .1.

Isochoric growth paths
In general, growth induces a change in volume. Producing a stress-free morphing by a compatible anisotropic growth field G which is also globally isochoric means keeping unchanged the volume of the tube:  Figure 4 shows the globally compatible and isochoric paths in the plane (γ ⊥ ,γ ) for thick tubes (left panel) and thin tubes (right panel). For thick tubes (t = .1), we show paths for α = 0 (blue), α = π/4 (red), and α = π/2 (purple). All these paths lie within the admissibility domain (pink shaded) in the parameters plane. When tubes are thin, it makes sense to identify a growth path which is also locally isochoric, that is, using Eqs. (3.21)-(3.24), pairs (γ ,γ ⊥ ) such thatγ or equivalently (γ 2 +γ 2 ⊥ )(γ 2 cos 2 α +γ 2 ⊥ sin 2 α) 3/2 = 2γ ⊥γ . (4.28) Equation (4.28) says that: isochoric growth paths are driven by the pairs (γ ,γ ⊥ ) such that it holds for α = 0, π/2, π/4, respectively. These latter are represented in Fig. 4 (right panel) by dashed green lines which are superposed to the corresponding globally isochoric paths evaluated fort → 0, according to the color code of the left panel. It is worth noting that for a thin tube the range of admissibility R a of the parameters covers all the plane.

Extension and torsion under remodeling paths
We can study remodeling paths corresponding to fixed values of the pair (γ ,γ ⊥ ) for different values of the fiber angle α. As observed in the previous Section, forγ =γ ⊥ , the torsion τ = 0 (see the red line in Fig. 5, left panel). Moreover, we can enhance torsion at different values of the fiber angle α switchingγ andγ ⊥ : for (γ ,γ ⊥ ) = (1, 2) we get the maximum value of τ at α < π/4 (solid purple line), whereas the same value can be attained for (γ ,γ ⊥ ) = (2, 1) at α > π/4 (dashed purple line). Likewise, switching from the first to the second pair, we switch from an extension ε which is maximum at α = 0 (solid purple line) to an extension ε which is maximum at α = π/2 (dashed purple line).

Finite elasticity with growth
As evidenced by the previous results, morphings driven by compatible growth fields are a sort of trivial universal solution of the three-dimensional elastic problem as they do not depend on the material response of the body. However, when a defect of compatibility is introduced in the modeling, material response is relevant and affects the solution of the problem. We assumed that the elastic strain energy density ϕ o associated with the mechanical state (χ , G) be a convex functionφ o of the elastic strain measure C e , and the corresponding Cauchy stress tensor is T = 2 (detF e ) −1 F e ∂φ o /∂C e F T e . Other stress measures, which are useful in discussing the elastic problem, are the reference stress S and the ground stress S o ; the relationships between them are that is, the reference stress S is the pull-back of the Cauchy stress by F as well as the pull-back of the ground stress by G . Typically, in finite elasticity, it is convenient to write down the balance equations of forces in terms of the reference stress S, to derive the constitutive representation of the stress from the elastic strain energy density in terms of the ground stress S o and to discuss the solution of the elastic problem in terms of the actual Cauchy stress T.
Assuming no bulk and boundary forces, the balance equations have the form div S = 0 in B and Sm = 0 on ∂B.
The material response of the body is described by the function ϕ o which depends only on the present value of the elastic deformation C e . As is quite common for soft materials, we also assume that the elastic deformation is isochoric, that is, J e = 1. As a consequence , the elastic strain energy determines the stress response but a constitutively indeterminate pressure term p which represents the reaction to the isochoric constraint. It holdsφ  Our isotropic choice implies that we are considering the active contribution of the fibers, whereas we neglect any reinforcing contributions of the fibers to the passive material response of the tube.

Growing bilayer tubes
From the solution of the compatibility equations, we know how to produce compatible growth/deformation fields within a tube, once the fiber main direction α and the two parametersγ andγ ⊥ have been chosen. When two of these tubes are glued and within each tube a growth pattern, corresponding to a compatible growth pattern for the tube, is assigned, we expect that the defect of compatibility between the different growth patterns allows to produce a morphing with elastic energy storage and that such storing is concentrated at the interface between the two tubes. We set up the problem within the framework above shortly reviewed of finite elasticity with growth and numerically solve a few characteristic problems to investigate our expectations. Let T a and T b be two fibered tubes which, glued together, produce a composite structure with an inner radius R i = R a equal to the inner radius of T a and an outer radius R e = R b equal to the outer radius of T b . Moreover, the outer radius of T a is equal to the inner radius of T b and denoted as R m (R i < R m < R e ). The passive material response of the two tubes is neo-Hookean, and μ a and μ b are the shear moduli in T a and T b , respectively. Both the tubes are uniformly fibered, and the helicoidal unit-vector fields e a (θ ) and e b (θ ) define the tangent fields to the fibers: e j (θ ) = cos α j e R (θ ) + sin α j e Z , j = a, b, (6.38) 3 Let us observe that it also holds: ∂Ce F T e , that is, the Cauchy stress can be viewed as the push-forward of the ground stressŜ o . with α j the angle between e j (θ ) and the circumferential direction. Finally, we assume that the growth fields and where E j = e j ⊗ e j , E j = e * j ⊗ e * j ( j = a, b) are compatible in the tubes T a and T b , respectively, that is, the fields (γ j , γ ⊥ j , γ j )( j = a, b) have the representation form discussed in Sect. 3 and are fully represented in terms of the parametersγ andγ ⊥ which take the valuesγ a andγ ⊥a for the tube T a and the valuesγ b and γ ⊥b for the tube T b .
In general, a free choice of the six parameters (γ j ,γ ⊥ j , α j ) for j = a, b would produce a compatible growth if the two tubes were not glued to, each other. In a composite structure, compatibility would be guaranteed by the choice which leaves free only three parameters over six. When fibers orientation and growth fields are not continuous at the interface between T a and T b , morphing is accompanied by elastic energy storage. The deprivation of compatibility can be controlled by releasing one by one the above conditions, as we'll show in the next Sections by a few case studies. In all the studies, we fixed the geometrical and material parameters as indicated in Table 1. The balance equation (5.33) has been solved with the finite element method, by implementing the model equations of Sect. 5 into the software COMSOL Multiphysics. We used a tetrahedral mesh with ∼ 200 k elements, having an average mesh size h=0.08 of the tube thickness, for a total of ∼ 1.1 M degrees of freedom; we adopted a mixed L2-P1 formulation, using second-order Lagrange shape functions for the displacement and linear ones for the pressure; the problem has been solved by the MUMPS direct method.

The compatible growth patterns
The first case studies correspond to a compatible growth pattern produced by an appropriate choice of the parameters (γ a ,γ ⊥a , α a ) and (γ b ,γ ⊥b , α b ) which satisfy all the three conditions (6.41). In this case, the bilayer tube behaves as a single tube: it twists, elongates, and inflates without any elastic energy expenditure with the torsion τ and the elongation ε taking the values (3.15.1,2) withγ =γ a andγ ⊥ =γ ⊥a .
Firstly, we set α a = α b = 0, that is, fibers are circumferential, longitudinal, and radial, andγ j andγ ⊥ j ( j = a, b) represent the circumferential and the longitudinal growth at R = R i ( j = a) and at R = R m ( j = b). In this case, a compatible growth pattern for the two tubes has the following representation form [see Eqs.  We numerically solve the elastic problem corresponding tō by the numerical implementation in the finite element code. As expected, the morphing corresponds to a pure dilation under zero torsion and extension, whereas we found zero stresses and zero elastic energy everywhere. Figure 6 (left panel) shows the deformed (green) tube corresponding toγ a = 1.2. The proximity between the reference (gray) tube and the green one allows to appreciate the large deformations induced by growth. Red arrows correspond to the radial displacement at the top end of the tube under the chosen compatible growth pattern, and the color code corresponds to the elastic energy which is everywhere null. Figure 6 (right panel) shows the solution corresponding to the choice α a = α b = π/2; in this case, growth produces extension combined with radial contraction of the tube which are both visible from the comparison between the gray and the green tube.
One more interesting case study corresponds to the choice α a = α b = ±π/4. The compatible surface growth pattern does not depend on the sign of the fiber angle α as Eqs. (3.15.1,2)- (3.18) show, and is the same in the two cases. Also the extension ε, which is an even function of α, is the same, whereas the torsion τ takes opposite signs in the two cases. We numerically solve the elastic problem corresponding tō The deformed configuration of the tube forγ a = 1.6 is shown in Fig. 7 (left) for α a = α b = −π/4. The right panel of the same figure shows how torsion τ and extension/shortening change forγ a taking values from 1 to 1.6. The blue line represents ε − 1 and shows the extension of the tube for growing values ofγ a (blue squares); the red line represents the torsion τ (red circles) which is negative, that is, clockwise. Atγ a = 1.6, we read on the vertical axis that ε − 1 = 0.2, that is, the final length of the tube is 1.2 times the original length. Likewise, we read τ = −0.45 which corresponds to a clockwise rotation which changes along the length of the tube and is equal to −0.45z = −0.45ε Z . Putting together the data, we conclude that the end of a tube of length L rotates of an angle equal to −0.45 · 1.2 · L. The specular situation corresponding to α a = α b = π/4 shows analogous patterns even if, in that case, torsion is positive and rotation counterclockwise.

A localized deprivation of compatibility
Our challenge is controlling the deprivation of compatibility and localizing the elastic energy storage by introducing a small perturbation in the compatible growth field above studied, that is, by relaxing just one of Fig. 7 Elastic solution corresponding to compatible growth patterns withγ = 1.6 and α a = α b = −π/4; the elastic energy is everywhere null, and torsionγ a = 1.6 is clockwise (left). Corresponding axial elongation (blue square line) and torsion (red circle line) forγ a ∈ (1, 1.6) (right) (color figure online) Fig. 8 Elastic solution corresponding to incompatible growth patterns withγ = 1.5 and α a = −α b = π/4; the elastic energy is almost everywhere null but a small area around the interface; torsion is dominated by epifibers (left). Dimensionless elastic energy density along the dimensionless radius; dashed vertical line shows the interface position (right) the three conditions (6.41). The choice we do and whose outcome we show in the following correspond to assuming α a = α b .
Actually, our choice as well as our studies on compatibility issues have been inspired by heart mechanics. Specifically, looking at the mechanics of the left ventricle, which we have been studying for ten years [4,10,11], it is well known that it is driven by blood pressure and the contraction of obliquely oriented myocardial fibers which induce a torsion of the organ about its long axis. Muscle contraction produces a sort of spring-like effect which extends, dilates, and rotates the ventricle while storing elastic energy, whereas the release of that elastic energy determines the recovery of the initial state [3,4].
So, let us choose α a = π 4 and α b = − π 4 , (6.46) and, to confirm our heart-inspiration, denote T a as the endotube and T b as the epitube. The different orientation of the fibers can be appreciated by looking at the red (endo) and blue (epi) highlighted fibers in the left panel of Fig. 8. Figure 8 shows different aspects of the solution, corresponding atγ = 1.5. Firstly, the left panel shows that torsion is clockwise, as in the left heart, due to the dominant effect of epifibers. However, the contrast offered by the endotube which would like to rotate counterclockwise determines an energy concentration around the interface between the two tubes which has been highlighted in the Figure. This energy storage aspect has been highlighted in the right panel of Fig. 8 where the dimensionless elastic energy density ϕ o /μ along the dimensionless radius has been represented; the red and blue lines correspond to the energy storage in the endotube and in the epitube, respectively. Interestingly, other than showing the peak of energy storage around the interface, it also shows an energy jump at the interface which holds even if we neglected any reinforcing contributions of the fibers to the passive material response of the tube. Indeed, the elastic energy depends also on the active material response by the elastic strain C e which is different in the two tubes, as fibers orientation is different. In this sense, the energy jump reveals the discontinuity in the active material response.
It sets also further issues, whose consideration is beyond the scope of the present work, which concern the possibility to reduce the energy peak and the energy jump by introducing different passive material response for the two tubes and/or further pills of compatibility deprivation.

Conclusions and future directions
The research question we set in the paper concerned the relationships between anisotropic growth, morphings, and storage of elastic energy in fibered soft tubes. We discussed the topic for compatible growth fields whose representation form in anisotropic fibered tubes has been explicitly characterized in [9]. Then, we investigated the effects of a deprivation of compatibility and the chances to control it to produce morphings with localized energy storage.
To do it, we glued together two tubes suffering a compatible growth to produce a composite tube which can admit a defect of compatibility at the interface between the two layers. The solution of the corresponding elastic problem showed that we produce a concentration of elastic energy at the interface between the two tubes.
The study has multiple scopes. On the one hand, it may drive the characterization of concentrated deprivation of compatibility due to growth patterns and explain morphing behaviors which can be found in Nature [8] and reproduced in synthetic structures. On the other hand, it may offer food for thought in additive manufacturing, whereas soft structures can be built layer after layer: an appropriate fiber structure of the layers can drive a controlled morphing of the structure under surface growth and small energy storage.
Moreover, after the identification of the compatible growth paths for a layered tube which allow to localize the major part of the energy stored during the growth at the interface between layers, we are asking how to drive growth in that way and to answer the question the modeling has to be set up within a dynamical context which includes the balance equations of forces and the balance equations of growth, as in [17][18][19][20].