A Ribbon Model for Nematic Polymer Networks

We present a theory of deformation of ribbons made of nematic polymer networks (NPNs). These materials exhibit properties of rubber and nematic liquid crystals, and can be activated by external stimuli of heat and light. A two-dimensional energy for a sheet of such a material has already been derived from the celebrated neo-classical energy of nematic elastomers in three space dimensions. Here, we use a dimension reduction method to obtain the appropriate energy for a ribbon from the aforementioned sheet energy. We also present an illustrative example of a rectangular NPN ribbon that undergoes in-plane serpentine deformations upon activation under an appropriate set of boundary conditions.


I. INTRODUCTION
Photoactivation of elastic materials promises to make an old dream come true: to convey energy from a distance so as to induce matter to do work with no direct, material contact with the energy source.In particular, photoactive nematic elastomers seem especially promising materials.These are rubber-like elastic solids whose constituting polymer chains incorporate photoactive molecules into nematic liquid crystal elastomers.
A typical example is the azobenzene molecule, which goes from a rod-like to a V-shaped conformation (a process called the trans-cis isomerization) when absorbing a photon of the appropriate frequency. 1Current wisdom has it that the change in shape induced by the (few) photoactive molecules present in the polymer chains (mostly constituted by nematogenic rod-like molecules unaffected by illumination) has the potential to alter the nematic orientational order reigning among the latter molecules, in a way comparable to what a change in temperature may achieve.Completely disordered molecules render the polymer network isotropic, a symmetry inherited by the macroscopic elastic response.Ordered molecules will instead grant an anisotropic response, which to first approximation is described by the neoclassical theory presented in the book by Warner and Terentjev [1] (the abundant literature that precedes it [2][3][4][5][6][7] is also of interest for the theoretician).
An interesting mechanical theory for photoactive nematic elastomers, which blends Warner and Terentjev's with the classical approach to the ordering nematic phase transition by Maier and Saupe [8], has recently been proposed in [9] (see also [10][11][12] for the relevant statistical mechanics antecedents).A vast body of reviews is available that also includes the specific topic of photoactivation of nematic elastomers [13][14][15][16][17][18], to which [19] should also be added, which is broader in scope and perhaps more germane to a mechanician's taste.General continuum theories for nematic elastomers are also available in the literature, mostly in the 3D Euclidean setting [20][21][22].As for specific applications of programmable liquid crystal elastomers, we refer the reader to a recent special issue [23].
Our development here will not be tied to a specific mechanism capable of inducing an ordering change in the nematogenic polymer chains that constitute the material.We shall be content to say that there is a scalar parameter S 0 that characterizes the order in the reference configuration and another scalar parameter, S, which characterizes the order in the current configuration, the difference between them being, however produced (either by light of by heat), the drive for the spontaneous deformation of the body.
The ordering of nematic molecules in both the reference and current configurations is further described by the director fields n 0 and n, which represent the average molecular orientation in the corresponding configurations. 2 We shall assume that the reference configuration is where the cross-linking takes place, so that both S 0 and n 0 are known at the start.S is induced by external stimuli, whose origin and nature will not play a specific role in our development.
The director n can be tied to the deformation of the body in several ways, including complete independence and complete enslaving.The nematic elastomers that we shall consider are of a special type; for them the cross-linking is so tight that the nematic director remains enslaved to the deformation. 3These material are called nematic polymer networks (NPNs). 4erhaps the most interesting manifestations of the ability of NPNs to produce changes in shape capable of doing work are achieved when they take the form of thin sheets.We represent one such sheet as a slab of thickness 2h which in its reference configuration S extends itself symmetrically on both sides of a planar surface S 0 .The director n 0 is blueprinted [26] on S 0 (in its own plane) and extended uniformly across the thickness of S, with the same scalar order parameter S 0 .
External stimuli can act on the molecular ordering so as to change S 0 into S, in a programmable way.The system is thus kicked out of equilibrium and a spontaneous deformation ensues, which makes the elastic free energy attain its minimum under the changed circumstances.
An elastic free-energy density (per unit volume), f e , was put forward in [2].This energy is delivered by a "trace formula", which was derived by assuming an anisotropic Gaussian distribution for the polymer chains in the rubber matrix. 5This formula features both the deformation f of the three-dimensional slab S and the (tensorial) measures of order in both the reference and current configurations (see Sect.II).For a sufficiently thin slab S, however, one would like to reduce f e to a function of the mapping y that changes the flat reference mid-surface S 0 into the orientable, curved surface S in the current configuration that can be regarded as the core of the deformed slab f (S), although it need not be its mid-surface.
For a NPN, for which f e ultimately depends only on f , such a dimension reduction was performed in [29] by extending a standard method of the theory of plates, known as the Kirchhoff-Love hypothesis [30].As expected, this method delivers a surface elastic free-energy density with two components, a stretching component f s scaling like h, and a bending component f b scaling like h 3 ; f s depends only on the two-dimensional stretching (or metric) tensor C := (∇y) T (∇y), while f b also depends on the invariant measures of curvature of S and the relative orientation of n in the frame of the principal directions of curvatures of S .
Not only do f s and f b scale differently with h, they also bear a different meaning in regard to the embedding of S in space.By Gauss' theorema egregium [31, p. 139], the Gaussian curvature K of S is fully determined by the metric tensor C, thus deserving the name of intrinsic curvature, whereas other measures of curvature of S , not determined by the metric, are called extrinsic and are affected by how S is embedded in space.Now, f s depends only on the intrinsic curvature, whereas f b also depends on the extrinsic ones.Studying the equilibria of S under the blended effect of both energy components f s and f b has proven a formidable task; a number of approximations have been proposed [32,33], but no satisfactory unified, general treatment of the blended surface energy has so far become available.
In this paper, we shall move yet a further step in the dimensional reduction cascade, by considering the ribbon limit for S 0 and its activated, deformed companion S .We shall describe S 0 as generated by a centerline r 0 and the imprinted director field n 0 , which in this limit reduces to a unit vector field defined on r 0 .As a consequence, S too will be described by a curve r in space, enriched with a unit vector field n, which is n 0 entrained by the deformation y of S 0 .In our model, in accord with a large body of literature (see, for example, Sect.5.2 of [34]), a NPN ribbon is a decorated curve retaining the essential geometric ingredients to represent a two-dimensional body, with one side much smaller than the other and retained only at the lowest order.This implies that the bend of the nematic director field vanishes both in the reference and deformed configurations.A different, more elaborate model of ribbon, contemplating non-vanishing nematic bend and accounting for higher-order terms in the ribbon's width, was proposed in [35] starting from a membrane Hamiltonian, which however appears to be divorced from the trace formula of three-dimensional elastomers, which is our main thrust here.
What is characteristic of our model compared to the varieties studied in classical elasticity [34,Chap. 5] is the material nature of the directors n 0 and n: they are linked to the nematic order of molecules present in the polymer network.Other studies on ribbons of activable elastomers populate the recent literature.We refer, in particular, to [36], where a ribbon model is derived via Γ-convergence from a plate energy [37] phrased in the language of non-Euclidean elasticity [38].Besides the method used, the major difference between these studies and ours lies in the assumption about the imprinted director n 0 , which is uniform throughout the cross-section of the parent sheet in our setting, whereas it is not uniform in those others, being there a possible source for kinematic incompatibility.
Our main purpose here is to derive the ribbon energy for a NPN from the plate energy obtained in [29] with no specific assumption on the ribbon's geometry.This goal is achieved in Sect.IV, after having set the necessary kinematic preliminaries in Sect.III.Section V is devoted to a special geometric setting: we study a rectangular ribbon, for which we find the explicit equilibrium planar solutions for a class of boundary conditions compatible with activable serpentine modes.We shall see how an activated serpentine shape of the ribbon is determined by the director field n 0 imprinted in its reference configuration.Finally, in Sect.VI, we collect our conclusions and comment on the implications of our work and its possible future extension.The paper is closed by three technical appendices, where details of our development are expounded for the ease of the demanding reader.

