scispace - formally typeset
Search or ask a question
Journal Article•DOI•

A computational framework for incompressible electromechanics based on convex multi-variable strain energies for geometrically exact shell theory

Rogelio Ortigosa1, Antonio J. Gil1•
Swansea University1
15 Apr 2017-Computer Methods in Applied Mechanics and Engineering (North-Holland)-Vol. 317, pp 792-816
TL;DR: In this paper, a new computational framework for the analysis of incompressible Electro Active Polymer (EAP) shells subjected to large strains and large electric fields is presented, based on a rotationless description of the kinematics of the shell, enhanced with extra degrees of freedom corresponding to the thickness stretch and the hydrostatic pressure.
Abstract: In this paper, a new computational framework for the analysis of incompressible Electro Active Polymer (EAP) shells subjected to large strains and large electric fields is presented. Two novelties are incorporated in this work. First, the variational and constitutive frameworks developed by the authors in recent publications (Gil and Ortigosa, 2016; Ortigosa and Gil, 2016; Ortigosa et al., 2016) in the context of three-dimensional electromechanics are particularised/degenerated to the case of geometrically exact shell theory. This formulation is computationally very convenient as EAPs are typically used as thin shell-like components in a vast range of applications. The proposed formulation follows a rotationless description of the kinematics of the shell, enhanced with extra degrees of freedom corresponding to the thickness stretch and the hydrostatic pressure, critical for the consideration of incompressibility. Different approaches are investigated for the interpolation of these extra fields and that of the electric potential across the thickness of the shell. Crucially, this allows for the simulation of multilayer and composite materials, which can display a discontinuous strain distribution across their thickness. As a second novelty, a continuum degenerate approach allows for the consideration of complex three-dimensional electromechanical constitutive models, as opposed to those defined in terms of the main strain measures of the shell. More specifically, convex multi-variable (three-dimensional) constitutive models, complying with the ellipticity condition and hence, satisfying material stability for the entire range of deformations and electric fields, are used for the first time in the context of shell theory.

Summary (3 min read)

Jump to: [1. Introduction] – [3.1. Shells kinematics] – [3.4. Tangent operators in incompressible electro-elasticity. Continuum degenerate shell formulation] – [3.4.1. Tangent operator of the internal energy e] – [4. Variational formulation of nearly and incompressible dielectric elastomer shells] – [5.2. Interpolation across the thickness of the shell] – [6. Numerical examples] – [6.1. Bending actuators] – [6.2. Helicoidal actuator] – [6.3. Hyperboloid piezoelectric polymer] and [7. Concluding remarks]

1. Introduction

  • Electro Active Polymers (EAPs) belong to a special class of smart materials with very attractive actuator and energy harvesting capabilities [5].
  • This rotationless approach, which complies with the principle of material frame indifference [21], avoids a well-known drawback associated with rotation-based formulations.
  • More specifically, convex multi-variable electromechanical constitutive models, satisfying the ellipticity condition and hence, material stability [27–29] for the entire range of deformations and electric fields, are used for the first time in the context of shell theory.
  • Section 3 presents the kinematical description of the proposed shell formulation.
  • Additionally, the 3 concept of multi-variable convexity is extended to the context of nonlinear shell theory.

3.1. Shells kinematics

  • This is the case for the majority of applications of EAPs, where they feature as thin shell-like components.
  • Notice that the spatial vector v does not have to be necessarily perpendicular to the plane Γ. Moreover, γ(ηα, s) in above equation (11) represents the thickness stretch [17], which accounts for possible deformations across the thickness of the shell.
  • Notice that this extra field γ might depend not only on the convective coordinates ηα, but also 7 upon s.

3.4. Tangent operators in incompressible electro-elasticity. Continuum degenerate shell formulation

  • As shown in Section 3.2, the kinematics of the shell leads to further geometrical non-linearities with respect to the continuum formulation.
  • As a result, these extra non-linearities will also be reflected in the tangent operators of the internal and Helmholtz’s energy functionals, e (20) and Φ (27), respectively4.

3.4.1. Tangent operator of the internal energy e

  • The tangent operator of both the isochoric and volumetric components of the internal energy, ê and U , respectively, can be defined for the continuum 4Refer to [4] for a comparison with the tangent operator of both the internal and Helmholtz’s energy functionals emerging in the continuum formulation.
  • Notice that these tensors introduce an additional geometrical nonlinearity, represented by the second terms on the right hand side of both tangent operators in equation (28), with respect to the tangent operators emerging in the continuum formulation, presented in Reference [4].

