Fig. 3.2.1 Finite Element Mesh Consisting of Triangular and Rectangular Element
Fig. 3.2.1 Finite Element Mesh Consisting of Triangular and Rectangular Element
Fig. 3.2.1 Finite Element Mesh Consisting of Triangular and Rectangular Element
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
10
(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)
(3.2.4)
.....
.....
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
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.
16
(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)
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.