**Mapping and Shape Functions**

**Introduction**

A finite element is called 'isoparametric' if the interpolation functions describing
its geometry are identical to the interpolation functions describing the distribution of displacements.

In many FEA software packages the coordinate system of an element is transformed into a coordinate
system called the natural coordinate system. In order to avoid mathematical difficulties
element distortion must be controlled. This approach simplifies the addition of new element
types to a software package.

Natural coordinates can be used to define variation of a quantity along an element, in the
above case, elastic displacement. Assuming elastic displacement, paricle displacement varies linearly between
points (nodes) 1 and 2.

This linear variation is called an interpolation function.

In (b): r = 0A/02

From (a): r = [x - ( x_{2} + x_{1} )/2] / [( x_{2} - x_{1} )/2] .. eqn (a1)

It is assumed that the displacement, u, varies linearly with r:

Eqn (c): u = a_{1} + a_{2} r

At node 1: u = u_{1} and r = -1, so u_{1} = a_{1} - a_{2}

At node 2: u = u_{2} and r = 1, so u_{2} = a_{1} + a_{2}

Adding gives: u_{1} + u_{2} = 2 a_{1}, or a_{1} = (u_{1} + u_{2})/2

Subtracting: u_{2} - u_{1} = 2 a_{2}, or a_{2} = (u_{2} - u_{1})/2

Substituting in (c) for a_{1} and a_{2}:

N

N

At any value of r, N_{1} + N_{} = (1 - r)/2 + (1 + r)/2 = 1

The co-ordinates of a point in the element (master) can be related to co-ordinates in the real element by re-arranging eqn (a1):

or: x = [(1 - r)/2]x_{1} + [(1 + r)/2]x_{2}

or: x = N_{1}x_{1} + N_{2}x_{2} ... This is a 'mapping relationship'

Having linear variation of displacement in an element which means constant stress is very limiting and it is much more useful to allow linear variation of stress within an element which requires higher order polynomial for interpolation, i.e. a quadratic 3 node linear element for the example above (and see in the section: Example on reducing errors ..).

at r = 0, u = u_{3}, so u_{3} = a_{1}

at r = 1, u = u_{2}, so u_{2} = a_{1} + a_{2} + a_{3}

Re-arranging gives:

Where N_{1}, N_{2} and N_{3} are the shape functions.

N_{1} = -r(1 - r)/2

N_{2} = r(1 + r)/2

N_{3} = (1 + r)(1 - r)

These are shown plotted below:

For any value of r between -1 and 1, the sum of the shape functions, in this case: N

This procedure means that 1 dimensional, 2D and 3D elements which in the real world may have 1 or more curved edges, can be mapped onto the straight line r; the square r, s; or the cube r, s, t.

David Grieve, 16th November 2012