Four Node Iso-Parametric Quadrilateral Element
Four Node Iso-Parametric Quadrilateral Element
(
= = ( ` ` `
( ) )
)
J
J
where |J| is the Jacobian (or determinant) and is defined by
x y x y
s t t s
=
J
- Since the Jacobian appears in the denominator in the above equation, it must not be
zero anywhere over the domain (1 s, t 1).
- The mapping is not valid if |J| is zero anywhere over the element.
Examples
(1) Check the validity of iso-parametric mapping for the element shown in Figure 2
3
1(0, 0) 2(1, 0)
3(2, 2)
4(0, 1)
x
y
Figure 2 Four-node quadrilateral element
- Nodal coordinates:
1 2 3 4
1 2 3 4
0, 1, 2, 0
0, 1, 2, 1
x x x x
y y y y
= = = =
= = = =
- Iso-parametric mapping:
1 1 2 2 3 3 4 4 2 3
1 1 2 2 3 3 4 4 3 4
1
2 (3 3 )
4
1
2 (3 3 )
4
x N x N x N x N x N N s t st
y N y N y N y N y N N s t st
= + + + = + = + + +
= + + + = + = + + +
- Jacobian:
3 1
1
[ ]
1 3 4
x y
t t
s s
x y s s
t t
(
(
+ + (
= = (
(
+ +
(
(
J
1 1 1 1
[(3 )(3 ) (1 )(1 )]
4 2 8 8
t s t s s t = + + + + = + + J
Thus, it is clear that |J| > 0 for 1 s 1 and 1 t 1.
Figure 3 Element shape obtained from the iso-parametric mapping
4
(2) Check the validity of iso-parametric mapping for the element shown in Figure 4.
1(0, 0)
2(1, 5)
3(6, 6) 4(0, 6)
x
y
Figure 4 Four-node quadrilateral element
- Nodal coordinates:
1 2 3 4
1 2 3 4
0, 1, 6, 0
0, 5, 6, 6
x x x x
y y y y
= = = =
= = = =
- Iso-parametric mapping:
1 1 2 2 3 3 4 4
1 1 2 2 3 3 4 4
1
(1 )(7 5 )
4
1
(17 5 7 5 )
4
x N x N x N x N x s t
y N y N y N y N y s t st
= + + + = + +
= + + + = + +
- Jacobian:
3 15 15
2 4 4
s t
= + J
Note that |J| = 0 at 3/2 15s/4 + 15t/4, or s t = 2/5. The mapping illustrated in Figure 5
clearly shows that the mapping is invalid.
5
Figure 5 Element shape obtained from the iso-parametric mapping
- In general the element geometry is invalid if the Jacobian is either zero or negative
anywhere in the element.
- Problem also arises when the Jacobian matrix is nearly singular either due to round-
off errors or due to badly shaped elements.
- To avoid problems due to badly shaped elements, it is recommended that the inside
angles in quadrilateral elements be > 15 and < 165 as illustrated in Figure 6.
> 15
o
< 165
o
Figure 6 Recommended ranges of internal angles in a quadrilateral element
Interpolation of Displacement
- Iso-parametric mapping method in Section 6.5.1 can be used for displacement
interpolation.
- A total of eight degrees-of-freedom is used.
6
1
1
2
1 2 3 4 2
1 2 3 4 3
3
4
4
0 0 0 0
[ ]{ }
0 0 0 0
u
v
u
N N N N v u
N N N N u v
v
u
v
(
= =
` `
(
)
)
N d
where the same shape functions in geometric interpolation can be used for displacement
interpolation:
1
2
3
4
1
(1 )(1 )
4
1
(1 )(1 )
4
1
(1 )(1 )
4
1
(1 )(1 )
4
N s t
N s t
N s t
N s t
= +
= + +
= +
- Jacobian Matrix
[ ] ,
x y
x y x y
s s
x y s t t s
t t
(
(
= = (
(
(
J J
where
1 1
2 2
1 2 3 4
3 3
4 4
1
[ ] [ 1 1 1 1 ]
4
x x
x x
x N N N N
t t t t
x x s s s s s
x x
= = + +
` `
) )
1
2
3
4
1
[ 1 1 1 1 ]
4
x
x
x
s s s s
x t
x
= + +
`
)
Similar expressions can be written for derivatives of y with respect to s and t. In general,
the determinant of the Jacobian matrix can be expressed as
7
1
2
1 2 3 4
3
4
0 1 1
1 0 1
1
[ ]
1 0 1 8
1 1 0
T
y t s t s
y t s s t
x x x x
y s t s t
y s s t t
+ + (
(
+ +
(
=
`
( +
(
+
)
J
- Displacement-strain relationship
/
/ 1 0 0 0
/
{ } / 0 0 0 1
/
/ / 0 1 1 0
/
xx
yy
xy
u x
u x
u y
v y
v x
u y v x
v y
(
(
= = =
` ` `
(
+ (
)
)
)
1
u
y y u
x
t s s
u
x x u
y
t s t
(
(
= ( ` `
( )
)
J
1
v
y y v
x
t s s
v
x x v
y
t s t
(
(
= ( ` `
( )
)
J
Writing the two equations together:
/ / / 0 0 /
/ / / 0 0 /
1
/ 0 0 / / /
/ 0 0 / / /
u x y t y s u s
u y x t x s u t
v x y t y s v s
v y x t x s v t
(
(
(
=
` `
(
(
) )
J
The strain can now be expressed as
/ / 0 0 / /
1 0 0 0
/ / 0 0 / /
1
0 0 0 1 [ ]
0 0 / / / /
0 1 1 0
0 0 / / / /
xx
yy
xy
y t y s u s u s
x t x s u t u t
y t y s v s v s
x t x s v t v t
(
(
(
(
(
=
` ` `
(
(
(
(
)
) )
A
J
The derivative of the trial solution with respect to s and t are easy to compute:
8
1
1
2
2
3
3
4
4
/ 1 0 1 0 1 0 1 0
/ 1 0 1 0 1 0 1 0
1
[ ]{ }
/ 0 1 0 1 0 1 0 1 4
/ 0 1 0 1 0 1 0 1
u
v
u u s t t t t
v u t s s s s
u v s t t t t
v v t s s s s
u
v
+ + (
(
+ +
(
=
` `
+ + (
(
+ +
)
)
G d
The displacement-strain matrix [B] can now be written as follows:
/
/
[ ] [ ][ ]{ } [ ]{ }
/
/
xx
yy
xy
u s
u t
v s
v t
= =
` `
)
)
A A G d B d
- The displacement-strain matrix [B] is not constant. Thus, the value of strain and
stress within an element changes as a function of position.
Finite Element Matrix Equation
- The element stiffness matrix is:
1 1
1 1
[ ] [ ] [ ][ ] [ ] [ ][ ]
T T
A
h dA h dsdt
=
k B C B B C B J (0.1)
- Different from the triangular element, the element stiffness matrix is not constant
within an element. Thus, analytical integration is not trivial. Numerical integration
ethod using Gaussian quadrature will be discussed in the next section.
- Work done by concentrated forces:
1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4
{ } { }
T
NF x y x y x y x y NF
W F u F v F u F v F u F v F u F v = + + + + + + + d Q
where
1 1 2 2 3 3 4 4
{ } [ ]
T
NF x y x y x y x y
F F F F F F F F = Q is the vector of applied nodal
forces.
- Work done by distributed load:
{ } [ ] { } { } { }
T T T
T T
S
W h dS = =
d N T d Q
9
4
1
2
3
4
1 2
3
s
s
s
s
S
S
S
S
T
x
T
y
Figure 7 Coordinates for evaluation of boundary integrals
Side 1-2:
(1 ) / 2 0 (1 ) / 2 0 0 0 0 0
[ ] , 1 1
0 (1 ) / 2 0 (1 ) / 2 0 0 0 0
s s
s
s s
+
(
=
(
+
N
Side 2-3:
0 0 (1 ) / 2 0 (1 ) / 2 0 0 0
[ ] , 1 1
0 0 0 (1 ) / 2 0 (1 ) / 2 0 0
s s
s
s s
+
(
=
(
+
N
Side 3-4:
0 0 0 0 (1 ) / 2 0 (1 ) / 2 0
[ ] , 1 1
0 0 0 0 0 (1 ) / 2 0 (1 ) / 2
s s
s
s s
+
(
=
(
+
N
Side 4-1:
(1 ) / 2 0 0 0 0 0 (1 ) / 2 0
[ ] , 1 1
0 (1 ) / 2 0 0 0 0 0 (1 ) / 2
s s
s
s s
+
(
=
(
+
N
- As an example, consider side 2-3, as shown in Figure 7. The iso-parametric
mapping for side 2-3 is
2 3 2 3
1 1 1 1
,
2 2 2 2
s s s s
x x x y y y
+ +
= + = +
3 2 3 2
2 3
3 2 3 2
1 1
or
2 2 2 2
or
2 2
dx x x x x
x x dx ds
ds
dy y y y y
dy ds
ds
= + = =
= =
The Jacobian of the transformation from S to s is defined as
side2-3
dS
J
ds
=
with reference to Figure 8 it can be developed as follows. From geometry
10
-1
(x
2
, y
2
)
1
s
dy
S
dS
(x
3
, y
3
)
dx
Side 2-3 of parent element
Figure 8 Evaluation of work done by applied pressure along side 2-3
2 2 2 2
3 2 3 2 23
1 1
( ) ( )
2 2
dS dx dy x x y y ds L ds = + = + =
Thus, the Jacobian for each side is equal to half the length of that side. Thus,
1 1
,side 2-3 23 23
side 2 3 1 1
0 0
0 0
(1 ) / 2 0
0 (1 ) / 2
1 1
{ } [ ] { } [ ] { }
(1 ) / 2 0 2 2
0 (1 ) / 2
0 0
0 0
x
T T
T
y
s
T s
h dS h L ds h L ds
T s
s
(
(
(
(
(
(
= = =
`
(
+
)
(
+
(
(
(
(
Q N T N T
Carrying out the integration gives
23
,side2-3
{ } [0 0 0 0]
2
T
T x y x y
hL
T T T T = Q
- Similar to the linear triangular element, this equation says that total pressure along a
side is divided equally among the two nodes along the side.
- All quantities needed for the potential energy have now been expressed in terms of
nodal unknowns. From the principle of minimum total potential energy,
[ ]{ } { }
NT T
= + k d Q Q
where [k] is 88 element stiffness matrix, {d} is 81 nodal displacement vector, and
{Q
NT
+Q
T
} is 81 applied nodal force vector.
11
Numerical Integration
- Numerical integration evaluates the integrals involved in the element stiffness matrix
and distributed force.
- In the finite element literature, the Gauss quadrature is usually preferred because it
requires fewer function evaluations as compared to other methods.
- In the Gauss quadrature, the integrand is evaluated at predefined points (called Gauss
points). The sum of this integrand values, multiplied by appropriate weights (called
Gauss weight) gives an approximation to the integral:
1
1
1
( ) ( )
n
i i
i
I f s ds w f s
=
where n = number of Gauss points, s
i
= Gauss points, w
i
= Gauss weights, and f(s
i
) =
function value at the Gauss point s
i
.
- The locations of Gauss points and weights are derived in such a way that with n
points, a polynomial of degree 2n1 is integrated exactly.
- Table 1 shows the locations of the Gauss points and corresponding weights.
Table 1 Gauss quadrature
NG
Integration Points
(s
i
)
Weights (w
i
)
Exact for poly-
nomial of degree
1 0.0 2.0 1
2 .5773502692 1.0 3
3
.7745966692
0.0
.5555555556(5/9)
.8888888889(8/9)
5
4
.8611363116
.3399810436
.3478546451
.6521451549
7
5
.9061798459
.5384693101
0.0
.2369268851
.4786286705
.5688888889
9
- Example: evaluate
1
7 6
1
(8 7 ) 2 I x x dx
= + =
where s
i
and t
j
= Gauss points, m = number of Gauss points in s direction, n = number of
Gauss points in t direction, w
i
and w
j
= Gauss weights.
- Total number of Gauss points = mn.
- Figure 9 shows few commonly used integration formulas.
s
t
s
t
s
t
(a) 11 (b) 22 (c) 33
Figure 9 Gauss integration points in two-dimensional parent elements
13
- The element stiffness matrix can be evaluated using 22 Gauss integration formulas:
1 1
2 2
1 1
1 1
[ ] [ ][ ][ ] [ ( , )][ ][ ( , )] ( , )
T T
i j i j i j i j
i j
h dsdt h w w s t s t s t
= =
=
k B C B J B C B J
Example: Interpolation using Quadrilateral Element
For a rectangular element shown in the figure, displacements at four nodes are given
by {u
1
, v
1
, u
2
, v
2
, u
3
, v
3
, u
4
, v
4
} = {0.0, 0.0, 1.0, 0.0, 2.0, 1.0, 0.0, 2.0}. Calculate
displacement and strain at point (s, t) = (1/3, 0).
x
y
1 (0,0) 2 (3,0)
3 (3,2) 4 (0,2)
s
t
1 (-1,-1) 2 (1,-1)
3 (1,1)
4 (-1,1)
First, the shape functions are given in the parent element as
1 2
3 4
1 1
(1 )(1 ), (1 )(1 )
4 4
1 1
(1 )(1 ), (1 )(1 )
4 4
N s t N s t
N s t N s t
= = +
= + + = +
Especially when (s, t) = (1/3, 0),
1 2 3 4
1 1 1 1
, , ,
6 3 3 6
N N N N = = = =
Interpolation of geometry [find physical coordinates corresponding to (s, t) = (1/3, 0)]:
4
1
4
1
1 1 1 1
0 3 3 0 2
6 3 3 6
1 1 1 1
0 0 2 2 1
6 3 3 6
I I
I
I I
I
x N x
y N y
=
=
= = + + + =
= = + + + =
Displacement interpolation (Same shape function for geometry interpolation: Iso-
parametric)
14
4
1
4
1
1 1 1 1
0 1 2 0 1
6 3 3 6
1 1 1 1 2
0 0 1 2
6 3 3 6 3
I I
I
I I
I
u N u
v N v
=
=
= = + + + =
= = + + + =
Derivatives of shape functions (Note shape functions depend on s and t):
1 1
2 2
3 3
4 4
1 1 1 1
(1 ) (1 )
4 4 4 6
1 1 1 1
(1 ) (1 )
4 4 4 3
1 1 1 1
(1 ) (1 )
4 4 4 3
1 1 1 1
(1 ) (1 )
4 4 4 6
N N
t s
s t
N N
t s
s t
N N
t s
s t
N N
t s
s t
= = = =
= = = + =
= + = = + =
= + = = =
Jacobian matrix:
1 1 1 1 3
0 3 3 0
4 4 4 4 2
1 1 1 1
0 0 2 2 0
4 4 4 4
1 1 1 1
0 3 3 0 0
6 3 3 6
1 1 1 1
0 0 2 2 1
6 3 3 6
x
s
y
s
x
t
y
t
= + + =
= + + =
= + + =
= + + =
1
2 3
0 0
[ ] , [ ] 3 2
0 1 0 1
x y
s s
x y
t t
(
( (
(
( (
= = = (
( (
(
(
J J
Derivative of shape functions:
[ ]
I I I
I
I
I
I I I I
N N x N y
N
N
s x s y s
x
s
N
N N x N y N
y
t x t y t t
= +
=
` `
= +
)
)
J
1
2
2
0
3
[ ] 3
0 1
I
I I I
I
I I I
N
N N N
x
s s s
N
N N N
y
t t t
(
(
= = =
` ` ` `
(
) ) )
)
J
15
Strain:
4 4
1 1
4 4
1 1
2
3
2 1 1 1 1 1
( 0 1 2 0)
3 4 4 4 4 2
1 1 1 1 2
0 0 1 2
6 3 3 6 3
?
I I
xx I I
I I
I I
yy I I
I I
xy
u N N
u u
x x s
v N N
v v
y y t
= =
= =
= = =
= + + =
= = =
= + + =
=