II. THIN SHEET ENERGY
Following in part [32], we recall in this section the plate-like theory for NPNs obtained in [29] through the dimension reduction method illustrated in [30]. 6This method was applied to the "trace formula" for the elastic free-energy density (per unit volume) that had been put forward for nematic elastomers [2,[40][41][42] (see also [1,Chapt. 6] for a comprehensive account of this theory) as an extension to anisotropic solids of the classical Gaussian theory for rubber elasticity (racapitulated in the landmark book [43]).
Two director fields feature in this theory, n 0 and n, the former defined in the reference configuration of the slab S and the latter defined in the current configuration f (S), where f is a diffeomorphism of S in three-dimensional Euclidean space E .In each configuration, the corresponding director represents the average orientation of the nematogenic molecules appended to the rubber polymeric matrix.They are more properly defined through the tensorial measures of material anisotropy that describe the end-to-end Gaussian distribution of polymer strands.These are the steplength tensors L 0 and L, in the reference and current configurations, respectively, which, following [7] and [44], we write as Here I is the identity (in three-dimensional space), A 0 and A are positive geometric parameters (representing the persistence length perpendicular to n 0 and n, respectively), S 0 and S are nematic scalar order parameters (related to the Maier-Saupe scalar order parameter, as shown in Appendix A).
The neo-classical theory of nematic elastomers expresses the elastic free-energy density f e (per unit volume in the reference configuration) as where F := ∇f is the deformation gradient and k > 0 is an elastic modulus.This is quite broadly known as the trace formula.
In nematic polymer networks, n is enslaved to F. In these materials, with which we are concerned in this paper, the director field n 0 is blueprinted in the elastic matrix [26] and conveyed by the deformation into n, which is thus delivered by In general, elastomers are incompressible, and so F must satisfy Both ( 3) and (4) will be enforced as constraints on all admissible deformations f of S.
With n 0 (and S 0 ) imprinted in the reference configuration at the time of cross-linking and n enslaved to the deformation, the only residual freedom lies with S, which can be changed by either thermal or optical stimuli.
It was shown in [29] that by use of ( 1) and ( 3) f e can be given the following form where C f := F T F is the right Cauchy-Green tensor associated with the deformation f and It is not difficult to show (see, for example, [32]) that F subject to ( 4) is minimized by where which shows how spontaneous deformations can be induced in these materials.For example, by heating the sample above the cross-linking temperature, we reduce the nematic order of the chains, so that S < S 0 .This in turn induces a spontaneous deformation so as to minimize the total elastic free energy: fibers along n 0 are shortened, whereas those in the plane orthogonal to n 0 are dilated.Clearly, the reverse behaviour is expected upon cooling.Thus, S can be regarded as the activation parameter of our theory, driven by external stimuli.For definiteness, we shall conventionally assume that both S 0 and S range in the interval (−1, 1). 7ere we are interested in thin sheets and in the appropriate dimension reduction of F (C f ) to be attributed to the mid surface S 0 of the slab S of thickness 2h.Formally, S 0 is a flat region in the (x 3 , x 1 ) plane of a (movable) Cartesian frame {e 1 , e 2 , e 3 } with e 2 fixed in space, and S is the set in E defined as The mapping y : S 0 → E describes the deformation of S 0 into the orientable surface S = y(S 0 ) in the deformed slab f (S); we shall assume that y is of class C 2 and that n 0 is a two-dimensional field imprinted on S 0 , so that n 0 • e 2 ≡ 0 (see Fig. 1). 8pplying (3) to the present setting, we obtain that In [29], we extended the classical Kirchhoff-Love hypothesis [30] to obtain a dimension reduction of F (C f ) in ( 6), that is, a method that convert f e in (5) into a surface energy-density fe (to be integrated over S 0 ).As standard in the theory of plates, such a surface energy is delivered by a polynomial in odd powers of h, conventionally truncated so as to retain the first two relevant ones, the first and the third power.The former is the stretching energy f s , accounting for the work done to alter distances and angles in S 0 , while the latter is the bending energy f b , accounting for the work done to fold S 0 .Thus, dropping the scaling constant 1 2 k A0 A , which has the physical dimensions of an energy per unit volume, we can write where (up to an inessential additive constant) Here C := (∇y) T (∇y) is the two-dimensional stretching (or metric) tensor, a 2 0 := n 0 • Cn 0 , H and K are the mean and Gaussian curvatures of S , defined as in terms of the (two-dimensional) curvature tensor ∇ s ν, where ν is a unit normal field to S , and The (scaled) total elastic free energy then reduces to the functional where A is the area measure.

III. RIBBON KINEMATICS
We first establish the kinematics of the deformation of a planar ribbon S 0 .We choose on S 0 a material line r 0 (s), where s is its arc-length coordinate.We call this line the directrix in consonance with the mathematical terminology of ruled surfaces.The centreline is endowed with an ordered orthonormal frame of directors {e 1 (s), e 2 , e 3 (s)} oriented such that, e 3 (s) = ∂ s r 0 (s), e 2 is a constant vector pointing out of the plane of the surface, and e 1 (s) := e 2 ×e 3 (s).Due to the orthonormality of the frame, we can associate with it a Darboux vector ω(s) such that it satisfies ∂ s e i = ω × e i , i = 1, 2, 3. Since r 0 (s) is a planar curve, ω can be fully described by just one component in the director frame, i.e. ω = ω 2 e 2 .
The planarity of the reference configuration allows us to describe it using the following parametrisation of a ruled surface,9 x(s, t) = r 0 (s) + tq 0 (s) , where q 0 (s) is a unit vector given by, and t is a coordinate along q 0 .We also define for later use, Hereafter, q 0 will be identified with the imprinted nematic director n 0 , and so it acquires a material significance, A planar strip on S 0 is represented in the (s, t) coordinates as the set {(s, t) : 0 ≤ s ≤ L, t − (s) ≤ t ≤ t + (s)}, where the functions t ± (s) are such that x(s, t ± ) corresponds to the edges of the strip (see Fig. 1).
Upon deformation, the centreline r 0 (s) ∈ S 0 is convected into a curve r(s) ∈ S .Similar to r 0 (s), we associate an ordered orthonormal director frame {d 1 (s), d 2 (s), d 3 (s)} with r(s) such that with v 3 (s) := |∂ s r| defined as the stretch of the centerline.The orthonormality of the directors d i , i = 1, 2, 3, allows us to associate with it a Darboux vector u(s) such that the bending strains about the corresponding directors.
The deformation of the reference configuration S 0 into the current configuration S is represented by the following map, y(s, t) = r(s) + ψ(s, t)q(s) , where q(s) is a unit vector given by, Equation ( 20) has an antecedent in equation ( 1) of [50], whose style and notation are often echoed here.
To obtain the deformation gradient induced by the mapping (20), we consider a curve in S 0 represented by (s(ξ), t(ξ)), parametrised by some parameter ξ.On differentiating ( 15) and ( 20) with respect to ξ, it follows that the curve and its image in the current configuration S satisfy, ẋ = (∂ s r 0 + t∂ s q 0 ) ṡ + q 0 ṫ , (22a) where the superimposed dot denotes differentiation with respect to ξ.The deformation gradient ∇y can be represented as, where vectors a and b are respectively the images of q 0 and q ⊥ 0 in the current configuration.Requiring that the identity ẏ = (∇y) ẋ be valid for all ( ṡ, ṫ), we obtain using ( 22) and ( 23), Using ∂ s r 0 = e 3 , and equations ( 16), (17), and ( 19), we give (24b) the following simpler form, With the kinematics of the deformation established, we next enforce the constraint of inextensibility on the material of the ribbon, and compute its consequences.

A. Incompressibility constraint
Nematic elastomers are subject to the constraint (4).Using ( 23), (24a), and (25), this constraint is also expressed as To solve (26) we first compute the following useful relations, where q ⊥ := d 2 × q.Making use of ( 27), we see that the most general solution of ( 26) is the following, Details of the computations leading to (28) are outlined in Appendix B. One consequence of ( 28) is the following useful relation, which we shall repeatedly invoke later.Correspondingly, we see that relations (28), upon using (29), simplify the representation for a and b from (24a) and ( 25) to the following, In light of ( 9), (18), and ( 23), we also see that (30a) implies that

B. Developability of the deformed configuration
The representation of the deformed configuration assumed in (20) is that of a ruled surface.However, since there exist ruled surfaces that are not developable, it needs to be determined whether or not (20) represents a developable surface.We show in this section that as a consequence of (28), the surface representation (20) is indeed developable.We also show that the normal to the surface coincides with the corresponding d 2 (s) along a given generator.
Whenever (26) is satisfied, the unit normal to the surface is defined as ν = a × b.Using (30), one immediately sees that ν can be written as, and so it turns out to be independent of t.
Next we show that the Gaussian curvature K = 0, confirming that the mapping under consideration results in a developable configuration.Resorting again to a generic curve (s(ξ), t(ξ)) we can write, where again a superimposed dot denotes differentiation with respect to ξ and ∇ s is the surface gradient on the deformed surface.By (22b), we see that the curvature tensor (which is symmetric) must be of the form an expression already obtained in [50] (see their equation ( 73)), implying that the Gaussian curvature of the deformed configuration vanishes identically.Inserting (34) into (33) and by use of ( 29), σ can be written as, which completes the determination of the curvature tensor.

C. Isometric limit
We now consider the limit in which the deformation map y is an isometry, 10 and compare the results with some of the standard strip models in the literature such as Wunderlich's [66] and Sadowsky's [67] models. 11The Cauchy-Green tensor of the deformation y as characterized by ( 28) can be written as, where a = |a| and b = |b|.It is then a simple matter to see, by use of (28b) and (30), that C = I 2 , where I 2 is the two-dimensional identity, if and only if in accord with equations ( 18) and ( 19) of [50], which further implies from ( 28) that α = α 0 and u 2 = ω 2 .
In this case, we can further characterize the generatrix r of the deformation y in (20).Letting κ > 0 and τ denote the curvature and torsion of r, for a = 1 we can also represent q as where T := r is the unit tangent to the directrix r.By computing q from this equation and combining the result with ( 29), we see that if ω 2 = 0 and α 0 = 0 then where N is the principal normal of r.By use of the Frenet-Serret equations, we easily convert this dichotomy into the following alternative. either where B := T × N is the binormal unit vector.Correspondingly, Combining ( 38), (40), and (41), we can write which, once inserted into (20), delivers a formula that differs only by a sign12 from equation (1.6) of [70].
In our approach, a material frame is preferred to the Frenet-Serret frame to describe the orientation in space of the ribbon, with the advantage of being also applicable at points where κ = 0.

IV. RIBBON ENERGY
In this section we carry out the further dimension reduction of the (scaled) per area energy fe of a thin sheet of NPN in ( 14) to obtain a one-dimensional energy for a ribbon.In (14), we set f s = hf 1 and f b = h 3 f 3 .Furthermore, it results from the kinematic analysis performed in the preceding section that f 1 and f 3 are given by The expression for f 3 in (43b) embodies a significant simplification of equation (11b) due to the fact that K = 0 and κ n = 0 in our current setting.Also, σ = 2H where H is the mean curvature.The expression on the right in (43a) corresponds to equation (11a).Also, it is worth recalling that in writing (43a), the unit vector q 0 has been identified with the director field n 0 imprinted on the reference configuration, as postulated in (18).While S 0 is determined at the time of cross-linking, we recall that S can be affected by external stimuli (such as light and heat); it is precisely the difference that can be induced (by external agents) between S 0 and S that drives the spontaneous deformation of a NPN sheet.We first compute the area element dA of the reference configuration S 0 which we would require for later use for dimensional reduction.The area element is given by, The (scaled) total elastic free energy F of the full surface S 0 in ( 14) can then be written as, where L is the length of the directrix r 0 in the reference configuration and f represents the (scaled) energy per unit arclength of the directrix r 0 .After some lengthy, but not difficult computations, we obtain the following representation of the reduced energy, where Details of this derivation are presented in Appendix C.Although a has a direct geometric meaning, which can be read off from (36), we find it convenient to express it via (28a) in terms of two other measures of deformation, namely, v 3 and the angle α that q makes with d 1 , In this way, equation (28b) is effectively incorporated.Since α 0 is known, by use of ( 48), for prescribed S 0 and S, F becomes a functional in the triple of functions (v 3 , α, u 1 ) subject to the boundary conditions appropriate to the specific problem at hand.For the minimizing triple, both (28c) and (28d) then also provide u 2 and u 3 , thus completing the description of the shape acquired by the actuated NPN ribbon.
Next we reduce the general energy function (46) to a simple case of rectangular geometry, and present an example where the ribbon undergoes in-plane serpentine deformations.

V. RECTANGULAR GEOMETRY
In this section, we consider a special NPN ribbon whose stress free natural configuration is a rectangle of width 2w and length L (see Fig. 2).We shall take the centerline r 0 to lie along the e 3 axis, so that s = x 3 and ω 2 = 0. We further let α 0 be a smooth function such that so that the short sides of the ribbon are both along e 1 .With such a geometric choice, Thus, the contributions to f in ( 46) that are quadratic in w vanish identically, and we obtain, where we have divided the numerator and the denominator inside the ln by cos α 0 .We briefly digress here to compare the above expression for a rectangular ribbon of finite width with that of Wunderlich's model, which considers isometric deformations.We set a = 1, v 3 = 1 and α = α 0 as per the conclusions reached in (37), along with S = S 0 = 0.Under these conditions, M = ∂ s α and V = 0 from (47d) and (47e).Substituting these in (51), and identifying tan α with η and u 2 1 = κ 2 , the energy density above reduces to, which differs from the Wunderlich's energy for a rectangular strip of finite width by an inessential additive constant.
It is then easy to see that at the lowest order in w, the above expression reduces to Sadowsky's energy for rectangular ribbons with small width. 13eturning to our main discussion, at the lowest order in w in expression (51), and using ( 48), (47d) and (47e) we obtain the following energy functional, which, in particular, is independent of ∂ s a.
The boundary conditions that we now consider will identify one special case, where the deformed ribbon is expected to be in the reference plane, while taking on a serpentine configuration.

A. Serpentining ribbon
Now, we fix the end-point at s = 0 of the ribbon and clamp the corresponding edge, so that frames {d 1 , d 2 , d 3 } and {e 1 , e 2 , e 3 } coincide for s = 0, while both the end-point at s = L and the corresponding edge are left free.Clearly, u 1 = 0 is compatible with all these boundary conditions and, furthermore, minimizes F in (53) for any α 0 , as u 1 is independent of both v 3 and α.Consequently, by (28c) and (28d), we conclude that also u 3 ≡ 0 and where we have used the relation u 2 = ∂ s ϑ, with ϑ being the angle between d 1 and e 1 .Expression (54) is to be integrated, once both α and v 3 are determined, with the initial condition ϑ(0) = 0.
Since neither ∂ s α nor ∂ s v 3 feature in F , the optimal α and v 3 are the minimizers of the integrand of F .Despite the apparent complexity of the task, finding these minimizers is boring, but not difficult.Only one such minimizing pair actually exists and is delivered explicitly by α = arctan(µ tan α 0 ) , (55a) where, for brevity, we have set It follows from ( 55) that a ≡ √ µ and which determines the deformed shape of the ribbon's centerline through the following quadratures, while its total length L * is given by As an illustration, we determine the deformed centerline of the ribbon when the imprinted nematic director q 0 is FIG. 3: Representation of α0 in (60) for a straight undeformed ribbon.The bends in the imprinted director q0 will result in bends of the ribbon's centerline upon activation.The width of the ribbons is exaggerated for clarity of presentation.
described by where n is an integer.Fig. 3 shows the decoration of the undeformed ribbon for n = 1, 2, 3, while the corresponding deformed centerline is depicted in Fig. 4 for both S > S 0 and S < S 0 ; they represent serpentining ribbons with n bends swerving on opposite sides of the undeformed ribbon.The planar deformed configurations corresponding to the three patterns of the imprinted director are shown in Fig. 4. The planar bending of the centerline is a vestige of the anisotropic stretching of the NPN sheet the ribbon was carved from.

VI. CONCLUSION
Motivated by a growing body of studies, both experimental and theoretical in nature, about the photoactivation of nematic elastomers (see, for example, [72]), we derived the elastic free-energy functional whose minimizers are the spontaneous equilibrium shapes of a NPN ribbon activated by a mismatch in molecular order induced by an external agent.Our starting point for this derivation was the plate theory that had been developed in [29] from the neo-classical formula for the energy density of nematic elastomers in three-dimensional space.
The explicit expression for the elastic energy-density (per unit length of the undeformed ribbon's centerline) is admittedly rather involved, but only because it is very general: the ribbon's centerline can be arbitrarily stretched, although the ribbon is itself an inextensible surface, and virtually no geometric limitations are imposed on the deformation of the ribbon.
The elastic free-energy functional we obtained was compared to the classical Sadowsky's and Wunderlich's functionals, valid only for isometric deformations, which represent limiting cases of deformations allowed in our model.This is where our development intersected the vast literature on the equilibria of a Möbius band, but with a major, essential difference: in our theory the straight lines that decompose the developable surfaces representing the ribbon in both reference and current configurations are material in nature, as they are determined by the imprinted nematic director.
To put our theory to the test, we considered the case of a rectangular ribbon, on which modulated director fields were imprinted at the time of cross-linking.The boundary conditions were compatible with spontaneous in-plane deformations, which we determined explicitly, showing how the serpentining shape of the deformed ribbon resonates with modulated director fields.Clearly, this is a quite simple application of our theory, where its potential of describing out-of-plane activated shapes of a ribbon is not fully deployed.To accomplish this latter task, which is presently being undertaken, we shall adapt to the present setting the theory put forward in [73], as we found intractable the Euler-Lagrange equations to the functional (53) derived from the classical formulae for the strains in terms of the representation of the director frame {d 1 , d 2 , d 3 } in a fixed, reference frame (see p. 301 of [74]).
In this paper, and in many others, the effect of an external stimulus is described in a simple, phenomenological way: we assume that it can be reduced to a uniform change in a scalar order parameter, which drives the ribbon out of equilibrium.This approach might not be sufficient to describe the details of photoactivation, which was the original inspiration of this work.Light is absorbed in solids in a way that depends on how it propagates, as described, for example, by Beer's law [75] (see also Chapt. 1 of [76]), which here becomes more stringent, as in general the activable mesogens are just a fraction of polymer strands.Implementing these details in our ribbon theory might make it applicable to real-life experiments, where the induced elastic behaviour also depends, possibly quite dramatically, on the direction of propagation of light (see, for example, [77,78]).
Upon substituting for f 1 and f 3 from (C2) and (43b), and using (35)  where V and M are defined in (47d) and (47e).A straightforward integration of the three integrals above delivers expression (46).

FIG. 1 :
FIG.1:(a)The reference configuration S0 of a planar ribbon, ruled by the directrix r0(s) and the generators in the direction of the unit vector q0(s).An orthonormal set of directors {e1(s), e2, e3(s)} are adapted to the centreline such that e3(s) is tangential to r0(s), and e2 is pointing out of the plane of the paper.The angle between q0 and e1 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 {d1(s), d2(s), d3(s)} is attached to the centreline with d3(s) oriented along the tangent vector, and d1(s) in the tangent plane at s.The angle between d1(s) and q0(s) is shown by α.

FIG. 4 :
FIG.4: Serpentining shapes of the ribbons depicted in Fig.3upon activation.Lengths are scaled to the length L of the undeformed ribbon.L * , the length of the deformed centerline, is the same for all n = 1, 2, 3.The number of bends equals n.