The Lagrange element, or nodal elment, was first defined with use of Lagrange interpolation polynomials in Courant's introduction of the finite elment method [1].
1D Lagrange elements order 1 The first element with nodes at 0 and .5. A second with nodes at .5 and 1.0 |
1D Lagrange elements order 2 with nodes at 0, .5, and 1.0 |
Order 2 Triangular Lagrange |
Order 3 Triangular Lagrange |
Element definition
The approximating basis functions are determined by the Kronecker delta function.
(1)The element can extended the polynomial definition to arbitrary order and requires only an affine transformation. The requirement of nodes is determined by the necessary order, see examples below for more details.
Examples
One dimension
One dimensional basis functions are the same as the Lagrange polynomials. For each segment there will be $k+1$ nodes where $k$ is the order of the polynomial space. The basis functions will be:
(2)Simplices
For simplices, each edge and face must have the correct number of nodes for the space, and the additional nodes place anywhere within the cell. Nodes on are typically put on a lattice over the element. A common way to tabulate the necessary basis functions is using a vandermonde matrix.
For example using a 2nd order triangular element represented functions will be of the form:
(3)which will require 6 nodes (or degrees of freedom). Since each edge needs 3 nodes, all the node will be on the edges. For the reference element we have the following nodes:
Vertices | (0,0) | (1,0) | (0,1) |
Edge-node | (.5,0) | (.5,.5) | (0,.5) |
Thus the vandermonde system will be:
(4)where $\delta_i$ is 1 in the $i$th component and 0 otherwise. Solving for $i=0...5$ produces the following basis functions:
$i$ | $\phi_i(x,y)$ |
---|---|
0 | $1 - 3x - 3y + 2x^2 + 4xy + 2y^2$ |
1 | $-x + 2x^2$ |
2 | $-y + 2y^2$ |
3 | $4x - 4x^2 - 4xy$ |
4 | $4xy$ |
5 | $4y - 4xy - 4y^2$ |
A similar process may be applied for tetrahedra. The basis polynomials go as:
Triangles | Tetrahedra | |||
---|---|---|---|---|
order | polynomial | number of nodes | polynomial | number of nodes |
0 | $1$ | 1 | $1$ | 1 |
1 | $1 + x + y$ | 1 + 2 | $1 + x + y + z$ | 1 + 3 |
2 | $1 + x + y + x^2 + xy + y^2$ | 1 + 2 + 3 | $1 + x + y + x^2 + xy + xz+ y^2 + yz + z^2$ | 1 + 3 + 5 |
k | … | 1 + 2 + … + $(k+1)$ | … | 1 + 3 + … + $(k+2)(k+1)/2$ |
Quadralateral and hexahedrals
Basis functions for quadralaterals and hexahedrals can be formed the same way as the simplex counterpart. For lower order function spaces the number of nodes required on each edge $k+1$ is more than required in total, for quadralateral of order 1 there are 4 nodes defined by the edges but only 3 needed to define the space. To resolve this additional higher order terms are added, for order 1 the $xy$ term is added.