Fig. 3.2.1 Finite Element Mesh Consisting of Triangular and Rectangular Element

Download as pdf or txt
Download as pdf or txt
You are on page 1of 9

Module 3: Element Properties

Lecture 2: Triangular Elements


9

The triangular element can be used to represent the arbitrary geometry much easily. On the other
hand, rectangular elements, in general, are of limited use as they are not well suited for representing
curved boundaries. However, an assemblage of rectangular and triangular element with triangular
elements near the boundary can be very effective (Fig. 3.2.1). Triangular elements may also be used
in 3-dimensionalaxi-symmetric problems, plates and shell structures. The shape function for
triangular elements (linear, quadratic and cubic) with various nodes (Fig. 3.2.2) can be formulated.
An internal node will exist for cubic element as seen in Fig. 3.2.2(c).

Fig. 3.2.1 Finite element mesh consisting of triangular and rectangular element

Fig. 3.2.2 Triangular elements


In displacement formulation, it is very important to approximate the variation of displacement in the
element by suitable function. The interpolation function can be derived either using the Cartesian
coordinate system or by the area coordinates.

10

3.2.1 Shape function usingCartesiancoordinates


Polynomials are easiest way of mathematical operation for expressing variation of displacement. For
example, the displacement variation within the element can be represented by the following function
in case of two dimensional plane stress/strain problems.
u = + x + y

(3.2.1)

v= + x+ y

(3.2.2)

where 0, 1, 2 .. are unknown coefficients. Thus the displacement vectors at any point P,in the
element (Fig.3.2.3) can be expressed with the following relation.

1
0

0
1

(3.2.3)

Or, {d}=[ ]{}

(3.2.4)

Fig. 3.2.3 Triangular element in Cartesian Coordinates


Similarly, for m node element having three degrees of freedom at each node, the displacement
function can be expressed as
. . . . .
(3.2.5)

.....
.....
Hence,in such case,

11

0
0
Where,

0
and,

. . .

0
0

(3.2.6)
1

. . . .

Now, for a linear triangular element with 2 degrees of freedom, eq. (3.2.3) can be written in terms of
the nodal displacements as follows.
1
0 0 0
1
0 0 0
0 0 0
1
(3.2.7)
0 0 0 1
0 0 0 1
0 0 0 1
Where, {d} is the nodal displacements. To simplify the above expression for finding out the shape
function, the displacements in X direction can be separated out which will be as follows:

=
= 1
(3.2.8)

1
To obtain the polynomial coefficients, {} the matrix of the above equation are to be inverted. Thus,

a 0 1 x1


a1 = 1 x 2

a 1 x 3
2

a1 a 2
1
=
b1 b 2
2A
c1 c 2

y1

y2

y3

-1

u
x y - x 3 y2
1
2 3
1
y 2 - y3
u 2 =
2A
x3 - x2

u 3

x 3 y1 - x1 y3
y3 - y1
x1 - x 3

a 3 u1

b 3 u 2

c3

u 3

Where, A is the area of the triangle and can be obtained as follows.


1
1
A
1
Now, eq. (3.2.1) can be written from the above polynomial coefficients.

x1 y 2 - x 2 y1 u1

y1 - y2 u 2


x 2 - x1
u 3

(3.2.9)

(3.2.10)

12

1
( x 2 y3 - x 3 y 2 ) + ( y 2 - y3 ) x + ( x 3 - x 2 ) y u1
2A
1
+
( x 3 y1 - x1y3 ) + ( y3 - y1 ) x + ( x1 - x 3 ) y u 2
2A
1
+
( x1y 2 - x 2 y1 ) + ( y1 - y 2 ) x + ( x 2 - x1 ) y u 3
2A
Thus, the interpolation function can be obtained from the above as:
1

( x 2 y 3 - x 3 y 2 ) + ( y 2 - y 3 ) x + ( x 3 - x 2 ) y

2A

N1

T
1
{N } = N 2 =
( x 3 y1 - x1y 3 ) + ( y 3 - y1 ) x + ( x1 - x 3 ) y

