Abstract
We establish that, for ideal unconstrained uniaxial nematic elastomers described by a homogeneous isotropic strain-energy density function, the only smooth deformations that can be controlled by the application of surface tractions only and are universal in the sense that they are independent of the strain-energy density are those for which the deformation gradient is constant and the liquid crystal director is either aligned uniformly or oriented randomly in Cartesian coordinates. This result generalizes the classical Ericksen’s theorem for nonlinear homogeneous isotropic hyperelastic materials. While Ericksen’s theorem is directly applicable to liquid crystal elastomers in an isotropic phase where the director is oriented randomly, in a nematic phase, the constitutive strain-energy density must account also for the liquid crystal orientation which leads to significant differences in the analysis compared to the purely elastic counterpart.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.Avoid common mistakes on your manuscript.
1 Introduction
Liquid crystal elastomers (LCEs) are advanced multifunctional materials that combine elasticity with orientational order [32]. Specifically, mechanical strains give rise to changes in liquid crystalline order and, conversely, changes in the orientational order generate mechanical stresses and strains. Because of their large reversible deformations and complex material responses in the presence of external stimuli, such as heat, light, electric or magnetic fields, LCEs are suitable for a wide range of applications in science, manufacturing, and medical research [4].
LCEs can be synthesized by various methods. Without special aligning mechanisms, polydomain samples are typically obtained where the material contains multiple subdomains, each of them with their own nematic alignment, termed the “director”. Monodomain LCEs can be achieved by a secondary cross-linking after electric or magnetic fields are applied or the material is mechanically deformed to induce the desired nematic orientation.
Sitting between the world of elastomers and that of liquid crystals, the analysis of LCEs can take advantage of developments in both these fields [22]. For unconstrained isotropic hyperelastic materials, it is well known that the only controllable deformations that can be maintained universally, i.e., independently of the material parameters, by application of surface tractions only, are those characterized by a constant deformation gradient in Cartesian coordinates, the so-called homogeneous deformations. This fundamental result due to Ericksen (1955) [8] (see also [29] for an alternative proof) has been central to the phenomenological study of many elastic materials for which constitutive parameters are derived from macroscopic experimental tests. The problem of controllable deformations for incompressible homogeneous isotropic hyperelastic materials is much more involved and was examined in [7, 17, 21]. We refer also to [28] for a review of these results. Extensions of the analysis to anisotropic elasticity are presented in [34, 35] and to anelasticity in [15, 33].
In addition, nematic elastomers can sustain large reversible deformations under small applied forces [9, 14, 16, 18, 19, 26, 27, 31, 36]. For ideal LCEs, the theoretical explanation is that the energy depending isotropically on the macroscopic deformation only through the relative strain of the microstructure is minimized by these materials through a force-free state, resulting in the so-called soft-elasticity phenomenon, where the microstructure consists of many homogeneously deformed parts, known as shear striping [1, 2, 5, 6]. Strain-energy densities depending also on the macroscopic strain to realistically account for semi-soft elasticity where the applied force required is small were proposed in [11–13] (see also [23, 24]).
In this paper, we consider ideal unconstrained uniaxial nematic LCEs characterized by a homogeneous isotropic hyperelastic strain-energy density and subject to surface tractions only, i.e., in the absence of body forces. We prove that, in Cartesian coordinates, assuming that the deformation gradient is sufficiently smooth, the deformation is always homogeneous, i.e., it has a constant deformation gradient, while the liquid crystal director is either uniformly aligned or oriented randomly. If we weaken the smoothness assumption, then multiple solutions with locally constant deformation gradients can exist, such as the classical shear striping. This result is a direct extension of Ericksen’s theorem for compressible elastic materials [8] to ideal LCEs. The case where LCEs are modelled as incompressible materials [1, 20, 24] is also of considerable interest, but will not be discussed here.
2 Ideal Liquid Crystal Elastomers
An ideal uniaxial nematic LCE can be described by a homogeneous strain-energy density function of the form [22]
where \(\textbf{F}\) represents the macroscopic deformation gradient from the cross-linking state, \(\textbf{n}\) is the nematic director (a unit vector for the localized direction of uniaxial nematic alignment) in the present configuration, and \(W(\textbf{A})\) is the strain-energy density of the isotropic polymer network, depending only on the (local) elastic deformation tensor \(\textbf{A}\). The tensors \(\textbf{F}\) and \(\textbf{A}\) satisfy the following relation [32],
where
is the nematically-induced (natural) deformation tensor defining a change of frame of reference from the isotropic to a nematic phase. In equation (3), \(a>0\) is a light- or temperature-dependent shape parameter, \(\nu \) represents the optothermal analogue to the Poisson ratio relating responses in directions parallel or perpendicular to the director \(\textbf{n}\), ⊗ denotes the tensor product of two vectors, and \(\textbf{I}=\mathrm{diag}(1,1,1)\) is the identity tensor, with \(\mathrm{diag}(\cdot , \cdot , \cdot )\) denoting the diagonal second order tensor. We assume both \(a\) and \(\nu \) to be spatially independent. The tensor \(\textbf{G}_{0}\) has a similar expression to that of \(\textbf{G}\), with \(\textbf{n}_{0}\) instead of \(\textbf{n}\), \(a_{0}\) instead of \(a\), and \(\nu _{0}\) instead of \(\nu \), corresponding to the reference cross-linking state. The director \(\textbf{n}\) is an observable (spatial) quantity, and may differ from \(\textbf{n}_{0}\) by a rotation.
The square of the natural deformation tensor, known as the step-length tensor, takes the equivalent form [10] (see also [25]),
where \(c\) is the effective step length of the polymeric chain and \(\textbf{Q}\) is the symmetric traceless order parameter tensor describing orientational order in nematic liquid crystals [3, 32]. In an isotropic phase, where the liquid crystal molecules are randomly oriented, \(\textbf{Q}=\textbf{0}\).
In practice, LCEs can be modeled as incompressible materials at the cost of imposing an extra constraint, namely \(\det \textbf{F}=1\). However, here we consider universal solutions for unconstrained LCEs, i.e., incompressibility is not enforced.
In equation (1), the elastic strain-energy density \(W\) is minimized by any deformation satisfying \(\textbf{A}\textbf{A}^{T}=\textbf{I}\), and the nematic strain-energy function \(W^{(nc)}\) is minimized by any deformation satisfying \(\left (\textbf{F}\textbf{G}_{0}\right )\left (\textbf{F}\textbf{G}_{0} \right )^{T}=\textbf{G}^{2}\). Without loss of generality, we assume that the strain-energy density described by equation (1) vanishes in the reference state where \(\textbf{F}=\textbf{I}\) and \(\textbf{G}=\textbf{G}_{0}\).
By (3), if \(\textbf{R}\) is a rigid-body rotation, then the following identity holds
Since nematic elastomers are weakly cross-linked, the nematic director can rotate freely, and the material displays isotropic mechanical properties. Then, the strain-energy density given by (1) satisfies the following conditions inherited from isotropic finite elasticity [22]:
- (C1):
-
Objectivity/Frame-indifference. The constitutive equation is unaffected by a superimposed rigid-body transformation (which involves a change of position after deformation). As \(\textbf{n}\) is defined with respect to the deformed configuration, it transforms when this configuration is rotated, whereas \(\textbf{n}_{0}\) does not. Material objectivity is guaranteed by defining strain-energy functions in terms of the scalar invariants.
This is because, by the material frame indifference of \(W\),
$$ W(\textbf{R}^{\mathsf{T}}\textbf{A})=W(\textbf{A}), $$(6)and by (2),
$$ \textbf{R}^{\mathsf{T}}\textbf{F}=\left (\textbf{R}^{\mathsf{T}} \textbf{G}\textbf{R}\right )\left (\textbf{R}^{\mathsf{T}}\textbf{A} \right )\textbf{G}_{0}^{-1}. $$(7)Then (1), (5), (6) and (7) imply
$$ W^{(nc)}(\textbf{R}^{\mathsf{T}}\textbf{F},\textbf{R}^{\mathsf{T}} \textbf{n})=W(\textbf{R}^{\mathsf{T}}\textbf{A})=W(\textbf{A})=W^{(nc)}( \textbf{F},\textbf{n}). $$(8) - (C2):
-
Isotropy. The constitutive equation is unaffected by a rigid-body transformation prior to deformation. As \(\textbf{n}\) is defined with respect to the deformed configuration, it does not change when the reference configuration is rotated, whereas \(\textbf{n}_{0}\) does. For isotropic materials, the strain-energy function is a symmetric function of the principal stretch ratios.
This is because, since \(W\) is isotropic, i.e.,
$$ W(\textbf{A})=W(\textbf{A}\textbf{R}), $$(9)and (2) holds, it follows that
$$ \textbf{F}\textbf{R}=\textbf{G}\left (\textbf{A}\textbf{R}\right )\left ( \textbf{R}^{\mathsf{T}}\textbf{G}_{0}^{-1}\textbf{R}\right ). $$(10)$$ W^{(nc)}(\textbf{F}\textbf{R},\textbf{n})=W(\textbf{A}\textbf{R})=W( \textbf{A})=W^{(nc)}(\textbf{F},\textbf{n}). $$(11)
Under the frame-indifference condition (C1), the LCE model defined by (1) can be expressed equivalently in terms of the scalar invariants [30], as follows:
where
and \(\{I_{1},I_{2},I_{3}\}\) are the principal invariants of the elastic Cauchy-Green tensors \(\textbf{A}\textbf{A}^{\mathsf{T}}\) and \(\textbf{A}^{\mathsf{T}}\textbf{A}\).
The following relations between \(\left \{I^{(nc)}_{1},I^{(nc)}_{2},I^{(nc)}_{3},I^{(nc)}_{4},I^{(nc)}_{5} \right \}\) and \(\{I_{1},I_{2},I_{3}\}\) are obtained:
where \(\det \textbf{G}=a^{(1-2\nu )/3}\),
and, using (18),
Note that the above expression for \(I_{2}\) holds since, by the Cayley-Hamilton theorem, we have
and multiplying the above equation by \(\left (\textbf{F}\textbf{G}_{0}^{2}\textbf{F}^{\mathsf{T}}\right )^{-1}\) gives
Hence,
By the isotropy condition (C2), the LCE model given by (1) can be written equivalently as
where \(\lambda _{1}^{2}\), \(\lambda _{2}^{2}\), \(\lambda _{3}^{2}\) are the principal eigenvalues of \(\textbf{F}\textbf{G}_{0}^{2}\textbf{F}^{\mathsf{T}}\) and \(\alpha _{1}^{2}\), \(\alpha _{2}^{2}\), \(\alpha _{3}^{2}\) are the principal eigenvalues of \(\textbf{A}^{\mathsf{T}}\textbf{A}\). The following relations between these principal eigenvalues and the corresponding principal invariants hold,
and
as usual.
3 Stresses and Stress-Free Configurations
For nematic LCEs, the director is ‘free’ to rotate, hence \(\textbf{F}\) and \(\textbf{n}\) are independent variables. Since \(\textbf{G}\) and \(\textbf{G}_{0}\) are symmetric, the Cauchy stress tensor for the nematic material with the strain-energy function described by (1) is calculated as follows [22],
where \(J=\det \textbf{F}=\det \textbf{G}\det \textbf{A}\det \textbf{G}_{0}^{-1} \) and
is the Cauchy stress tensor from isotropic finite elasticity. Equivalently,
where
are scalar functions of the principal invariants \(\{I_{1}, I_{2}, I_{2}\}\) of \(\textbf{B}=\textbf{A}\textbf{A}^{T}\).
The first Piola-Kirchhoff stress tensor for the nematic material is equal to
where \(\mathrm{Cof}\ \textbf{F}=J\textbf{F}^{-T}\) and
is the first Piola-Kirchhoff stress tensor from isotropic elasticity.
The corresponding second Piola-Kirchhoff stress tensor is
where \(\textbf{S}=\textbf{A}^{-1}\textbf{P}\) is the second Piola-Kirchhoff from isotropic elasticity.
Since the relation (29) between the Cauchy stress tensor \(\textbf{T}\) and the Cauchy-Green tensor \(\textbf{B}\) is not invertible in general, more than one deformation gradient \(\textbf{A}\) may induce the same stress tensor \(\textbf{T}\). Thus the Cauchy stress tensor \(\textbf{T}^{(nc)}\) in a nematic LCE may also be generated by more than one deformation gradient \(\textbf{F}\). Such deformations corresponding to the same stress tensor can also alternate in LCE materials to produce inhomogeneous patterns. In this case, for geometric compatibility, any two different gradient tensors \(\textbf{F}\) and \(\widehat{\textbf{F}}\) corresponding to two alternating phases must be rank-one connected, i.e., \(\mathrm{rank}\left (\textbf{F}-\widehat{\textbf{F}}\right )=1\).
For example, assuming that \(\textbf{T}=\textbf{0}\) for rigid-body rotations, i.e., when \(\textbf{A}=\textbf{R}\), then \(\textbf{T}^{(nc)}=\textbf{0}\) when \(\textbf{F}=\textbf{G}\textbf{R}\textbf{G}_{0}^{-1}\). In particular, in the absence of elastic deformations, such that \(\textbf{A}=\textbf{I}\) and \(\textbf{G}=\textbf{G}_{0}\), the undeformed configuration, with \(\textbf{F}=\textbf{I}\), is stress free, i.e., \(\textbf{T}^{(nc)}=\textbf{0}\). However, since it is possible for a stress-free state to be generated by more than one deformation, when these deformations are geometrically compatible, they can also alternate in the same material, producing an inhomogeneous pattern.
More generally, if \(\textbf{T}=p\textbf{I}\), where \(p\) is constant, while \(\textbf{A}=c \textbf{R}\), where \(c\) is constant and \(\textbf{R}\) is a rigid-body rotation, then \(\textbf{T}^{(nc)}=p\left (\det \textbf{G}^{-1}\det \textbf{G}_{0} \right )\textbf{I}\) for \(\textbf{F}=c \textbf{G}\textbf{R}\textbf{G}_{0}^{-1}\). As it is possible for the same stress to be generated by different deformations, if these are geometrically compatible, then they can also alternate in the same material to produce inhomogeneous patterns.
Our result concerning controllable deformation in ideal nematic elastomers is established in the next section.
4 Controllable Deformations
Theorem 1
For an ideal unconstrained uniaxial nematic LCE characterized by equation (1), such that \(\textbf{G}_{0}\) is constant in Cartesian coordinates, a deformation with gradient \(\textbf{F}(\textbf{X})\) of differentiability class \(\mathcal{C}^{2}\), such that \(\det \textbf{F}(\textbf{X})>0\), can be maintained for all \(W^{(nc)}\) by the application of surface tractions only (without body forces) if and only if both \(\textbf{F}\) and \(\textbf{G}\) are constant in Cartesian coordinates.
Proof
We set the notation \(\overline{\textbf{F}}=\textbf{F}\textbf{G}_{0}\), so that \(\overline{\textbf{F}}=\textbf{G}\textbf{A}\) (see Fig. 1). Since \(\textbf{G}_{0}\) is constant, this amounts to a change of coordinates from the cross-linking state, with coordinates \((X_{1},X_{2},X_{3})\), to a virtual isotropic state, with coordinates \((\overline{X}_{1},\overline{X}_{2},\overline{X}_{3})\), such that
Henceforth, all calculations are carried out within this new system of coordinates.
The first Piola-Kirchhoff stress takes the form
In the absence of body forces, the equation of elastostatic equilibrium is
Equivalently,
By the standard argument, if this holds for arbitrary \(W^{(nc)}\), then
and
In particular, if \(i=j\), then
For \(i=1\), \(\partial I^{(nc)}_{1}/\partial \overline{\textbf{F}}=2 \overline{\textbf{F}}\), hence
If the components of the deformation gradient are expressed in Cartesian coordinates as
then \(\mathrm{Div}\ \overline{\textbf{F}}=\textbf{0}\) takes the form
and \(\mathrm{Grad}\ I^{(nc)}_{1}=\textbf{0}\) can be written as
Therefore,
Noting that the right-hand side of the above equality is a sum of squares, it follows that
Hence \(\overline{\textbf{F}}\) is constant in the Cartesian coordinates \((\overline{X}_{1},\overline{X}_{2},\overline{X}_{3})\), and since \(\textbf{G}_{0}\) is constant, \(\textbf{F}\) is also constant in the Cartesian coordinates \((X_{1},X_{2},X_{3})\).
Next, we show that \(\textbf{G}\) is constant as well. For \(i=4\),
Taking any finite deformation with \(\overline{\textbf{F}}\) constant, we choose the Cartesian coordinates along its principal directions, such that \(\overline{\textbf{F}}=\mathrm{diag}\left (\overline{F}_{11}, \overline{F}_{22}, \overline{F}_{33}\right )\). Then the first equation in \(I^{(nc)}_{4}\) implies:
while the second equation implies
From the above equations, we infer that \(n_{i}\) is independent of \(\overline{X}_{i}\), for \(i=1,2,3\), and the following conditions are satisfied simultaneously:
By the isotropy condition (C2), we can apply any rotation \(\textbf{R}\) prior to the deformation, such that \(\left \{\overline{F}_{11}, \overline{F}_{22}, \overline{F}_{33} \right \}\) permute along the diagonal, while \(\textbf{n}\) does not rotate. After replacing \(\mathrm{diag}\left (\overline{F}_{11}, \overline{F}_{22}, \overline{F}_{33}\right )\) first with \(\mathrm{diag}\left (\overline{F}_{22}, \overline{F}_{33}, \overline{F}_{11}\right )\) and second with \(\mathrm{diag}\left (\overline{F}_{33}, \overline{F}_{11}, \overline{F}_{22}\right )\), we also obtain, respectively:
and
We distinguish the following two cases. First, if at least two of the three diagonal values \(\left \{\overline{F}_{11}, \overline{F}_{22}, \overline{F}_{33} \right \}\) are different from each other, then \(\mathrm{Grad}\ \textbf{n}=\textbf{0}\), i.e., \(\textbf{n}\) is uniform in the Cartesian coordinates \((\overline{X}_{1},\overline{X}_{2},\overline{X}_{3})\), and since \(\textbf{G}_{0}\) is constant, \(\textbf{n}\) is also uniform in the Cartesian coordinates \((X_{1},X_{2},X_{3})\) (see also [34]). Thus \(\textbf{G}\) is constant. Second, if \(\overline{F}_{11}=\overline{F}_{22}=\overline{F}_{33}=\overline{F}\), then \(\overline{\textbf{F}}=\overline{F}\textbf{I}\), where \(\overline{F}\) is a constant scalar in the Cartesian coordinates \((\overline{X}_{1},\overline{X}_{2},\overline{X}_{3})\). Hence, the current configuration is an isotropic state, i.e., \(\textbf{G}\) is constant, in the Cartesian coordinates \((\overline{X}_{1},\overline{X}_{2},\overline{X}_{3})\). Because \(\textbf{G}_{0}\) is constant, it follows that \(\textbf{G}\) is constant also in the Cartesian coordinates \((X_{1},X_{2},X_{3})\), i.e., either \(\textbf{n}=\textbf{n}_{0}\) or both the reference and current configurations are isotropic.
Conversely, when \(\textbf{F}\) and \(\textbf{G}\) are constant in the Cartesian coordinates, given that \(\textbf{G}_{0}\) is constant, the equation of elastostatic equilibrium in the absence of body forces is automatically satisfied. □
We note that, when equal and opposite homogeneous shear deformations generate alternating shear stripes in LCEs as in soft (or semi-soft) elastic phenomena, the deformation gradient in two adjacent stripes are geometrically compatible (see, e.g., [22, Chap. 6]). Therefore, by assuming that, the deformation gradient is piecewise of differentiability class \(\mathcal{C}^{2}\), the result of the above theorem can be extended to shear-striping patterns, which can also be maintained universally by the application of surface tractions only in ideal nematic elastomers. This is formally presented in the next result.
Corollary 1
For an ideal unconstrained uniaxial nematic LCE characterized by equation (1), such that \(\textbf{G}_{0}\) is constant in Cartesian coordinates, a deformation with gradient \(\textbf{F}(\textbf{X})\) which is piecewise of differentiability class \(\mathcal{C}^{2}\), such that \(\det \textbf{F}(\textbf{X})>0\), can be maintained for all \(W^{(nc)}\) by the application of surface tractions only (without body forces) if and only if both \(\textbf{F}\) and \(\textbf{G}\) are piecewise constant in Cartesian coordinates. In two adjacent subdomains, the deformation gradient, \(\textbf{F}_{+}\) and \(\textbf{F}_{-}\), must be rank-one connected, i.e., \(\mathrm{rank}\left (\textbf{F}_{+}-\textbf{F}_{-}\right )=1\), for geometric compatibility.
5 Conclusion
For unconstrained uniaxial nematic LCEs described by a homogeneous isotropic strain-energy density function, we proved that the only deformations that are independent of the material parameters and can be maintained by the application of traction forces on the boundary of the body, assuming that the deformation gradient is sufficiently smooth, are the homogeneous deformations, i.e., those for which the deformation gradient is constant and the director is either uniform or randomly oriented in the Cartesian coordinates. This result is consistent with the classical Ericksen’s theorem for nonlinear homogeneous isotropic hyperelastic materials, which is directly applicable to ideal LCEs in an isotropic phase. However, in a nematic phase, the constitutive strain-energy density is a function of five invariants instead of three, to account also for the liquid crystal orientation.
Since LCEs are to a large extent incompressible materials, and we know from the elastic case that this extra constraint guarantees the existence of additional universal solutions, there are also important inhomogeneous deformations that can be controlled by the application of surface tractions and are independent of the constitutive parameters. We refer to [1, 20, 24] where such deformations have been examined.
Data Availability
This paper has no data.
References
Carlson, D.E., Fried, E., Sellers, S.: Force-free states, relative strain, and soft elasticity in nematic elastomers. J. Elast. 69, 161–180 (2002). https://doi.org/10.1023/A:1027377904576
Conti, S., DeSimone, A., Dolzmann, G.: Soft elastic response of stretched sheets of nematic elastomers: a numerical study. J. Mech. Phys. Solids 50, 1431–1451 (2002). https://doi.org/10.1016/S0022-5096(01)00120-X
de Gennes P.G., Prost, J.: The Physics of Liquid Crystals, 2nd edn. Clarendon, Oxford (1993)
de Jeu, W.H. (ed.): Liquid Crystal Elastomers: Materials and Applications Springer, New York (2012)
DeSimone, A., Dolzmann, G.: Material instabilities in nematic elastomers. Phys. D, Nonlinear Phenom. 136(1–2), 175–191 (2000). https://doi.org/10.1016/S0167-2789(99)00153-0
DeSimone, A., Teresi, L.: Elastic energies for nematic elastomers. Eur. Phys. J. E 29, 191–204 (2009). https://doi.org/10.1140/epje/i2009-10467-9
Ericksen, J.L.: Deformations possible in every isotropic, incompressible, perfectly elastic body. Z. Angew. Math. Phys. 5, 466–489 (1954)
Ericksen, J.L.: Deformation possible in every compressible isotropic perfectly elastic materials. J. Math. Phys. 34, 126–128 (1955)
Finkelmann, H., Kundler, I., Terentjev, E.M., Warner, M.: Critical stripe-domain instability of nematic elastomers. J. Phys. II 7, 1059–1069 (1997). https://doi.org/10.1051/jp2:1997171
Finkelmann, H., Greve, A., Warner, M.: The elastic anisotropy of nematic elastomers. Eur. Phys. J. E 5, 281–293 (2001). https://doi.org/10.1007/s101890170060
Fried, E., Sellers, S.: Free-energy density functions for nematic elastomers. J. Mech. Phys. Solids 52(7), 1671–1689 (2004). https://doi.org/10.1016/j.jmps.2003.12.005
Fried, E., Sellers, S.: Orientational order and finite strain in nematic elastomers. J. Chem. Phys. 123(4), 043521 (2005). https://doi.org/10.1063/1.1979479
Fried, E., Sellers, S.: Soft elasticity is not necessary for striping in nematic elastomers. J. Appl. Phys. 100, 043521 (2006). https://doi.org/10.1063/1.2234824
Golubović, L., Lubensky, T.C.: Nonlinear elasticity of amorphous solids. Phys. Rev. Lett. 63, 1082–1085 (1989). https://doi.org/10.1103/PhysRevLett.63.1082
Goodbrake, C., Yavari, A., Goriely, A.: The anelastic Ericksen problem: universal deformations and universal eigenstrains in incompressible nonlinear anelasticity. J. Elast. 142(2), 291–381 (2020). https://doi.org/10.1016/j.jmps.2019.103782
Higaki, H., Takigawa, T., Urayama, K.: Nonuniform and uniform deformations of stretched nematic elastomers. Macromolecules 46, 5223–5231 (2013). https://doi.org/10.1021/ma400771z
Klingbeil, W.W., Shield, R.T.: On a class of solutions in plane finite elasticity. Z. Angew. Math. Phys. 17, 489–511 (1966). https://doi.org/10.1007/BF01595984
Kundler, I., Finkelmann, H.: Strain-induced director reorientation in nematic liquid single crystal elastomers. Macromol. Rapid Commun. 16, 679–686 (1995). https://doi.org/10.1002/marc.1995.030160908
Kundler, I., Finkelmann, H.: Director reorientation via stripe-domains in nematic elastomers: influence of cross-link density, anisotropy of the network and smectic clusters. Macromol. Chem. Phys. 199, 677–686 (1998)
Lee, V., Bhattacharya, K.: Universal deformations of incompressible nonlinear elasticity as applied to ideal liquid crystal elastomers. J. Elast. (2023). https://doi.org/10.1007/s10659-023-10018-9
Marris, A.W., Shiau, J.F.: Universal deformations in isotropic incompressible hyperelastic materials when the deformation tensor has equal proper values. Arch. Ration. Mech. Anal. 36, 135–160 (1970). https://doi.org/10.1007/BF00250814
Mihai, L.A.: Stochastic Elasticity: A Nondeterministic Approach to the Nonlinear Field Theory. Springer, Cham, Switzerland (2022). https://doi.org/10.1007/978-3-031-06692-4
Mihai, L.A., Goriely, A.: Likely striping in stochastic nematic elastomers. Math. Mech. Solids 25(10), 1851–1872 (2020). https://doi.org/10.1177/1081286520914958
Mihai, L.A., Goriely, A.: Instabilities in liquid crystal elastomers. Mater. Res. Soc. Bull. 46, 784–794 (2021). https://doi.org/10.1557/s43577-021-00115-2
Mihai, L.A., Mistry, D., Raistrick, T., Gleeson, H.F., Goriely, A.: A mathematical model for the auxetic response of liquid crystal elastomers. Philos. Trans. R. Soc. A 380, 20210326 (2022). https://doi.org/10.1098/rsta.2021.0326
Petelin, A., Čopič, M.: Observation of a soft mode of elastic instability in liquid crystal elastomers. Phys. Rev. Lett. 103, 077801 (2009). https://doi.org/10.1103/PhysRevLett.103.077801
Petelin, A., Čopič, M.: Strain dependence of the nematic fluctuation relaxation in liquid-crystal elastomerss. Phys. Rev. E 82, 011703 (2010). https://doi.org/10.1103/PhysRevE.82.011703
Saccomandi, G.: Universal solutions and relations in finite elasticity. In: Hayes, M., Saccomandi, G. (eds.) Topics in Finite Elasticity, pp. 95–130. Springer, Wien (2001)
Shield, R.T.: Deformations possible in every compressible, isotropic, perfectly elastic material. J. Elast. 1, 91–92 (1971)
Spencer, A.J.M.: Theory of invariants. In: Eringen, A.C. (ed.) Continuum Physics 1, pp. 239–253. Academic Press, New York (1971)
Talroze, R.V., Zubarev, E.R., Kuptsov, S.A., Merekalov, A.S., Yuranova, T.I., Plate, N.A., Finkelmann, H.: Liquid crystal acrylate-based networks: polymer backbone-LC order interaction. React. Funct. Polym. 41, 1–11 (1999). https://doi.org/10.1016/S1381-5148(99)00032-2
Warner, M., Terentjev, E.M.: Liquid Crystal Elastomers. Oxford University Press, Oxford (2007)
Yavari, A., Goriely, A.: The anelastic Ericksen problem: universal eigenstrains and deformations in compressible isotropic elastic solids. Proc. R. Soc. A 472(2196), 20160690 (2016). https://doi.org/10.1098/rspa.2016.0690
Yavari, A., Goriely, A.: Universal deformations in anisotropic nonlinear elastic solids. J. Mech. Phys. Solids 156, 104598 (2021). https://doi.org/10.1016/j.jmps.2021.104598
Yavari, A., Goriely, A.: The universal program of nonlinear hyperelasticity. J. Elast., 1–56 (2022). https://doi.org/10.1007/s10659-022-09906-3
Zubarev, E.R., Kuptsov, S.A., Yuranova, T.I., Talroze, R.V., Finkelmann, H.: Monodomain liquid crystalline networks: reorientation mechanism from uniform to stripe domains. Liq. Cryst. 26, 1531–1540 (1999). https://doi.org/10.1080/026782999203869
Acknowledgement
The authors would like to thank the Isaac Newton Institute for Mathematical Sciences for support and hospitality during the programme Uncertainty Quantification and Stochastic Modelling of Materials when work on this paper was undertaken.
Funding
This work was supported by EPSRC Grant Number EP/R014604/1.
Author information
Authors and Affiliations
Contributions
Both the authors contributed equally to the design, analysis and presentation of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Mihai, L.A., Goriely, A. Controllable Deformations of Unconstrained Ideal Nematic Elastomers. J Elast 156, 95–106 (2024). https://doi.org/10.1007/s10659-023-10038-5
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10659-023-10038-5