Original Research ARTICLE
Modeling Defects, Shape Evolution, and Programmed Auto-Origami in Liquid Crystal Elastomers
- 1Liquid Crystal Institute, Kent State University, Kent, OH, USA
- 2Department of Chemical Engineering, Columbia University, New York, NY, USA
Liquid crystal elastomers represent a novel class of programmable shape-transforming materials whose shape-change trajectory is encoded in the material’s nematic director field. Using three-dimensional non-linear finite element elastodynamics simulation, we model a variety of different actuation geometries and device designs: thin films containing topological defects, patterns that induce formation of folds and twists, and a bas-relief structure. The inclusion of finite bending energy in the simulation model reveals features of actuation trajectory that may be absent when bending energy is neglected. We examine geometries with a director pattern uniform through the film thickness encoding multiple regions of positive Gaussian curvature. Simulations indicate that heating such a system uniformly produces a disordered state with curved regions emerging randomly in both directions due to the film’s up/down symmetry. In contrast, applying a thermal gradient by heating the material first on one side breaks up/down symmetry and results in a deterministic trajectory, producing a more ordered final shape. We demonstrate that a folding zone design containing cut-out areas accommodates transverse displacements without warping or buckling, and demonstrate that bas-relief and more complex bent/twisted structures can be assembled by combining simple design motifs.
Liquid crystal elastomers (LCEs) undergo reversible shape transformations under any stimulus that changes their degree of nematic order, including heat, illumination, or change of chemical environment (White and Broer, 2015). LCEs are cross-linked polymers with liquid crystal mesogens either incorporated along the polymer main chain or attached as side-chains. The degree of orientational order of the mesogens within an LCE depends on temperature and may display a nematic–isotropic transition on heating. Due to coupling between the mesogens and the polymer, such a phase transition may give rise to macroscopic mechanical strains of several hundred percent.
Complex shape transformation trajectories may be encoded in the material by patterning of the nematic director when the polymer is cross-linked (Liu et al., 2015). The director field is typically controlled by forming samples between flat substrates with a surface anchoring pattern (de Haan et al., 2014). The anchoring pattern can be identical on both substrates, such that for a thin enough sample, the director field is uniform through the material’s thickness. Alternatively, the anchoring pattern may differ between the substrates by a prescribed twist angle, where the sense of director twist between the substrates is controlled by the addition of chiral dopant (Sawa et al., 2011). If the material is cross-linked in the nematic phase and the cross-link density is sufficiently high, the director field is “blueprinted” in the material; this programmed pattern then gives rise to locally anisotropic deformation when the scalar nematic order parameter changes, e.g., on heating. As the resulting programmable shape transformations are reversible, these materials are attractive candidates for applications as conformable soft actuators (Spillmann et al., 2007), to control surface textures (McConney et al., 2013), or to open/close valves or apertures (Sánchez-Ferrer et al., 2011; Modes et al., 2013; Schuhladen et al., 2014).
In previous work, we collaborated with experimenters to explore how even relatively simple director patterns give rise to complex actuation behavior. For instance, a flat, straight LCE ribbon with director twist through the thickness can form either helicoid or spiral structures, depending on sample properties including director twist geometry and sample aspect ratio; remarkably, these chiral structures undergo transitions from right- to left-handed with temperature (Sawa et al., 2011, 2013). A larger sample with multiple twisted domains arrayed in stripes gives rise to actuation via accordion folds, while a checkerboard pattern of twisted domains produces a buckled texture (de Haan et al., 2013). In each case, we used finite element method simulations to investigate the relationship between director structure and resulting shape actuation. We also examined the soft elastic response of LCE materials with lower cross-link density; these materials show soft elastic response characterized by formation of striped director texture under strain, producing an extended plateau in the stress-strain curve (Mbanga et al., 2010).
Here, we explore a variety of programmed shape transformations in LCEs with more complex blueprinted director fields, inspired in many cases by recent experiments. An explanation of simulation methods is provided below in the final section.
Accurate simulation tools are needed for design/engineering of future devices, to explore dependence of shape transformation geometries on blueprinted microstructure, and to test how known design motifs combine. Simulations can in particular provide insight into observed differences between analytical predictions – which often neglect bending energy – and experimental observations. Here, we examine several geometries that have been studied experimentally and go on to explore several novel designs.
Modes et al. (2010, 2011) predicted that an initially flat LCE film containing a +1 disclination in the director field, with the director oriented in the azimuthal direction, deforms into the shape of a cone on heating and an anticone (saddle) on cooling. Recent experimental work by McConney et al. (2013) and Ware et al. (2015) demonstrated this behavior and also examined higher order topological defects. Analytical predictions were made in the approximation that bending energy is neglected. We use finite element simulation to examine expected deformation trajectories in the case where bending energy is included.
Figure 1 shows a finite element simulation of an initially square LCE film with an azimuthal +1 defect. The film is represented by a 3D tetrahedral mesh with an aspect ratio of 50:1. As predicted, we find formation of a cone on heating, and formation of a saddle on cooling, with respect to the cross-link temperature. We note that on heating, the sample initially remains flat in a metastable state until at a high enough temperature the cone begins to emerge outward, spontaneously breaking up–down symmetry. This behavior is driven by the inclusion of a finite bending energy in the simulation model.
Figure 1. (A) Finite element simulations – LCE with an azimuthal +1 defect, and aspect ratio 50:1, deforms into a cone on heating and an anticone (saddle) on cooling. (B) Cone height as a function of ΔS, the reduction of nematic scalar order parameter on heating. The sample initially remains flat on heating, in a metastable state. The cone emerges only when heated sufficiently that ΔS > 0.2.
To further investigate effects of finite bending energy, we carry out another simulation of a film containing a +1 defect using a disk-shaped sample; its shape evolution is shown in Figure 2. The transient structure shows formation of a “sombrero” shape with a distinct rim region, which gradually flattens and snaps through to reach the expected final cone shape. The observed transient shapes may provide insight into structures observed in experiment (de Haan et al., 2012; McConney et al., 2013).
Figure 2. Finite element simulation of an LCE disk with a +1 azimuthal defect, on heating. We observe a transient “sombrero” shape with an upward-curving rim. The rim deforms and snaps through, leaving a cone shaped center surrounded by nearly flat annulus. Eventually, the sample reaches the theoretically predicted cone shape (Modes et al., 2010).
We also investigate shape transformations in LCE films containing higher order defects. Shown in Figure 3 is a film with a −2 defect. The resulting actuated shape has a protruding center region surrounded by six dimples. Because the top and bottom of the film are symmetric, each of these regions of positive Gaussian curvature has equal probability to pop upwards or downwards. To enhance the probability of finding a symmetric final state, we biased the transformation by applying a temperature gradient at the start of the simulation, heating the sample from one side before stabilizing at uniform temperature.
Figure 3. Left: structure of a −2 defect. Right: finite element simulation of an LCE disk containing a −2 defect, shape transformation on heating. Up/down symmetry is broken by heating from one side, producing a symmetric structure. Color indicates local elastic strain energy.
Figure 4 (center) shows the final state of an initially flat LCE film containing a −4 defect, again heated from one side to achieve a symmetric structure. Figure 4 (right) shows the structure that forms when the sample is heated uniformly and different regions pop upwards or downwards at random, producing a disordered and asymmetric shape. These images can be compared to recent experiments by McConney et al. (2013) on high-order defects.
Figure 4. Left: structure of a −4 defect. Center: finite element simulation of an LCE disk with a −4 defect, shape transformation on heating, heated initially from above. Right: if heated uniformly, each region of positive Gaussian curvature may emerge either upward or downward, producing a complex, disordered final shape.
In order to achieve deterministic rather than random shape evolution, it is necessary to break this up–down symmetry. Instead of applying a thermal gradient, a better option would be to introduce a material gradient. For instance, the cross-link density may be varied through the thickness of the film, by introducing a UV absorbing component in the mixture and cross-link via illumination from one side. Alternatively, the director field on top and bottom substrates can be altered so that they are not identical.
Inducing a Fold or Twist
“Auto-origami” refers to spontaneous deformation of a programmable material whose microstructure encodes a sequence of folds, stretches, and shears to produce a trajectory of motion from one stable conformation to another, triggered by an external stimulus such as a change of temperature. A key goal in design of shape-changing auto-origami devices is the formation of folds, actuated, e.g., on heating, induced by a localized gradient in the nematic director field between the top and bottom surfaces. For instance, director twist through the thickness of the material can be used to induce a fold (Ware et al., 2015). One challenge is that in addition to driving a folding deformation, an imposed director gradient of this type may also induce a more complex deformation that gives rise to unintended twist or buckling of the sample.
To localize the stresses associated with bending deformation to the hinge region, we designed a structure containing holes. In Figure 5, we show a bending joint composed of two straight (non-LCE) elastic side walls connected by an array of short LCE beams, each of which has the director field with homeotropic anchoring on top, planar anchoring on bottom parallel to the beam’s long axis, with a smooth director gradient in between. On heating, when each beam bends, the accompanying transverse deformation is accommodated by the adjacent voids. As a result, the overall structure bends smoothly without any transverse shrinkage or buckling. This design could function well as an adhesive tape which, when adhered to paper or fabric, induces a fold at a specified temperature. Bending tapes that respond at different temperatures could induce a series of folds that occur in sequence, e.g., on gradual warming.
Figure 5. Finite element simulation of an actuator that bends on heating; color indicates strain energy density. Non-LCE side walls are connected by short LCE beams, each of which has homeotropic director orientation along the top surface and planar orientation at the bottom surface with a smooth gradient in between. Open spaces along the bending zone accommodate transverse strain of the LCE beams without lateral shrinkage or buckling of the straight side walls.
We also demonstrate combining two different types of director gradients along a beam in order to construct a more complex shape evolution. Figure 6 shows a beam mounted at one end on a substrate, with two different director domains. The left 30% of the sample has a splay geometry with homeotropic anchoring on the bottom and planar anchoring on the top; on heating, it curves to form a C-shaped bend. The remaining part of the sample has planar anchoring on both sides with a 90° twist between top and bottom surfaces, like those we previously modeled (Sawa et al., 2011, 2013), and on heating twists to form a helix. Time evolution on heating produces the trajectory shown in Figure 6. This structure can be reconfigured reversibly by changing temperature.
Figure 6. Finite element simulation of shape evolution of an LCE beam on heating. Left 30% of beam has splayed nematic director, which induces a C-shaped bend, while the remaining beam has a 90° director twist through the thickness, giving rise to a helical twist.
Actuators in the form of bas-relief or raised patterns are desirable for applications including haptic displays and for dynamic control of surface texture. A variety of director motifs have been proposed to pattern surfaces to achieve desired surface structures including director fields with (McConney et al., 2013) and without (Mostajeran et al., 2015) topological defects.
Here, we examine a simple design motif to drive emergence of a bas-relief pattern on heating. To create a raised region in the shape of a valentine heart, we model a square LCE film with different anchoring conditions on the top and bottom surfaces. On the top surface, we impose homeotropic anchoring inside a heart-shaped region, and uniform planar anchoring outside, and on the bottom surface, we reverse the two anchoring conditions, as shown in Figure 7. In between, we assume a smooth variation of the nematic director without defects. In practice, it may be necessary to impose a small pre-tilt angle in the homeotropic anchored regions in order to prevent defect formation. We model shape transformation of this system on heating and find formation of a raised heart-shaped bas-relief structure. Unlike the patterns driven by topological defects of the type discussed above, this design lacks up–down symmetry and thus actuates deterministically with no randomness in the final state.
Figure 7. (A) Anchoring patterns on top and bottom surfaces with a combination of planar and homeotropic domains. (B) Finite element simulation of actuation on heating of the resulting LCE film; the heart-shaped domain emerges.
To create a three-dimensional non-linear finite element model of our system, we start by meshing our sample in an unstructured tetrahedral mesh using the open-source application Salome (available at http://www.salome-platform.org/). We write a Hamiltonian for the system with both potential and kinetic energy terms as
Here, the first term in the sum over elements is the elastic strain energy, where εij is the Green-Lagrange strain, assumed piecewise constant over each volume element, and Cijkl is a tensor of elastic constants. The second term in the sum over elements represents coupling between the strain and the nematic order tensor Qij, where represents the “blueprinted” nematic order, defined in each element, at the moment when the sample was cross-linked. The parameter α represents the strength of coupling between nematic order and strain and governs how much the material deforms when the state of nematic order changes, while Vn represents the volume of the element measured in the reference state. The sum over nodes represents calculation of the kinetic energy, where mp is the effective mass of each vertex, calculated using the lumped mass approximation, and vp is the node’s velocity.
When the sample is heated, we assume that the blueprinted nematic director, described as piecewise constant within each element, remains fixed in the body reference frame but the scalar order parameter drops. Thus the nematic order tensor within each element takes the form , where S(T) is the scalar order parameter as a function of temperature and ni represents the components of the nematic director. Any change in the degree of scalar order thus changes Qij and drives deformation.
The system evolves forward in time via explicit time stepping in a manner reminiscent of molecular dynamics. At the start of the simulation, all velocities are initially set to 0. At each time step, the effective force on each node in the mesh is calculated as described below, and the equations of motion for the nodes are integrated using the velocity Verlet method (Press et al., 2007).
The force calculation proceeds as follows. The displacement field is linear interpolated within each tetrahedral element using the positions of the corner nodes in the current state and in the element’s undeformed, stress-free reference state, which we take to be the shape of the element when the material is cross-linked. From the gradients of the displacement field, we calculate the components of the Green-Lagrange strain tensor and from that, calculate the potential energy. In this manner, we write the potential energy of the system as a function of the node positions. We calculate the effective force on each of the four nodes by taking the negative derivative of the potential energy within each element with respect to the node positions. The effective force on each node receives a contribution from each element of which it is a member.
Interestingly, our finite element algorithm closely resembles a conventional molecular dynamics simulation (Rapaport, 2004), where nodes in the finite element mesh move instead of atoms. Node interactions take the form of a many body potential, calculated for each element, with energy and forces depending on the node positions both in their current state and in their initial, stress-free reference state. Like a molecular dynamics algorithm, this calculation is entirely local and can thus be easily parallelized for fast execution on a cluster or GPU-equipped processor. If no dissipative forces or mechanical constraints are added, the calculated dynamics conserves the sum of kinetic and potential energy, momentum, and angular momentum. In practice, we add dissipative forces to allow the sample to reach mechanical equilibrium; motion may thus be under- or over-damped. The capability to model momentum-conserving elastodynamics, and not just quasistatic relaxation, allows us to model snap-through shape transitions (Shankar et al., 2013; Smith et al., 2014) and other complex dynamical behavior.
Within this model, if two parts of the sample intersect, they pass through each other, ghostlike, without any collision. To avoid this unphysical behavior in any simulation, we add short-range pairwise repulsive interactions between the surface nodes using the Weeks–Chandler–Andersen potential, which is the Lennard-Jones potential truncated at its minimum and shifted such that the potential vanishes at the cutoff (Ahmed and Sadus, 2009). We then use a molecular dynamics-type algorithm to detect collisions between such surface nodes and calculate the resulting forces.
While less computationally efficient than a two-dimensional shell model (Chung et al., 2015), our fully three-dimensional model can describe material systems with director or cross-link density gradients across the sample thickness. With GPU implementation we can achieve speeds of up to 100 time steps per CPU second for a mesh of more than 105 elements. More details of the model may be found in Mbanga (2012) and Gimenez-Pinto (2014).
Devices engineered using patterned LCEs have a somewhat limited design space as the nematic director field may be controlled only via anchoring patterns at the sample’s surface. For this reason, an arbitrary three-dimensional director field cannot be realized. An experimental method to blueprint LCE materials with fully three-dimensional voxel-by-voxel control of the director orientation would provide a broader class of available structures, but has not yet been developed to the best of our knowledge.
In contrast, other materials with anisotropic shape evolution have recently been introduced (Gladman et al., 2016) that can be printed in successive layers, providing a greater degree of control over the full three-dimensional orientation field. Finite element modeling will play a key role in engineering both types of devices. Future developments will include numerical methods to address the so-called inverse problem, that is, to design a director field that gives rise to a shape actuation that closely follows a prescribed trajectory and achieves a target final shape.
Finite element simulation is a valuable design tool to explore a wide variety of device designs for materials undergoing programmable shape change. The three-dimensional nature of our model allows study of systems where bending energy plays an important role, and allows modeling of samples where the nematic director is not uniform through the film thickness.
While quasistatic relaxation is sufficient to model a slow shape evolution process, some shape transitions involve mechanical instabilities and may show rapid snap-through behavior. These rapid shape transitions can be modeled via underdamped, momentum-conserving elastodynamics simulation. We plan to further explore such systems in future work.
RS in the primary investigator of this work and developed the simulation method used in this paper. VG-P implemented the simulation for several different director structures, and performed data analysis and visualization of the results. AK implemented the simulation algorithm in the CUDA programming language for execution on a GPU-enabled computer, producing dramatically faster computational speed, and carried out design and shape-change simulations for several of the designs listed here. He also developed improved data visualization methods.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors thank Tim White, Mike McConney, Qihuo Wei, and Badel Mbanga for helpful discussions. Work supported by the National Science Foundation (NSF-DMR 1106014, NSF-DMR-1409658, and NSF-CMMI 1436565). Computer resources provided by the Ohio Supercomputer Center and by the Kent State College of Arts and Sciences.
Ahmed, A., and Sadus, R. J. (2009). Phase diagram of the Weeks-Chandler-Andersen potential from very low to high temperatures and pressures. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 80, 061101. doi: 10.1103/PhysRevE.80.061101
Chung, H., Choi, J., Yun, J.-H., and Cho, M. (2015). Light and thermal responses of liquid-crystal-network films: a finite element study. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 91, 042503. doi:10.1103/PhysRevE.91.042503
de Haan, L. T., Gimenez-Pinto, V., Konya, A., Nguyen, T.-S., Verjans, J. M. N., Sánchez-Somolinos, C., et al. (2013). Accordion-like actuators of multiple 3D patterned liquid crystal polymer films. Adv. Funct. Mater. 24, 1251–1258. doi:10.1002/adfm.201302568
de Haan, L. T., Sánchez-Somolinos, C., Bastiaansen, C. M. W., Schenning, A. P. H. J., and Broer, D. J. (2012). Engineering of complex order and the macroscopic deformation of liquid crystal polymer networks. Angew. Chem. Int. Ed. Engl. 51, 12469–12472. doi:10.1002/anie.201205964
Gimenez-Pinto, V. K. (2014). Modeling Liquid Crystal Polymeric Devices. Available at: http://rave.ohiolink.edu/etdc/view?acc_num=kent1398039813
McConney, M. E., Martinez, A., Tondiglia, V. P., Lee, K. M., Langley, D., Smalyukh, I. I., et al. (2013). Topography from topology: photoinduced surface features generated in liquid crystal polymer networks. Adv. Mater. 25, 5880–5885. doi:10.1002/adma.201301891
Modes, C. D., Warner, M., Sanchez-Somolinos, C., de Haan, L. T., and Broer, D. (2013). Angular deficits in flat space: remotely controllable apertures in nematic solid sheets. Proc. R. Soc. A Math. Phys. Eng. Sci. 469, 631. doi:10.1098/rspa.2012.0631
Mostajeran, C., Ware, T. H., and White, T. J. (2015). Encoding Gaussian Curvature in Glassy and Elastomeric Liquid Crystal Polymer Networks. 6. Available at: http://arxiv.org/abs/1510.00803
Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P. (2007). Numerical Recipes 3rd Edition: The Art of Scientific Computing. Available at: http://dl.acm.org/citation.cfm?id=1403886
Rapaport, D. C. (2004). The Art of Molecular Dynamics Simulation. Available at: http://dl.acm.org/citation.cfm?id=2222483
Sánchez-Ferrer, A., Fischl, T., Stubenrauch, M., Albrecht, A., Wurmus, H., Hoffmann, M., et al. (2011). Liquid-crystalline elastomer microvalve for microfluidics. Adv. Mater. 23, 4526–4530. doi:10.1002/adma.201102277
Sawa, Y., Urayama, K., Takigawa, T., Gimenez-Pinto, V., Mbanga, B. L., Ye, F., et al. (2013). Shape and chirality transitions in off-axis twist nematic elastomer ribbons. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 88, 022502. doi:10.1103/PhysRevE.88.022502
Sawa, Y., Ye, F., Urayama, K., Takigawa, T., Gimenez-Pinto, V., Selinger, R. L. B., et al. (2011). Shape selection of twist-nematic-elastomer ribbons. Proc. Natl. Acad. Sci. U.S.A. 108, 6364–6368. doi:10.1073/pnas.1017658108
Schuhladen, S., Preller, F., Rix, R., Petsch, S., Zentel, R., and Zappe, H. (2014). Iris-like tunable aperture employing liquid-crystal elastomers. Adv. Mater. 26, 7247–7251. doi:10.1002/adma.201402878
Shankar, M. R., Smith, M. L., Tondiglia, V. P., Lee, K. M., McConney, M. E., Wang, D. H., et al. (2013). Contactless, photoinitiated snap-through in azobenzene-functionalized polymers. Proc. Natl. Acad. Sci. U.S.A. 110, 18792–18797. doi:10.1073/pnas.1313195110
Smith, M. L., Shankar, M. R., Backman, R., Tondiglia, V. P., Lee, K. M., McConney, M. E., et al. (2014). “Designing light responsive bistable arches for rapid, remotely triggered actuation,” in SPIE Smart Structures and Materials + Nondestructive Evaluation and Health Monitoring, eds N. C. Goulbourne and H. E. Naguib (San Diego, CA: International Society for Optics and Photonics), 90580F.
Spillmann, C. M., Naciri, J., Martin, B. D., Farahat, W., Herr, H., and Ratna, B. R. (2007). Stacking nematic elastomers for artificial muscle applications. Sensors Actuators A Phys. 133, 500–505. doi:10.1016/j.sna.2006.04.045
Keywords: liquid crystals, polymers, shape-programmable materials, finite element simulation, non-linear elasticity
Citation: Konya A, Gimenez-Pinto V and Selinger RLB (2016) Modeling Defects, Shape Evolution, and Programmed Auto-Origami in Liquid Crystal Elastomers. Front. Mater. 3:24. doi: 10.3389/fmats.2016.00024
Received: 09 February 2016; Accepted: 02 May 2016;
Published: 01 June 2016
Edited by:Oliver Bäumchen, Max-Planck-Institute for Dynamics and Self-Organization, Germany
Reviewed by:Dong Chen, Zhejiang University, China
Shashi Thutupalli, Princeton University, USA
Aurora Hernandez-Machado, Universitat de Barcelona, Spain
Copyright: © 2016 Konya, Gimenez-Pinto and Selinger. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Robin L. B. Selinger, email@example.com