2A

x
y
x
y
y
y
x
x
x
y
+
+
(
)
(
)
(
)

2 1
1
2
2
1
1 2

2A

u=

(3.2.11)

(3.2.12)
Such three node triangular element is commonly known as constant strain triangle (CST) as its strain
is assumed to be constant inside the element. This property may be derived from eq. (3.2.1) and
eq.(3.2.2). For example, in case of 2-D plane stress/strain problem, one can express the strain inside
the triangle with the help of eq.(3.2.1) and eq.(3.2.2):

u (a 0 + a1x + a 2 y)
=
= a1
x
x
v (a 3 + a 4 x + a 5 y)
ey =
=
= a5
y
y
v u
g xy =
+
= a2 + a4
x y
ex =

(3.2.13)

CST is the simplest element to develop mathematically. As there is no variation of strain inside the
element, the mesh size of the triangular element should be small enough to get correct results. This
element produces constant temperature gradients ensuring constant heat flow within the element for
heat transfer problems.
3.2.2 Higher Order Triangular Elements
Higher order elements are useful if the boundary of the geometry is curve in nature. For curved case,
higher order triangular element can be suited effectively while generating the finite element mesh.
Moreover, in case of flexural action in the member, higher order elements can produce more
accurate results compare to those using linear elements. Various types of higher order triangular

13
elements are used in practice. However, most commonly used triangular element is the six node
element for which development of shape functions are explained below.
3.2.2.1Shape function for six node element
Fig. 3.2.4 shows a triangular element with six nodes. The additional three nodes (4, 5, and 6) are
situated at the midpoints of the sides of the element. A complete polynomial representation of the
field variable can be expressed with the help of Pascal triangle:

x, y 0 1 x 2 y 3 x2 4 xy 5 y 2

(3.2.14)

Fig. 3.2.4 (a) Six node triangular element (b) Lines of constant values of the area coordinates
Using the above field variable function, one can reach the following expression using interpolation
function and the nodal values.
6

f( x, y) = Ni ( x, y)fi
i=1

(3.2.15)

Here, the every shape function must be such that its value will be unity if evaluated at its related
node and zero if evaluated at any of the other five nodes. Moreover, as the field variable
representation is quadratic, each interpolation function will also become quadratic. Fig. 3.2.4(a)
shows the six node element with node numbering convention along with the area coordinates at three
corners. The six node element with lines of constant values of the area coordinates passing through
the nodes is shown in Fig. 3.2.4(b). Now the interpolation functions can be constructed with the help
of area coordinates from the above diagram. For example, the interpolation function N1 should be
unity at node 1 and zero at all other five nodes. According to the above diagram, the value of L1 is 1

14
at node 1 and at node 4 and 6. Again, L1 will be 0 at nodes 2, 3 and 5. To satisfy all these
conditions, one can propose following expression:

1
N1 ( x, y) = N1 (L1 ,L2 ,L3 ) = L1 L1 -

(3.2.16)

Evaluating the above expression, the value of N1 is becoming at node 1 though it must become
unity. Therefore, the above expression is slightly modified satisfying all the conditions and will be as
follows:

1
N1 = 2L1 L1 - = L1 (2L1 -1)
(3.2.17)

2
Eq. (3.2.17) assures the required conditions at all the six nodes and is a quadraticfunction, asL1 is a
linear function of x and y. The remaining five interpolationfunctions can also be obtained in similar
fashion applying the required nodal conditions. Thus, the shape function for the six node triangle
element can be written as given below.
N1 = L1 (2L1 -1)
N 2 = L2 (2L2 -1)
N3 = L3 (2L3 -1)

(3.2.18)

N 4 = 4L1L2
N5 = 4L2L3
N6 = 4L3L1

Such six node triangular element is commonly known as linear strain triangle (LST) as its strain is
assumed to vary linearly inside the element. In case of 2-D plane stress/strain problem, the element
displacement field for such quadratic triangle may be expressed as

u x, y 0 1 x 2 y 3 x2 4 xy 5 y 2
v x, y 6 7 x 8 y 9 x2 10 xy 11 y 2

