Bending and Stretching in a Narrow Ribbon of Nematic Polymer Networks

We study the spontaneous out-of-plane bending of a planar untwisted ribbon composed of nematic polymer networks activated by a change in temperature. Our theory accounts for both stretching and bending energies, which compete to establish equilibrium. We show that when equilibrium is attained these energy components obey a complementarity relation: one is maximum where the other is minimum. Moreover, we identify a bleaching regime: for sufficiently large values of an activation parameter (which measures the mismatch between the degrees of order in polymer organization in the reference and current configurations), the ribbon’s deformation is essentially independent of its thickness.


Introduction
A nematic elastomer differs from a standard elastomer in that it hosts nematogenic rod-like molecules within its polymer strands.If the temperature is sufficiently low, the nematogenic components develop a tendency to be aligned alike (the very signature of a nematic liquid crystal), thus introducing a degree of anisotropy in the spatial organization of polymer strands.This suffices to alter the elastic response of these materials.Since the degree of nematic order depends on temperature, a change in the latter can induce a spontaneous deformation of the body.
The nematic order is described by a scalar, here denoted Q, and a unit vector n (or director ), where the former has a bearing on the spatial organization of polymer strands and the latter represents the average orientation of the nematogenic constituents.
Nematic elastomers can be envisioned as a fluid-solid mixture, where the nematic component is fluid and the cross-linked polymer matrix is solid [1].Generally, fluid and solid components enjoy a certain degree of mutual independence.A full spectrum of materials exist, whose mechanical behavior can be classified in accordance with the degree of freedom allowed to n within the polymer matrix [2,3].
Nematic polymer networks (NPNs) are nematic elastomers of a special kind: their cross-linking is so tight that the nematic director is completely enslaved to the matrix deformation.Only these materials will form the object of our study.
We shall assume that a NPN is cross-linked in the nematic phase, which will constitute the stress-free, reference configuration of the bodies considered here.Deformations will be induced by a change in temperature, which drives the reference configuration out of equilibrium. 1 A scalar activation parameter will describe the mismatch between nematic order in the reference and current configurations; it will act as a control parameter.
We are interested in the slimiest of bodies, a ribbon.Our continuum theory is based on the three-dimensional theory of nematic elastomers revolving around the celebrated trace formula for the volume elastic energy density, which has long been studied [5][6][7][8] (as also illustrated in Sect.6.1 of [9] under the suggestive heading of neo-classical theory).
For the readers interested in broadening their background on nematic elastomers who are intimidated by the vast available literature, a recommended reference is the influential book by Warner and Terentjev [9]; the theoretical literature that preceded and prepared for it [5,[10][11][12][13][14] is also of interest.General continuum theories have also been proposed [15][16][17], some also very recently.Applications are plenty; a collection can be found in a book [18] and a recent special issue [19].Finally, some valuable guidance can be gained from the reviews [20][21][22][23][24][25][26].
This paper is an outgrowth of a previous study [27], where we derived our theory for NPN ribbons from a reduction procedure from three to two space dimensions proposed in [28] with the intent of enclosing both stretching and bending contents in the surface elastic energy of thin NPN sheets.To make our development self-contained, we recall in Sect. 2 the fundamentals of our theory; the articulation in subsections will aid the experienced reader to gather the essential information.In Sect.3, we specify the kinematics of the simple problem that we consider here.Section 4, which is the heart of the paper, is where we state and solve the variational problem for the optimal shape of a thermally activated ribbon.A plethora of activated equilibria are described in Sect. 5 for a whole range of activation and geometric parameters.Finally, in Sect.6, we collect two qualitative conclusions of our study: one identifies a regular pattern in the relative distribution of stretching and bending energies at equilibrium, the other suggests the existence of an activation regime where the ribbon's deformation is insensitive to its thickness.

Preliminaries
To make our paper self-contained, we summarize in this section some results already obtained in [28] and [27].The reader familiar with this material may jump ahead to the following sections.
Our theory is based on the trace formula for the elastic volume energy density (in the reference configuration) of nematic elastomers in the form put forward by [29], where F is the deformation gradient and the shear elastic modulus k is given by k in terms of the number density of polymer strands n s , the Boltzmann constant k B , and the absolute current temperature T .In (1), L 0 and L are the step-length tensors of the material in the reference and current configurations, respectively.They are generally different because different are the arrangements of polymer chains in the two configurations, if a change in temperature has occurred, altering the orientational order of the nematogenic molecules that comprise the polymer strands.Following [14,30], we represent them as where I is the identity (in three-dimensional space) and n 0 and n are the nematic directors in reference and current configurations.Moreover, and their twins, are expressions derived within the classical statistical mechanics model for a polymer chain thought of as consisting in freely jointed nematogenic monomers, each of length m (see Chapt. 6 of [9] and also [31,32]).In (4), Q 0 (and Q) represents the Maier-Saupe nematic scalar order parameter describing the degree of orientation of monomers in the reference (and current) configuration of the material.2Q (and, correspondingly, Q 0 ) is defined as where P 2 is the second Legendre polynomial, u ∈ S 2 is the unit vector along a single monomer, and the brackets whose end-values represent the limiting cases of u's distributed isotropically in the plane orthogonal to n and u's aligned with n, respectively.
Here, we shall assume that Q (and Q 0 ) ranges in (0, 1).In this interval, by (4b), S is a monotonically increasing function of Q, which also increases when the temperature T is decreased below the isotropic-to-nematic transition.Thus, a decrease in temperature results in an increase in S.This explains why the mechanism of thermal activation can be effectively described as having S = S 0 .A temperature increase makes S < S 0 , whereas a temperature decrease makes S > S 0 .
When, as in our case, the scalar order parameters S 0 and S are prescribed, the second term in (1) is not affected by the deformation and can be omitted, thus reducing (1) to the following bare trace formula of [5] (also discussed in [35]), which depends only on F and n and will be adopted from now on, The degree of cross-linking in the material is responsible for the mobility of the nematic director relative to the network matrix [3].A nematic polymer network (NPN) represents the tightest end of the cross-linking spectrum, where n is enslaved to the deformation. 3For these materials, which are the only ones considered here, we assume that which says that the nematic director n 0 imprinted in the polymer network at the cross-linking time is conveyed by the solid matrix of the body.We further assume that the material is incompressible,4 so that F is subject to the constraint det F = 1.
(8) By use of (3), ( 2), (4b), and ( 7), ( 6) can be given the following form [28], where C f := F T F is the right Cauchy-Green tensor associated with the deformation f and To apply our theory to the spontaneous deformation of a ribbon out of its plane (which is the theme of the following sections), we must first learn how to reduce the volume free-energy density in (9) to a surface free-energy density to be attributed to a thin sheet, thought of as a two-dimensional body.This task was accomplished in [28]; in the following subsection, we recall the main results of this study.

Dimension reduction
The method adopted in [28] is classical and consists in performing an expansion of f e retaining up to the cubic terms in the sheet's thickness, thus identifying stretching and bending contents in the surface energy by the power in the thickness they scale with.
We identify the undeformed body with a slab S of thickness 2h and midsurface S 0 in the {e 3 , e 1 } plane of a Cartesian frame.We further assume that n 0 is imprinted in S so that it does not depend on the x 2 (out-of-plane) coordinate and n 0 • e 2 ≡ 0.Moreover, we represent the three-dimensional deformation f as f (x, x 2 ) = y(x) + φ(x, x 2 )ν, (12) where x varies in S 0 , x 2 ranges in the interval [−h, h], ν is the normal to the midsurface S = y(S 0 ) of the deformed slab f (S) (see Fig. 1), and φ(x, x 2 ) is As shown in [28], the constraint of incompressibility (8) where H and K are the mean and Gaussian curvature of S , respectively, defined as H(y(x)) := 1 2 tr(∇ s ν) and K(y(x)) := det(∇ s ν), (15) in terms of the two-dimensional curvature tensor ∇ s ν at the point y(x) on S .Moreover, under assumption (7) and for h sufficiently small, the nematic director n in the present configuration f (S) is delivered by (see [4]) Since n 0 • e 2 ≡ 0, (∇y)n 0 • ν ≡ 0, and ∇ s ν is a symmetric tensor mapping the local tangent plane to S into itself, (16) shows that n • ν = 0 everywhere within f (S), but while n 0 is uniform across the thickness of S, n is not so across the thickness of f (S).However, by ( 16) and ( 13), on S n is given by and so it appears to be conveyed by the deformation of S 0 , in accordance with the three-dimensional constraint (7).Accordingly, to mimic (8), we shall assume that det(∇y) = 1, which makes S 0 an inextensible surface.Laborious calculations, building upon ( 14), (17), and (18), established in [28] that the surface elastic energy density for the deformation of S 0 into S is given by where are the (dimensionless) stretching and bending free-energy contents, respectively.In (20), C := (∇y) T (∇y) is the right Cauchy-Green tensor associated with the deformation y, a 2 0 := n 0 • Cn 0 , and which couples n with the geometry of S .As shown in detail in Appendix A, in the absence of activation, that is, for S = S 0 , f s attains its minimum subject to (18) for C = I 2 , where I 2 is the identity in two-dimensional space.Thus, in the membrane approximation, the surface elastic energy of an inactive sheet is minimized by isometries.
Finally, dropping in ( 19) the higher order terms in h, we write the total elastic free-energy functional in the form where now A denotes the area measure and e 0 h is a characteristic surface energy density precisely as e 0 in ( 10) is a characteristic volume energy density.

General theory of NPN ribbons
We reproduce here only the essential aspects of the theory of NPN ribbons, obtained through a dimension reduction procedure developed in [27].The procedure reduces the surface energy density of a sheet of NPNs, derived in [28] and summarized in Sect.2, to a line energy density valid for ribbons, whose reference and deformed configurations are assumed to be ruled surfaces.Consider a ribbon made of NPNs, whose stress-free natural configuration S 0 , and deformed configuration S are depicted in Fig. 2. We should think of Fig. 2: (a) The reference configuration S 0 of a planar ribbon, ruled by the directrix r 0 (s) and the generators in the direction of the unit vector q 0 (s).An orthonormal set of directors {e 1 (s), e 2 , e 3 (s)} are adapted to the centreline such that e 3 (s) is tangential to r 0 (s), and e 2 is pointing out of the plane of the paper.The angle between q 0 and e 1 is denoted by α 0 .(b) A current, non-planar, configuration S of the ribbon with the directrix denoted by r(s), and the direction of the generators given by q(s).An orthonormal triad {d 1 (s), d 2 (s), d 3 (s)} is attached to the centreline with d 3 (s) oriented along the tangent vector, and d 1 (s) in the tangent plane at S .The angle between d 1 (s) and q 0 (s) is shown by α. them as follows.S 0 is a stress-free configuration of the inactive ribbon (when S = S 0 ).Once activated, so that S = S 0 , S 0 is no longer stress-free and the ribbon falls out of equilibrium; a deformation then restores equilibrium by morphing S 0 into S .
We choose a material line r 0 ∈ S 0 in the natural configuration, endowed with a triad of orthonormal vectors {e 1 , e 2 , e 3 }, oriented as shown in Fig. 2.Here e 2 is independent of s, while {e 3 , e 1 } is a movable frame in the plane orthogonal to e 2 .Letting ω denote the Darboux vector associated with the frame {e 1 , e 2 , e 3 }, so that ∂ s e i = ω × e i , we easily see that ω = ω 2 e 2 .We refer to r 0 as the directrix, or the centerline, of the ribbon.We parameterize the natural configuration as a ruled surface with the following representation, where q 0 (s) is a unit vector given by q 0 (s) = cos α 0 (s)e 1 + sin α 0 (s)e 3 with Here s and t are material coordinates along the centerline r 0 and q 0 , respectively, and α 0 is a smooth function of s such that q 0 (0) and q 0 (L) coincide, respectively, with the short edges of the ribbon at s = 0 and s = L. Formally, in the (s, t) coordinates, the ribbon is represented by the set {(s, t) : 0 s L, t − (s) t t + (s)}, where t ± (s) are functions describing the two long sides of the ribbon.We will identify the unit vector q 0 with the imprinted nematic director n 0 , thereby imparting material character to the former, When the ribbon deforms, the centerline r 0 (s) and the unit vectors q 0 (s) are convected to r(s) ∈ S and q(s) ∈ S , where S denotes a deformed configuration of the ribbon.Similar to r 0 (s), the deformed centerline r(s) is endowed with an ordered orthonormal director frame {d 1 (s), d 2 (s), d 3 (s)}, where d 3 (s) is constrained to lie along the tangent (see Fig. 2).The kinematics of the centerline are captured by the following two relations, where v 3 (s) is identified as the stretch of the centerline, and u(s) is the Darboux vector associated with the frame the bending strains about the corresponding directors.
The deformed configuration S of the ribbon is parameterized as where q(s) is a unit vector given by and ψ(s, t) is a smooth function.Letting q ⊥ 0 be the unit vector in the {e 3 , e 1 } plane defined by we represent the deformation gradient ∇y as where a and b are the images of q 0 and q ⊥ 0 in the current configuration.In the present setting, (18) amounts to the condition It was proved in [27] that (30) is equivalently represented by the following relations, where a and v 3 are independent deformation measures.Equations ( 31) have important consequences (see again [27]), namely, where q ⊥ := d 2 × q, and The latter two equations express a and b in the orthonormal frame {q, q ⊥ } that spans the local tangent plane at y(x) to the current configuration of the ribbon.Moreover, in light of (33a), (29), and ( 24), (17) implies that which confers a material character to the rulings of the surface representing via (26) the current configuration of the ribbon.It also follows from the inextensibility constraint (30) that the unit normal ν to S can be written as from which we derived in [27] the following formula for the curvature tensor, where so that, in particular, the Gaussian curvature K of S vanishes identically, whereas the mean curvature H = σ/2.
In [27], using the kinematics described above along with relations (31), we obtained a dimension reduction of the sheet elastic energy in (22) to an energy defined over the centerline of a generic ribbon with a planar natural configuration.The resulting energy functional is given by where e 0 is the characteristic volume energy density introduced in (10), and

A Narrow Rectangular Ribbon
Our focus in this article is on a simple ribbon with a rectangular geometry, as shown in Fig. 3.We assume the centerline r 0 in the reference configuration to be a straight line such that it lies along e 3 , so that ω 2 ≡ 0. We further impose a condition on α 0 such that, which ensures that the shorter edge of the ribbon is along e 1 , and that the rulings (23) cover the entire material surface of the natural configuration.As a consequence of this choice, we have, To obtain an energy density for a narrow ribbon, we consider the leading order term in w in the expansion of (39), under which the energy functional takes the following form, Here F is renormalized by the scaling energy 2wL 2 e 0 .The total length L of the centerline in the reference configuration is taken to be the length scale, and consequently, all lengths from this point onward will be stated in units of L. Furthermore, f 1 and f 3 are given by, Although a has a transparent geometric meaning, defining F in (43) we preferred to express it in terms of two other (independent) measures of deformation, namely, v 3 and the angle α that n makes with d 1 .They are related to a through the equation which makes (31b) identically satisfied.Similarly (31c) can be reduced to the following Our objective here is to minimize F in (43) under constraints (46) and (31d) with an appropriate set of boundary conditions leading to out-of-plane bending of the ribbon.S 0 is the scalar order parameter frozen in the reference configuration at the time of cross-linking, while S is the current scalar order parameter induced by a change in temperature. 6Since both S and temperature T are held fixed as the spontaneous deformation unfolds, the minimizers of F are not affected by its being scaled to a quantity, e 0 , depending on both S and T .

The Variational Problem
We seek solutions with out-of-plane deformation of a narrow rectangular NPN ribbon on which the following boundary conditions are prescribed, i.e., we fix the end-to-end position vector, as well as the relative rotation of d 1 between the two ends.Using (25), constraints (47) can be given the following integral forms, We restrict our interest to the cases where the centerline deforms in the {e 2 , e 3 } plane only, meaning that the deformed centerline can be accorded the following representation, In line with this assumption of planar deformations, we further restrict the set of possible deformations with the following conditions on the strains, which makes (31d) satisfied for any α.Boundary condition (48b) is identically satisfied under this assumption, as the only non-zero component left of u is u 1 .Equations ( 50), (31c), (31d), and (41) together imply the following for non-zero u 1 , Next, we parameterize d 3 and d 2 using an angle ϑ as follows, It is then easy to deduce using the above and (25) 2 the following parametrization for the only non-zero component of the strain vector, where a prime denotes differentiation with respect to s.Using (51) and (53), the functions f 1 and f 3 in (44) reduce to, where is the effective activation parameter. 7With this definition of ν, the centerline of the ribbon would be expected to expand when ν > 1.The function f 1 , representing the stretching energy density in the ribbon, attains its minimum at v * 3 = √ ν (see also Appendix A).The coupling between bending and stretching in the ribbon under consideration is clearly delineated in the expression for f 3 , where the stretch v 3 appears in the denominator.This functional dependence of f 3 on v 3 effectively reduces the bending rigidity of the ribbon along the length in a variable fashion.
The energy functional (43) augmented with the constraint (48a) can then be written as follows in the unknown functions ϑ(s) and v 3 (s), where Λ i = Λ • e i , i ∈ {2, 3}, are the Cartesian components of a constant vector multiplier Λ employed to enforce (48a).Using standard variational procedure, we compute the variation in delivering the Euler-Lagrange (EL) equation for ϑ along with the corresponding boundary conditions, Similarly, the EL equation corresponding to v 3 is given by a straightforward partial derivative in v 3 of the integrand of (56), There are no boundary conditions associated with v 3 .

Equilibrium equations
Here we obtain a complete system of first order ODEs representing our equilibrium problem.The state space of our system is described by the vector {ϑ, u 1 , v 3 , x 2 , x 3 , Λ 2 , Λ 3 }, with 7 fields.The evolution of ϑ, u 1 , x 2 , x 3 with s can be easily obtained using (53), (58a), (25) 1 , (49), and (52b).The two multipliers Λ 2 and Λ 3 are constant, so their evolution is trivial.The only field whose evolution needs to be determined is v 3 , as the corresponding EL equation ( 59) is algebraic in v 3 .To that end, we differentiate (59) and use it to resolve for u 1 and v 3 using (58a) (and (53) to eliminate ϑ ).As a result, we arrive at the following set of first order ODEs, These 7 ODEs need to be complemented with 7 boundary conditions.Two of these are given by (58b), and the rest are the following 5, where the last condition ensures that the boundary condition on v 3 (0) is consistent with the parent algebraic equation (59) from which (60c) was derived.Our system is now complete, and can be viewed as a two point boundary value problem with ν as a control parameter.
Physical intuition suggests that for ν = 1, that is, when S = S 0 , no spontaneous deformation of the ribbon should take place, as the material is not activated.In the absence of deformation, the ribbon remains in the straight configuration, where ϑ ≡ 0, u 1 ≡ 0, and v 3 ≡ 1, which solve (60) with Λ 2 = Λ 3 = 0 and x 2 = 0, x 3 = s.It is indeed a simple matter to show that the configuration is a solution of (60) for all h > 0 and ν > 0. This is the trivial solution of our equilibrium problem and we would like to obtain solutions to (60) emanating from it.
The system of governing equations (60) is invariant under the transformations ϑ → −ϑ, Thus, for every equilibrium configuration (x 3 (s), x 2 (s)), its mirror image (x 3 (s), −x 2 (s)) is also a valid equilibrium configuration.All equilibrium configurations come in mirror-symmetric pairs.

Bifurcation analysis
Consider the following solution to (60) which differs from the trivial solution (62) by an infinitesimal amount, where ε is a small non-zero perturbation parameter.Using (58b) and (61), we deduce the following boundary conditions on the perturbations, Substituting ( 64) into (60), and equating equal powers of ε, we obtain, θ = ũ1 , (65a) x 2 = θ , (65d) The boundary conditions on x3 in (64) imply that ṽ3 vanishes (from (65c) and (65e)), and so does Λ3 (from (65g)).Using (65a) and (65b) we arrive at the following second-order linear ODE for θ, For ν < 1, the above ODE has no solutions that satisfy the boundary conditions on θ stated in (64).The implication being that all such configurations are likely to be locally stable, as no equilibrium solution emanates from it.This was to be expected, as ν < 1 implies that S > S 0 , and so the activated configuration is more ordered than the reference configuration, entailing an extension of the ribbon along d 1 (see Fig. 3), and by the surface inextensibility constraint, a contraction along d 3 , which is incompatible with boundary condition (47a). 8n the other hand, for ν ≥ 1, which is the only case that interests us here, equation (66) admits non-trivial solutions consistent with the boundary conditions in (64) when ν = ν n , where Consequently, x2 is the jet (eigenfunction) generated by x2 = sin(nπs).Expression (67) indicates that for a given branch, identified by an integer value of n, the separation between two adjacent values of ν n must increase with increasing value of h.In the next section, we quantify this effect, illuminating the interplay of two tendencies in the mechanical behavior of an activated ribbon, which for brevity we call "membrane-like" and "plate-like".The former tendency is naively expected to prevail for small values of h, where the stretching energy is dominant and the bending energy would be regarded as a perturbation; the latter, on the contrary, is expected to play a role for larger values of h.We shall see how the activation parameter can be dysfunctional towards such naive expectations.

Activated Equilibria
The governing set of equations ( 60), accompanied by the boundary conditions (58b) and (61), were solved numerically using a parameter continuation method [40], implemented in AUTO-07P [41].The bifurcation parameter was chosen to be the effective activation parameter ν in (55), which was varied from 1.0 to 2.0, and the trivial solution (62) was given as the base solution from which the continuation was carried out.We resolve the degeneracy induced by mirror symmetry by selecting a single member in each symmetric solution pair.Each equilibrium profile shown below should ideally be accompanied by its mirror image.
Fig. 4: Plots in panels (4a) to (4d) show bifurcation diagrams of the ribbon with increasing thickness.The curves are generated numerically using AUTO-07P, a continuation and bifurcation software for ODEs, while the circles denote the analytical predictions of the bifurcation points made using equation (67).The thick curve in all the plots corresponds to the first mode of deformation with one crest (see Fig. 5a).The spacing between bifurcation points increases with thickness, as predicted by equation (67).
Our main objective behind performing the forthcoming simulations is to highlight the interplay between bending and stretching of the ribbon.As we saw in the expressions for the bending and stretching energies in (54), the two mechanisms of deformation are coupled via the stretch v 3 appearing in the denominator in the expression for f 3 .The following plots reveal and quantify this non-trivial coupling between the two competing modes of deformation.
We begin by tracking the evolution of the bifurcation diagrams of the ribbon as its thickness is increased.The bifurcation plots, i.e. the L 2 norm of the function x 2 (s) vs ν, generated from the numerical simulations for increasing values of the thickness h are shown in Fig. 4. The bifurcation points (ν n , 0) as predicted by the analytical expression (67) are shown in the plots using circles, and are in good agreement with the numerical solutions.Figures 4a to  4d show the distances between the bifurcation points increasing with thickness, confirming our assertion made at the end of Sect.4.2.This increasing gap between bifurcation points is an indication of the gradual transformation of the ribbon's behavior from membrane-like to plate-like.Fig. 5: Plots in panels (5a) and (5b), respectively, correspond to the first and the second modes of deformation of a ribbon of thickness h = 0.08 (in units of L), for different values of parameter ν.The total energy F in (43) of each configuration is shown in the legend below each panel.For a given value of ν, the energy of the first mode is very close to the energy of the second mode, however, the latter remains higher than the former in all cases, and both are lower than the energy h(1 + ν 2 ) of the trivial solution (62).
A series of configurations of the ribbon centerline for different values of the effective activation parameter ν are shown in Fig. 5.The configurations in Fig. 5a lie on the first bifurcating branch (the thickened curve) in Fig. 4b, while configurations in Fig. 5b lie on the second bifurcating branch (adjacent to the thickened curve) in Fig. 4b.The total energy of each configuration is stated in a legend below the two plots in Fig. 5.A curious observation is that for a given value of ν the energies of the first and the second mode of deformation are very close, although the latter is always higher than the former.As expected, the maximum displacement is larger for the first mode than for the second mode (a rough estimate is that the former is twice as efficient as the latter).For the first mode, the maximum displacement also dramatically increases with activation, ranging from roughly 20% to above 40% of the ribbon's length, as ν ranges from 1.2 to 2.
Next we explore the effect of thickness on the maximum displacement of the ribbon in the first mode.We plot various configurations in Fig. 6a for different values of ν, ranging from 1.2 to 2.0, and thickness h ranging from 0.01 to 0.18 (in L units).We find that increasing the thickness of the ribbon reduces its maximum displacement for any given value of ν.However, as was perhaps to be expected, this effect is much more pronounced for smaller values of ν, as can be easily observed from Fig 6a .So much so that for ν = 2 the deformed shapes of the ribbon corresponding to different thicknesses are nearly indiscernible from one another.We may say that high values of the activation parameter obliterate the distinction between "membrane-like" and "plate-like" behavior of the ribbon.The dependence of displacement on thickness is quantified in Fig. 6b, which shows the maximum displacement of the ribbon centerline plotted against its thickness for several values of the effective activation parameter.
Finally, we track the distribution of the bending energy vs the stretching energy along the arc-length coordinate in Fig. 7.The gray curve represents the bending energy, while the blue curve represents the stretching energy along the arc-length coordinate.We observe that the stretching energy is always quantitatively dominant (by at least one order of magnitude) over the bending energy.However, there appears to be a sort of complementarity relation between the two.Boundary conditions dictate that the bending energy vanishes at the ribbon's end-points, whereas the stretching energy can be freely distributed at equilibrium.Now, although dominant in magnitude, the latter chooses to be concentrated where the bending energy must be absent and to be Fig.7: Several plots depicting the interplay of bending energy (f 3 h 3 ) and stretching energy (f 1 h) of a rectangular NPN ribbon for different values of h and ν plotted against s.In most panels, the two energies complement one another such that the local maximum of one corresponds to the local minimum of the other.This behavior slightly changes in panels (7i), (7j), (7k), and (7l), in that a new local maximum of the stretching energy forms in the center, where a local maximum of the bending energy lies.minimum where the bending energy is maximum.Thus, the stretching energy is primarily concentrated at the extreme ends of the ribbon, while the bending energy is concentrated in the middle.The stretching energy attains a local extremum at the center of the ribbon in all cases.This extremum is the absolute minimum for most of the plots shown.However, for plots in panels (7i), (7j), (7k), and (7l) this extremum morphs curiously into a local maximum.
On similar lines, we plot in Fig. 8 the stretch v 3 and the curvature −u 1 for various values of ν and h.We observe that the curvature and the stretch both attain their respective maxima at the center of the ribbon, with their minima located at the two ends.The coincidence of maximum curvature and maximum stretch was perhaps to be expected, as the bending rigidity in (54) decreases with increasing v 3 .We also observe that the variation in the stretch  along the length is relatively much smaller, and v 3 does not differ much from the value v * 3 = √ ν, shown by the dashed line on each plot, at which the function f 1 in (54) attains its minimum.The stretching measure v 3 is moderate (corresponding to a length increment of nearly 20%) for small ν and significant (corresponding to a length increment of nearly 40%) for large ν.
We see from Fig. 9 that, for given ν, an increase in h results into a decrease in the average stretch defined by v3 := 1 0 v 3 ds.
(68) Fig. 9: The values of the average stretch v3 , as defined by (68), are compared with the minimum v * 3 of the stretching energy for the different values of h and ν considered in our analysis.The gray datapoint " " has been added for ν = 1.2 at v3 = 1, corresponding to the trivial solution in (62), which, as shown in Fig. 4d, is the only energy minimizer for that value of ν.
Comparing this latter with v * 3 , we see that v3 ≤ v * 3 for all values of h and ν that we explored.Moreover, for given h, v3 and v * 3 get closer as ν increases and, for given ν, they become farther apart as h increases (see Fig. 9).Another telling feature is revealed by Fig. 9: for the smallest value of h, v3 falls nicely over v * 3 , marking the expected prevalence of a membrane-like behavior, which is soon lost as h is increased.
Both Fig. 9 and the plots in Figs.8j, 8k, and 8l show clearly that, for large values of ν, h has little effect on the deformation of the ribbon, in complete accordance with Fig. 6.We shall designate this as the bleaching regime.Our analysis justifies approximating v 3 with v * 3 in the bleaching regime: there the membrane-like behavior prevails irrespective of the value of h and the bending energy acts as a selection criterion that singles out a deformation out of myriads with the same minimum stretching energy.This, however, is no longer the case for smaller values of the activation parameter (say, near ν = 1), where the independent measures of deformation (v 3 , u 1 ) are closely interlaced and the mechanical behavior of the ribbon cannot be classified neither as "pure-membrane" nor as "pure-plate".

Conclusions
In this paper, which is the ideal companion of [27], we studied the outof-plane deformations of a narrow ribbon consisting of a nematic polymer network, a nematic elastomer where the degree of cross-linking is sufficiently high to justify the enslavement of the nematic director to the polymer matrix deformation.
The material is activated by a change in temperature: this produces a mismatch, measured by a control parameter ν > 1, between the degrees of order in the polymer chains organization in the reference and current configurations.This drives a spontaneous deformation of the ribbon.In our theory, both stretching and bending energies are contemplated, scaling with different powers of the ribbon's thickness h.
We deliberately kept simple the deformation pattern, choosing boundary conditions compatible with the absence of twist, as our primary aim was exploring the interplay between stretching and bending components of the elastic energy while both ν and h were independently varied.
Two major conclusions were reached here.First, stretching and bending energies appear to obey a complementarity relation: one is maximum where the other is minimum.Second, our findings point to the existence of a bleaching regime: for ν sufficiently large, the deformation of the ribbon is essentially independent of h.This latter feature may be relevant to applications where displacements induced by activation are desired to be large in not too thin ribbons.
Clearly, envisioning twist in the admissible class of deformations might have a bearing on our present conclusions, but not-we believe-to the point of obliterating them.It would be desirable to address this issue in the future.

Fig. 1 :
Fig. 1: Schematic representation of the deformation of a thin sheet.S 0 is the planar midsurface of the slab S of thickness 2h representing the reference configuration of the body; it lies in the plane {e 3 , e 1 } of Cartesian frame {e 1 , e 2 , e 3 }.S is the image under the mapping y of S 0 ; it is an oriented surface with unit normal ν.

Fig. 3 :
Fig. 3: Schematic of natural and deformed configurations of a narrow NPN ribbon of length L and width 2w.The width of the ribbon has been exaggerated for presentation purposes.

Fig. 6 :
Fig.6: Panel (6a): Configurations of a thin NPN ribbon with different values of ν and thickness h, where each color represents a fixed value of ν.In all cases, increasing thickness results in decreasing the maximum displacement of the ribbon.However, this effect appears to be more pronounced for smaller values of ν.Panel (6b): Maximum displacement of the first mode of deformation of an NPN ribbon plotted against its normalized thickness for various values of ν.

Fig. 8 :
Fig. 8: Several plots showing the distribution of the curvature −u 1 > 0 and the stretch v 3 of a rectangular NPN ribbon for different values of h and ν against s.The location of the maxima and minima of the curvature and the stretch on the centerline coincides for all cases.The values of v 3 hover close to the local minimum at v * 3 = √ ν of the stretching energy function f 1 in (54).