4. Variational formulation of nearly and incompressible dielectric elastomer shells

  • The objective of this Section is to present the variational framework for the proposed shell formulation.
  • An iterative6 Newton-Raphson process is usually preferred to converge 5The expression of the external virtual work DW ext[δu0, δv, δγ] is well known and, hence, omitted.
  • 6The letter k will indicate iteration number.
  • These recursive relationships (carried out at every Gauss point of the domain) between the Helmholtz’s energy Φ̂ and the internal energy.

5.2. Interpolation across the thickness of the shell

  • The interpolation of the uniparametric functions J a(s) is carried out via element-wise (e) continuous (or discontinuous) Lagrange polynomial interpolants of degree pJ .
  • When considering continuous (or discontinuous) interpolants, this will be denoted as Continuum-Based-Continuous (CBC) (or Continuum-Based-Discontinuous (CBD)) approach, both described as J a(s) = ns∑ e=1 pJ+1∑ b=1 J abe N bJ e(s), (47) where J abe represents a degree of freedom, N bJ e(s) its associated shape function and ns the number of elements in the discretisation of s.
  • The CBC approach has been used for the fields {ϕ, γ, p} and the CBD approach has been specifically used for the field γ when discontinuous strains are expected across the thickness.
  • In addition, CBC and CBD approaches have been compared against a truncated Taylor series expansion, as that in [43], denoted as Taylor-Expansion (TE) approach.

6. Numerical examples

  • The objective of this section is to demonstrate the applicability of the proposed shell formulation via a series of numerical examples, in which convex multi-variable electromechanical constitutive models, defined in the context of continuum formulations [1–4], will be considered.
  • In all the examples, a reconstruction of the continuum associated with the shell has been carried out at a post-processing level.
  • This reconstruction, based on the mapping x in equation (11), enables to show results not only in the mid surface of the shell but also across its thickness.

6.1. Bending actuators

  • This example considers the actuation device with geometry depicted in Figure 3. 6.1.1. Results for bending actuator configuration 1 Objective 2: The second objective is to test the performance of the formulation in scenarios characterised by the presence of discontinuities of the electric field distribution across the thickness of the shell.
  • Interestingly, Figures 7g−l show the purely mechanical and electrical contributions of the Cauchy stress tensor.
  • Regarding objective 2 and objective 3, the same conclusion as those obtained in the previous example are obtained and hence, omitted for brevity.

6.2. Helicoidal actuator

  • Regarding the boundary conditions, the degrees of freedom associated with the displacements of the mid surface of the shell and the director field d at X3 = 0m are completely constrained.
  • An electric charge per unit undeformed area of +ω0 and −ω0 is applied in both electrodes .
  • The value of the material parameters chosen for this particular example are shown in Table 2.
  • The first objective of this example is to demonstrate the applicability of the proposed formulation to scenarios where the reference configuration of the shell is curved, as that described by the cylindrically parametrised geometry (in the reference configuration) in equation (55), also known as Objective 1.
  • Figure 13 shows the contour plot of various stress and electric-like fields for a fixed value of the applied electric charge ω0.

6.3. Hyperboloid piezoelectric polymer

  • The hyperboloid with geometry described in Figure 15, presented in the context of pure elasticity in Reference [12], has been considered.
  • The material is transversely anisotropic, with the preferred axis of anisotropy N tangent to the surface of the hyperboloid as depicted in Figure 15.
  • The objective of this following example is to demonstrate the applicability of the proposed shell formulation to piezoelectric materials, where deformations can create a distribution of electric field in the material.
  • 8The area expansion has been computed as 1/γ, with γ the thickness stretch.
  • 34 35 Figures 16 displays contour plot of the (mechanically induced) electric field E3 for different values of the applied surface force q. Finally, Figure 17 shows the contour plot distribution of H22, σ33, p, ϕ, E1 and D03 for a given value of the applied surface force q.

