Skip to main content

Numerical study of liquid crystal elastomers by a mixed finite element method

  • C. LUO (a1) and M. C. CALDERER (a2)

Liquid crystal elastomers present features not found in ordinary elastic materials, such as semi-soft elasticity and the related stripe domain phenomenon. In this paper, the two-dimensional Bladon–Terentjev–Warner model and the one-constant Oseen–Frank energy expression are combined to study the liquid crystal elastomer. We also impose two material constraints, the incompressibility of the elastomer and the unit director norm of the liquid crystal. We prove existence of minimiser of the energy for the proposed model. Next we formulate the discrete model, and also prove that it possesses a minimiser of the energy. The inf-sup values of the discrete linearised system are then related to the smallest singular values of certain matrices. Next the existence and uniqueness of the Lagrange multipliers associated with the two material constraints are proved under the assumption that the inf-sup conditions hold. Finally numerical simulations of the clamped-pulling experiment are presented for elastomer samples with aspect ratio 1 or 3. The semi-soft elasticity is successfully recovered in both cases. The stripe domain phenomenon, however, is not observed, which might be due to the relative coarse mesh employed in the numerical experiment. Possible improvements are discussed that might lead to the recovery of the stripe domain phenomenon.

Hide All
[1]Adams R. & Fournier J. (1975) Sobolev Spaces, Vol. 65, Academic Press, New York.
[2]Ball J. (1976) Convexity conditions and existence theorems in nonlinear elasticity. Arch. Ration. Mech. Anal. 63 (4), 337403.
[3]Bladon P., Terentjev E. & Warner M. (1993) Transitions and instabilities in liquid crystal elastomers. Phys. Rev. E 47 (6), 38383840.
[4]Bochkanov S. & Bystritsky V. ALGLIB [online]. URL:
[5]Brezzi F. & Fortin M. (1991) Mixed and Hybrid Finite Element Methods, Springer, Berlin.
[6]Calderer M., Liu C. & Yan B. (2006) A mathematical theory for nematic elastomers with non-uniform prolate spheroids. Advances in applied and computational mathematics, 245–259.
[7]Cesana P. & DeSimone A. (2009) Strain-order coupling in nematic elastomers: equilibrium configurations. Math. Models Methods Appl. Sci. 19, 601630.
[8]Chapelle D. & Bathe K. (1993) The inf-sup test. Comput. Struct. 47, 537537.
[9]Ciarlet P. (1988) Mathematical Elasticity, Vol 1, North-Holland, The Netherlands.
[10]Clarke F. (1975) Generalized gradients and applications. Trans. Am. Math. Soc. 205, 247262.
[11]Clarke F. (1976) A new approach to Lagrange multipliers. Math. Oper. Res. 1 (2), 165174.
[12]Clarke S., Hotta A., Tajbakhsh A. & Terentjev E. (2001) Effect of crosslinker geometry on equilibrium thermal and mechanical properties of nematic elastomers. Phys. Rev. E 64 (6), 61702.
[13]Conti S., DeSimone A. & Dolzmann G. (2002) Semisoft elasticity and director reorientation in stretched sheets of nematic elastomers. Phys. Rev. E 66 (6), 061710.
[14]Conti S., DeSimone A. & Dolzmann G. (2002) Soft elastic response of stretched sheets of nematic elastomers: aA numerical study. J. the Mech. Phys. Solids 50 (7), 14311451.
[15]de Gennes P. & Prost J. (1995) The Physics of Liquid Crystals, Oxford University Press, USA.
[16]DeSimone A. (1999) Energetics of fine domain structures. Ferroelectrics 222 (1–4), 533542.
[17]DeSimone A. & Dolzmann G. (2000) Material instabilities in nematic elastomers. Phys. D: Nonlinear Phenom. 136 (1–2), 175191.
[18]DeSimone A. & Teresi L. (2009) Elastic energies for nematic elastomers. Eur. Phys. J. E: Soft Matter and Biol. Phys. 29 (2), 191204.
[19]Ericksen J. (1991) Liquid crystals with variable degree of orientation. Arch. Ration. Mech. Anal. 113 (2), 97120.
[20]Evans L. (1998) Partial Differential Equations. American Mathematical Society, USA.
[21]Frank F. (1958) I. Liquid crystals. On the theory of liquid crystals. Discuss. Faraday Soc. 25, 1928.
[22]Hardt R., Kinderlehrer D. & Lin F. (1986) Existence and partial regularity of static liquid crystal configurations. Commun. Math. Phys. 105 (4), 547570.
[23]Hébert M., Kant R. & De Gennes P. (1997) Dynamics and thermodynamics of artificial muscles based on nematic gels. J. Physique 7 (7), 909919.
[24]Hu Q., Tai X. & Winther R. (2009) A saddle point approach to the computation of harmonic maps. SIAM J. Numer. Anal. 47 (2), 15001523.
[25]Kundler I. & Finkelmann H. (1995) Strain-induced director reorientation in nematic liquid single crystal elastomers. Macromol. Rapid Commun. 16 (9), 679686.
[26]Küpfer J. & Finkelmann H. (1994) Liquid crystal elastomers: Influence of the orientational distribution of the crosslinks on the phase behaviour and reorientation processes. Macromol. Chem. Phys. 195 (4), 13531367.
[27]Le Tallec P. (1981) Compatibility condition and existence results in discrete finite incompressible elasticity. Comput. Methods Appl. Mech. Eng. 27 (2), 239259.
[28]Logg A & Wells G. (2010) DOLFIN: Automated finite element computing. ACM Trans. Math. Softw. 37 (2), 128. URL:
[29]Luo C. (2010) Modeling, Analysis and Numerical Simulation of Liquid Crystal Elastomer, PhD Thesis, University of Minnesota, USA.
[30]Majumdar A. & Zarnescu A. (2010) Landau–De Gennes theory of nematic liquid crystals: The Oseen–Frank limit and beyond. Arch. Ration. Mech. Anal. 196 (1), 227280.
[31]Rudin W. (1991) Functional analysis. In: International Series in Pure and Applied Mathematics, McGraw-Hill, USA.
[32]Tallec P. L. (1994) Numerical methods for nonlinear three-dimensional elasticity. In: Handbook of Numerical Analysis, Vol. 3, pp. 465622. North-Holland, The Netherlands.
[33]Verwey G., Warner M. & Terentjev E. (1996) Elastic instability and stripe domains in liquid crystalline elastomers. J. Phys. II France 6, 12731290.
[34]Warner M. & Terentjev E. (2007) Liquid Crystal Elastomers, Oxford University Press, USA.
[35]Zubarev E., Kuptsov S., Yuranova T., Talroze R. & Finkelmann H. (1999) Monodomain liquid crystalline networks: Reorientation mechanism from uniform to stripe domains. Liq. Cryst. 26 (10), 15311540.
Recommend this journal

Email your librarian or administrator to recommend adding this journal to your organisation's collection.

European Journal of Applied Mathematics
  • ISSN: 0956-7925
  • EISSN: 1469-4425
  • URL: /core/journals/european-journal-of-applied-mathematics
Please enter your name
Please enter a valid email address
Who would you like to send this to? *



Full text views

Total number of HTML views: 0
Total number of PDF views: 86 *
Loading metrics...

Abstract views

Total abstract views: 164 *
Loading metrics...

* Views captured on Cambridge Core between September 2016 - 18th January 2018. This data will be updated every 24 hours.