(3.2.19)

So the element strain can be derived from the above displacement field as follows.

u
= a1 + 2a 3 x + a 4 y
x
v
ey =
= a 8 + a10 x + 2a11y
y
v u
g xy =
+
= a 2 + a 4 x + 2a 5 y + a 7 + 2a 9 x + a10 y
x y
ex =

(3.2.20)

The above expression shows that the strain components are linearly varying inside the element.
Therefore, this six node element is called linear strain triangle. The main advantage of this element is
that it can capture the variation of strains and therefore stresses of the element.

15
3.2.3Construction of Shape Function by Degrading Technique
Sometimes, the geometry of the structure or its loading and boundary conditions are such that the
stresses developed in few locations are quite high. On the other hand, variations of stresses are less
in some areas and as a result, refinement of finite element mesh is not necessary. It would be
economical in terms of computation if higher order elements are chosen where stress concentration
is high and lower order elements at area away from the critical area. Fig.3.2.5 shows graphical
representationswhere various order of triangular elements are used for generating a finite element
mesh.

Fig. 3.2.5 Triangular elements with different number of nodes


Fig. 3.2.5contains four types of element. Type 1 has only three nodes, type 2 element has five nodes,
type 3 has four nodes and type 4 has six nodes. The shape function for 3-node and 6-node triangular
elements has already been derived. The shape functions of 6-node element can suitably be degraded
to derive shape functions of other two types of triangular elements.
3.2.3.1 Five node triangular element
Let consider a six node triangular element as shown in Fig. 3.2.6(a) whose shape functions and nodal
displacements are (N1, N2, N3, N4, N5, N6) and (u1, u2, u3, u4, u5, u6) respectively. Similarly, for a five
node triangular element as shown in Fig. 3.2.6(b), the shape functions and nodal displacements are
considered as (N1, N2, N3, N4, N5) and (u1, u2, u3, u4, u5) respectively. Thus, the
displacementat any point in a six node triangular element will become

16

u = N1u1 + N 2 u 2 + N3u 3 + N 4 u 4 + N5u 5 + N 6 u 6

(3.2.21)
Where, N1, N2, , N6 are the shape functions and is given in eq.(3.2.18).If there is no node between
1 and 3, the displacement along line 1-3is considered to vary linearly. Thus the displacement at an
assumed node6 may be written as

u '6 =

u1 + u 3
2

(3.2.22)

Substituting, the value of u5 for u5 in eq.(3.2.21) the following expression will be obtained.

u = N1u1 + N 2 u 2 + N 3u 3 + N 4 u 4 + N 5 u 5 + N 6

u1 + u 3
2

(3.2.23)

Fig. 3.2.6 Degrading for five node element


Thus, the displacement function can be expressed by five nodal displacements as:

N
N
u = N1 + 6 u1 + N 2 u 2 + N3 + 6 u 3 + N 4 u 4 + N 6 u 5

2
2

(3.2.24)

However, the displacement function for the five node triangular element can be expressed as

u = N1u 1 + N 2 u 2 + N 3 u 3 + N 4 u 4 + N 5 u 5

(3.2.25)
Comparing eq.(3.2.24) and eq.(3.2.25) and observing node 6 of six node triangle corresponds to
node 5 of five node triangle, we can write

N6
N
, N '2 = N2, N3 = N3 + 6 , N4 = N4 and N5 = N5
2
2
Hence, the shape function of a five node triangular element will be
N1 = N1 +

(3.2.26)

17

N1 = N1 +

N6
4L L
= L1 (2L1 -1) + 1 3 = L1 (1- 2L 2 )
2
2

N 2 = L 2 (2L 2 -1)
N
4L L
N3 = N3 + 6 = L3 (2L3 -1) + 1 3 = L3 (1- 2L 2 )
2
2
N 4 = 4L1L 2

(3.2.27)

N5 = 4L 2 L3
Thus, for a five node triangular element, the above shape function can be used for finite element
analysis.

You might also like