7. Concluding remarks

  • This paper has provided a computational approach to formulate incompressible EAPs shells undergoing large strains and large electric field scenarios.
  • The proposed formulation, based upon a rotationless kinematical description of the shell, stems from the variational and constitutive framework proposed by the authors in previous publications [1–4], degenerated in this paper to the case of a nonlinear shell theory.
  • Two approaches have been considered for the interpolation of the electric potential across the thickness of the shell.
  • Specifically, the continuumbased-continuous (CBC) approach described in Section 5.2 and the Taylor expansion approach (TE) in [43].
  • A comparison of the results rendered by both approaches has been presented.

Did you find this useful? Give us your feedback

Figures (14)
Citations
More filters
Journal Article•DOI•
A reduced mixed finite-element formulation for modeling the viscoelastic response of electro-active polymers at finite deformation:

[...]

Dana Bishara1, Mahmood Jabareen1•
Technion – Israel Institute of Technology1
01 May 2019-Mathematics and Mechanics of Solids
TL;DR: In this article, a parameter identification procedure has been held for characterizing the widely used dielectric elastomer VHB, which has been performed using various experimental setups and parameters.
Abstract: In this work, a parameter identification procedure has been held for characterizing the widely used dielectric elastomer VHB. The calibration procedure has been performed using various experimental...

18 citations

Journal Article•DOI•
A mixed variational framework for the design of energy–momentum integration schemes based on convex multi-variable electro-elastodynamics

[...]

Marlon Franke1, Rogelio Ortigosa2, Alexander Janz1, Antonio J. Gil3, Peter Betsch1 •
Karlsruhe Institute of Technology1, University of Cartagena2, Swansea University3
01 Jul 2019-Computer Methods in Applied Mechanics and Engineering
TL;DR: The objective of the present paper is the introduction of new mixed variational principles for EM time integrators in electromechanics, hence bridging the gap between the previous work presented by the authors and Betsch et al. (2018), opening up the possibility to a variety of new Finite Element implementations.
Abstract: In Ortigosa et al. (2018), the authors presented a new family of time integrators for large deformation electromechanics. In that paper, definition of appropriate algorithmic expressions for the discrete derivatives of the internal energy and consideration of multi-variable convexity of the internal energy was made. These two ingredients were essential for the definition of a new energy–momentum (EM) time integrator in the context of large deformation electromechanics relying on materially stable (ellipticity compliant) constitutive models. In Betsch et al. (2018), the authors introduced a family of EM time integrators making use of mixed variational principles for large strain mechanics. In addition to the displacement field, the right Cauchy–Green deformation tensor , its co-factor and its Jacobian were introduced as unknown fields in the formulation. An elegant cascade system of kinematic constraints was introduced in this paper, crucial for the satisfaction of the required conservation properties of the new family of EM time integrators. The objective of the present paper is the introduction of new mixed variational principles for EM time integrators in electromechanics, hence bridging the gap between the previous work presented by the authors in Ortigosa et al. (2018) and Betsch et al. (2018), opening up the possibility to a variety of new Finite Element implementations. The following characteristics of the proposed EM time integrator make it very appealing: (i) the new family of time integrators can be shown to be thermodynamically consistent and second order accurate; (ii) piecewise discontinuous interpolation of the unknown fields (except displacements and electric potential) has been carried out, in order to yield a computational cost comparable to that of standard displacement–potential formulations. Finally, a series of numerical examples are included in order to demonstrate the robustness and conservation properties of the proposed scheme, specifically in the case of long-term simulations.

10 citations

Journal Article•DOI•
A complete direct approach to nonlinear modeling of dielectric elastomer plates

[...]

Elisabeth Hansy-Staudigl1, Michael Krommer1, Alexander Humer•
Vienna University of Technology1
01 Nov 2019-Acta Mechanica
TL;DR: In this article, a direct approach to nonlinear modeling of thin plates, which are made of incompressible dielectric elastomers, is presented, where the effect of electrostatic forces is incorporated by the purely electrical contribution to the augmented Helmholtz free energy.
Abstract: In this paper, we present a complete direct approach to nonlinear modeling of thin plates, which are made of incompressible dielectric elastomers. In particular, the dielectric elastomers are assumed to exhibit a neo-Hookean elastic behavior, and the effect of electrostatic forces is incorporated by the purely electrical contribution to the augmented Helmholtz free energy. Our approach does not involve any extraction-type procedure from the three-dimensional energy to derive the plate augmented free energy, but directly postulates the form of this energy for the structural plate problem treated in this paper. Results computed within the framework of this novel approach are compared to results available in the literature as well as to our own three-dimensional finite element solutions. A very good agreement is found.

9 citations

Journal Article•DOI•
A new stabilisation approach for level-set based topology optimisation of hyperelastic materials

[...]

Rogelio Ortigosa1, Rogelio Ortigosa2, Jesús Martínez-Frutos1, Antonio J. Gil2, D. Herrero-Pérez1 •
University of Cartagena1, Swansea University2
01 Jan 2019-Structural and Multidisciplinary Optimization
TL;DR: The robustness and applicability of this new methodological approach are thoroughly demonstrated through an ample spectrum of challenging numerical examples, ranging from benchmark two-dimensional (plane stress) examples to larger scale three-dimensional applications.
Abstract: This paper introduces a novel computational approach for level-set based topology optimisation of hyperelastic materials at large strains This, to date, is considered an unresolved open problem in topology optimisation due to its extremely challenging nature Two computational strategies have been proposed to address this problem The first strategy resorts to an arc-length in the pre-buckling region of intermediate topology optimisation (TO) iterations where numerical difficulties arise (associated with nucleation, disconnected elements, etc), and is then continued by a novel regularisation technique in the post-buckling region In the second strategy, the regularisation technique is used for the entire loading process at each TO iteration The success of both rests on the combination of three distinct key ingredients First, the nonlinear equilibrium equations of motion are solved in a consistent incrementally linearised fashion by splitting the design load into a number of load increments Second, the resulting linearised tangent elasticity tensor is stabilised (regularised) in order to prevent its loss of positive definiteness and, thus, avoid the loss of convexity of the discrete tangent operator Third, and with the purpose of avoiding excessive numerical stabilisation, a scalar degradation function is applied on the regularised linearised elasticity tensor, based on a novel regularisation indicator field The robustness and applicability of this new methodological approach are thoroughly demonstrated through an ample spectrum of challenging numerical examples, ranging from benchmark two-dimensional (plane stress) examples to larger scale three-dimensional applications Crucially, the performance of all the designs has been tested at a post-processing stage without adding any source of artificial stiffness Specifically, an arc-length Newton-Raphson method has been employed in conjunction with a ratio of the material parameters for void and solid regions of 10− 12

8 citations

Cites background from "A computational framework for incom..."

  • ...The reader is referred to Ortigosa et al (2015); Ortigosa and Gil (2017) for a detailed derivation of C in (11)....

    [...]

Journal Article•DOI•
Isogeometric MITC shell

[...]

Yongzhen Mi1, Xiang Yu2•
Shanghai Jiao Tong University1, Institute of High Performance Computing Singapore2
15 Apr 2021-Computer Methods in Applied Mechanics and Engineering
TL;DR: In this paper, an isogeometric formulation of the Reissner-Mindlin shell, using the Mixed Interpolation of Tensorial Components (MITC) technique to alleviate shear locking and membrane locking, is proposed.
Abstract: This paper proposes an isogeometric formulation of the Reissner–Mindlin shell, using the Mixed Interpolation of Tensorial Components (MITC) technique to alleviate shear locking and membrane locking. Instead of over each Non-Uniform Rational B-Spline (NURBS) element, kinematics of the shell is directly formulated on the entire NURBS patch, with its displacements decoupled to the motions of control points and of director vectors tied to the control points. Since control points are in general not located on the surface, those control-point director vectors are interpolated from normal vectors at a set of pre-defined integration points on the patch. The assumed in-plane membrane strain field and the assumed transverse shear strain field are built as linear combinations of NURBS basis functions with degrees lower than those employed by the displacement interpolation, and tied to the original covariant strain fields at a set of well-selected tying points. Integration points for the definition of control-point director vectors and tying points of the assumed covariant strain fields are provided by the optimal and reduced macro-element quadrature rules, respectively. In this way, the classical MITC technique is fully incorporated into the isogeometric framework with arbitrary polynomial orders or patch configurations. Thanks to the high smoothness of NURBS basis functions, the geometry-preserving nature of isogeometric transformation, and the locking-resistive capability of MITC technique, the isogeometric MITC shell formulation shows superior convergence behavior, minimal geometrical error, and good robustness against patch distortion. In particular, its convergence property is insensitive to NURBS polynomial order or control point density. These advantages are demonstrated through a number of well-established benchmark problems, validating the proposed isogeometric MITC formulation being an accurate and efficient solution for linear analysis of shell structures .

6 citations

References
More filters
Journal Article•DOI•
High-Speed Electrically Actuated Elastomers with Strain Greater Than 100%

[...]

Ronald E. Pelrine1, Roy D. Kornbluh1, Qibing Pei1, Jose P. Joseph1•
SRI International1
04 Feb 2000-Science
TL;DR: It is shown that prestraining the film further improves the performance of electrical actuators made from films of dielectric elastomers coated on both sides with compliant electrode material.
Abstract: Electrical actuators were made from films of dielectric elastomers (such as silicones) coated on both sides with compliant electrode material. When voltage was applied, the resulting electrostatic forces compressed the film in thickness and expanded it in area, producing strains up to 30 to 40%. It is now shown that prestraining the film further improves the performance of these devices. Actuated strains up to 117% were demonstrated with silicone elastomers, and up to 215% with acrylic elastomers using biaxially and uniaxially prestrained films. The strain, pressure, and response time of silicone exceeded those of natural muscle; specific energy densities greatly exceeded those of other field-actuated materials. Because the actuation mechanism is faster than in other high-strain electroactive polymers, this technology may be suitable for diverse applications.

2,969 citations

Journal Article•DOI•
Convexity conditions and existence theorems in nonlinear elasticity

[...]

John M. Ball1•
Heriot-Watt University1
01 Dec 1976-Archive for Rational Mechanics and Analysis

2,329 citations

"A computational framework for incom..." refers background in this paper

  • ...The aforementioned frameworks are based on the use of convex multi-variable strain energies, which represent an extension of the concept of polyconvexity [5] to the field of nonlinear electro-mechanics....

    [...]

Book•
Nonlinear Continuum Mechanics for Finite Element Analysis

[...]

Javier Bonet1, Richard D. Wood1•
Swansea University1
28 Sep 1997
TL;DR: Bonet and Wood as discussed by the authors provide a complete, clear, and unified treatment of nonlinear continuum analysis and finite element techniques under one roof, providing an essential resource for postgraduates studying non-linear continuum mechanics and ideal for those in industry requiring an appreciation of the way in which their computer simulation programs work.
Abstract: Designing engineering components that make optimal use of materials requires consideration of the nonlinear characteristics associated with both manufacturing and working environments. The modeling of these characteristics can only be done through numerical formulation and simulation, and this requires an understanding of both the theoretical background and associated computer solution techniques. By presenting both nonlinear continuum analysis and associated finite element techniques under one roof, Bonet and Wood provide, in this edition of this successful text, a complete, clear, and unified treatment of these important subjects. New chapters dealing with hyperelastic plastic behavior are included, and the authors have thoroughly updated the FLagSHyP program, freely accessible at www.flagshyp.com. Worked examples and exercises complete each chapter, making the text an essential resource for postgraduates studying nonlinear continuum mechanics. It is also ideal for those in industry requiring an appreciation of the way in which their computer simulation programs work.

1,859 citations

Journal Article•DOI•
Electrostriction of polymer dielectrics with compliant electrodes as a means of actuation

[...]

Ronald E. Pelrine1, Roy D. Kornbluh1, Jose P. Joseph1•
SRI International1
01 Jan 1998-Sensors and Actuators A-physical
TL;DR: Electrostrictive polymer (EP) dielectric actuators have been shown to produce 5 to 20 times the effective actuation pressure of conventional air-gap electrostatics at the same electric field strength as mentioned in this paper.
Abstract: The electrostriction of elastomeric polymer dielectrics with compliant electrodes is potentially useful as a small-scale, solid-state actuator technology. Electrostrictive polymer (EP) materials are capable of efficient and fast response with high strains (> 30%), good actuation pressures (up to 1.9 MPa), and high specific energy densities (up to 0.1 J g−1). In this article, the mechanism of electrostriction is shown to be due to the electrostatic attraction of free charges on the electrodes. Although EP actuators are electrostatics based, they are shown to produce 5–20 times the effective actuation pressure of conventional air-gap electrostatics at the same electric field strength. The thin uniform dielectric films necessary for fabrication of EP actuators have been fabricated by techniques such as spin coating, casting, and dipping. A variety of materials and techniques have been used to produce the compliant electrodes, including lift-off stenciling techniques for powdered graphite, selective wetting of ionically conductive polymers, and spray coating of carbon blacks and fibrils in polymeric binders. Prototype actuators have been demonstrated in a variety of configurations such as stretched films, stacks, rolls, tubes, and unimorphs. Potential applications of the technology in areas such as microrobots, sound generators, and displays are discussed in this article.

1,436 citations

Journal Article•DOI•
Analysis of thick and thin shell structures by curved finite elements

[...]

Sohrabuddin Ahmad1, Bruce M. Irons2, O. C. Zienkiewicz2•
University of Engineering and Technology, Lahore1, Swansea University2
01 Jul 1970-International Journal for Numerical Methods in Engineering
TL;DR: In this article, a general formulation for the curved, arbitrary shape of thick shell finite elements is presented along with a simplified form for axisymmetric situations, which is suitable for thin to thick shell applications.
Abstract: A general formulation for the curved, arbitrary shape of thick shell finite elements is presented in this paper along with a simplified form for axisymmetric situations. A number of examples ranging from thin to thick shell applications are given, which include a cooling tower, water tanks, an idealized arch dam and an actual arch dam with deformable foundation. A new process using curved, thick shell finite elements is developed overcoming the previous approximations to the geometry of the structure and the neglect of shear deformation. A general formulation for a curved, arbitrary shape of shell is developed as well as a simplified form suitable for axisymmetric situations. Several illustrated examples ranging from thin to thick shell applications are given to assess the accuracy of solution attainable. These examples include a cooling tower, tanks, and an idealized dam for which many alternative solutions were used. The usefulness of the development in the context of arch dams, where a ‘thick shell’ situation exists, leads in practice to a fuller discussion of problems of foundation deformation, etc., so that practical application becomes possible and economical.

1,205 citations

Related Papers (5)
A new framework for large strain electromechanics based on convex multi-variable strain energies: Finite Element discretisation and computational implementation

[...]

15 Apr 2016-Computer Methods in Applied Mechanics and Engineering
Rogelio Ortigosa, Antonio J. Gil
A new framework for large strain electromechanics based on convex multi-variable strain energies: Variational formulation and material characterisation

[...]

15 Apr 2016-Computer Methods in Applied Mechanics and Engineering
Antonio J. Gil, Rogelio Ortigosa
A computational framework for polyconvex large strain elasticity

[...]

01 Jan 2015-Computer Methods in Applied Mechanics and Engineering
Javier Bonet, Antonio J. Gil, Rogelio Ortigosa
Electrostatic Forces and Stored Energy for Deformable Dielectric Materials

[...]

01 Jul 2005-Journal of Applied Mechanics
Robert M. McMeeking, Chad M. Landis
Numerical modelling of non-linear electroelasticity

[...]

07 May 2007-International Journal for Numerical Methods in Engineering
Duc Khoi Vu, Paul Steinmann, Gunnar Possart
Frequently Asked Questions (2)
Q1. What have the authors contributed in "A computational framework for incompressible electromechanics based on convex multi-variable strain energies for geometrically exact shell theory" ?

In this paper, a new computational framework for the analysis of incompressible Electro Active Polymer ( EAP ) shells subjected to large strains and large electric fields is presented. Two novelties are incorporated in this work. First, the variational and constitutive frameworks developed by the authors in recent publications [ 1–4 ] in the context of three-dimensional electromechanics are particularised/degenerated to the case of geometrically exact shell theory. The proposed formulation follows a rotationless description of the kinematics of the shell, enhanced with extra degrees of freedom corresponding to the thickness stretch and the hydrostatic pressure, critical for the consideration of incompressibility. More specifically, convex multi-variable ( three-dimensional ) constitutive models, complying with the ellipticity condition and hence, satisfying material stability for the entire range of deformations and electric fields, Corresponding author: r. ortigosa @ swansea. Different approaches are investigated for the interpolation of these extra fields and that of the electric potential across the thickness of the shell. 

Moreover, the kinematics of the shell allows for the possibility of compression and stretch across the thickness of the shell [ 17 ], crucial for the consideration of incompressible behaviour. Two approaches have been considered for the interpolation of the electric potential across the thickness of the shell.