ME6603-FINITE ELEMENT ANALYSIS

ME6603 –FINITE ELEMENT ANALYSIS

CONTENTS

UNIT – I

Fundamental Concepts

Syllabus
General Methods of the Finite Element Analysis
General Steps of the Finite Element Analysis
Boundary Conditions
Consideration during Discretization process
Rayleigh – Ritz Method (Variational
Approach) Problems (I set)
Weighted Residual method
Problems (II set)
Matrix Algebra
Matrix Operation
Gaussian Elimination Method
Problems (III set)
Advantages of Finite Element Method
Disadvantages of Finite Element Method
Applications of Finite Element Analysis

UNIT – II
One Dimension Problems
Syllabus
One Dimensional elements
Bar, Beam and Truss
Stress, Strain and Displacement
Types of Loading
Finite Element Modeling
Co – Ordinates
PEC / MECH / FEA / NOL Page 1
ME6603-FINITE ELEMENT ANALYSIS

Natural Co – Ord in at e ( ε)
Shape function
Polynomial Shape function
Stiffness Matrix [K]
Properties of Stiffness Matrix
Equation of Stiffness Matrix for One dimensional bar element
Finite Element Equation for One dimensional bar element
The Load (or) Force Vector {F}
Problem (I set)
Trusses
Stiffness Matrix [K] for a truss element
Finite Element Equation for Two noded Truss element
Problem (II set)
The Galerkin Approach
Types of beam
Types of Transverse Load
Problem (III set)

UNIT – III
Two Dimension Problems – Scalar variable Problems

Syllabus
Two dimensional elements
Plane Stress and Plane Strain
Finite Element Modeling
Constant Strain Triangular (CST) Element
Shape function for the CST element
Displacement function for the CST element
Strain – Displacement matrix [B] for CST element
Stress – Strain relationship matrix (or) Constitutive matrix [D] for two
dimensional element
Stress – Strain relationship matrix for two dimensional plane stress problems

PEC / MECH / FEA / NOL Page 2

ME6603-FINITE ELEMENT ANALYSIS

Stress – Strain relationship matrix for two dimensional plane strain proble
ms Stiffness matrix equation for two dimensional element (CST element)
Temperature Effects
Galerkin Approach
Linear Strain Triangular (LST) element
Problem (I set)
Scalar variable problems
Equation of Temperature function (T) for one dimensional heat conduction
Equation of Shape functions (N1 & N2) for one dimensional heat conduction
Equation of Stiffness Matrix (K) for one dimensional heat conduction
Finite Element Equations for one dimensional heat conduction
Finite element Equation for Torsional Bar element
Problem (II set)

UNIT – IV
AXISYMMETRIC CONTINUUM

Syllabus
Elasticity Equations
Axisymmetric Elements
Axisymmetric Formulation
Equation of shape function for Axisymmetric element
Equation of Strain – Displacement Matrix [B] for Axisymmetric
element Equation of Stress – Strain Matrix [D] for Axisymmetric
element Equation of Stiffness Matrix [K] for Axisymmetric element
Temperature Effects
Problem (I set)

PEC / MECH / FEA / NOL Page 3

ME6603-FINITE ELEMENT ANALYSIS

UNIT – V
ISOPARAMETRIC ELEMENTS FOR TWO DIMENSIONAL CONTINUUM

Syllabus
Isoparametric element
Superparametric element
Subparametric element
Equation of Shape function for 4 noded rectangular parent element
Equation of Stiffness Matrix for 4 noded isoparametric quadrilateral element
Equation of element force vector
Numerical Integration (Gaussian Quadrature)
Problem (I set)

PEC / MECH / FEA / NOL Page 4

ME6603-FINITE ELEMENT ANALYSIS

ME 1401 INTRODUCTION OF FINITE ELEMENT ANALYSIS

Unit – I Fundamental Concepts

Syllabus
Historical background – Matrix approach – Application to the continuum –
Discretisation – Matrix algebra – Gaussian elimination – Governing equations for
continuum – Classical Techniques in FEM – Weighted residual method – Ritz
method

 General Methods of the Finite Element Analysis

1. Force Method – Internal forces are considered as the unknowns of the problem.
2. Displacement or stiffness method – Displacements of the nodes are considered
as the unknowns of the problem.

 General Steps of the Finite Element Analysis

Discretization of structure > Numbering of Nodes and Elements > Selection of
Displacement function or interpolation function > Define the material behavior by
using Strain – Displacement and Stress – Strain relationships > Derivation of
element stiffness matrix and equations > Assemble the element equations to
obtain the global or total equations > Applying boundary conditions > Solution
for the unknown displacements > computation of the element strains and stresses
from the nodal displacements > Interpret the results (post processing).

 Boundary Conditions
It can be either on displacements or on stresses. The boundary conditions on
displacements to prevail at certain points on the boundary of the body, whereas
the boundary conditions on stresses require that the stresses induced must be in
equilibrium with the external forces applied at certain points on the boundary of
the body.

 Consideration During Discretization process

Types of element > Size of element > Location of node > Number of elements.

PEC / MECH / FEA / NOL Page 5

ME6603-FINITE ELEMENT ANALYSIS

 Rayleigh – Ritz Method (Variational Approach)

It is useful for solving complex structural problems. This method is possible only
if a suitable functional is available. Otherwise, Galerkin’s method of weighted
residual is used.

 Problems (I set)
1. A simply supported beam subjected to uniformly distributed load over entire
span. Determine the bending moment and deflection at midspan by using
Rayleigh – Ritz method and compare with exact solutions.
2. A bar of uniform cross section is clamed at one end and left free at another end
and it is subjected to a uniform axial load P. Calculate the displacement and stress
in a bar by using two terms polynomial and three terms polynomial. Compare
with exact solutions.

 Weighted Residual method

It is a powerful approximate procedure applicable to several problems. For non –
structural problems, the method of weighted residuals becomes very useful. It has
many types. The popular four methods are,
1. Point collocation method,
Residuals are set to zero at n different locations X i, and the weighting function wi
is denoted as (x - xi).

  (x  xi) R (x; a1, a2, a3… an) dx = 0

2. Subdomain collocation method,

 1 forxinD1
w1 = 
0 forxnotinD1
3. Least square method,
 [R (x; a1, a2, a3… an)]2 dx = minimum.

4. Galerkin’s method.
wi = Ni (x)
 Ni (x) [R (x; a1, a2, a3… an)]2 dx = 0, i = 1, 2, 3, …n.
PEC / MECH / FEA / NOL Page 6
ME6603-FINITE ELEMENT ANALYSIS

 Problems (II set)

1. The following differential equation is available for a physical phenomenon.
2
d y
+ 50 = 0, 0  x 10. Trial function is y = a1x (10-x). Boundary conditions
dx 2
are y (0) = 0 and y (10) = 0. Find the value of the parameter a1 by the following
methods, (1) Point collocation method, (2) Subdomain collocation method, (3)
Least square method and (4) Galerkin’s method.

2. The differential equation of a physical phenomenon is given by

d2y
2
- 10x2 = 5. Obtain two term Galerkin solution b using the trial functions:
dx
N1(x) = x(x-1); N2(x) = x2(x-1); 0  x 1. Boundary conditions are y (0) = 0 and
y (1) = 0.

 Matrix Algebra
Equal matrix: Two matrixes are having same order and corresponding elements
are equal.
Diagonal matrix: Square matrix in which all the elements other than the diagonal
are zero.
Scalar matrix: Square matrix in which all the elements are equal.
Unit matrix: All diagonal elements are unity and other elements are zero.

 Matrix Operation
Scalar multiplication, Addition and Subtraction of matrices, Multiplication of
matrices, Transpose of a matrix, Determinant of a matrix, inverse of a matrix,
Cofactor or Adjoint method to determine the inverse of a matrix, Row reduction
method (Gauss Jordan method) to determine the inverse of a matrix, Matrix
differentiation and matrix integration.
 Gaussian Elimination Method
It is most commonly used for solving simultaneous linear equations. It is easily
adapted to the computer for solving such equations.

PEC / MECH / FEA / NOL Page 7

ME6603-FINITE ELEMENT ANALYSIS

 Problems (III set)

1. 3x+y-z = 3, 2x-8y+z = -5, x-2y+9z = 8. Solve by using Gauss – Elimination
method.
2. 2a+4b+2c = 15, 2a+b+2c = -5, 4a+b-2c = 0. Solve the equations by using
Gauss – Elimination method.

 Advantages of Finite Element Method

1. FEM can handle irregular geometry in a convenient manner.
2. Handles general load conditions without difficulty
3. Non – homogeneous materials can be handled easily.
4. Higher order elements may be implemented.

 Disadvantages of Finite Element Method

1. It requires a digital computer and fairly extensive
2. It requires longer execution time compared with FEM.
3. Output result will vary considerably.

 Applications of Finite Element Analysis

Structural Problems:
1. Stress analysis including truss and frame analysis
2. Stress concentration problems typically associated with holes, fillets or
other changes in geometry in a body.
3. Buckling Analysis: Example: Connecting rod subjected to axial
compression.
4. Vibration Analysis: Example: A beam subjected to different types of
loading.
Non - Structural Problems:
1. Heat Transfer analysis:
Example: Steady state thermal analysis on composite cylinder.
2. Fluid flow analysis:
Example: Fluid flow through pipes.

PEC / MECH / FEA / NOL Page 8

ME6603-FINITE ELEMENT ANALYSIS

Unit – II One Dimension Problems

Syllabus
Finite element modeling – Coordinates and shape functions- Potential energy
approach – Galarkin approach – Assembly of stiffness matrix and load vector – Finite
element equations – Quadratic shape functions – Applications to plane trusses

 One Dimensional elements

Bar and beam elements are considered as One Dimensional elements. These
elements are often used to model trusses and frame structures.

 Bar, Beam and Truss

Bar is a member which resists only axial loads. A beam can resist axial, lateral
and twisting loads. A truss is an assemblage of bars with pin joints and a frame is
an assemblage of beam elements.

 Stress, Strain and Displacement

Stress is denoted in the form of vector by the variable x as σx, Strain is denoted in
the form of vector by the variable x as e x, Displacement is denoted in the form of
vector by the variable x as ux.

 Types of Loading
(1) Body force (f)
It is a distributed force acting on every elemental volume of the body. Unit
is Force / Unit volume. Ex: Self weight due to gravity.
(2) Traction (T)
It is a distributed force acting on the surface of the body. Unit is
Force / Unit area. But for one dimensional problem, unit is Force / Unit length.
Ex: Frictional resistance, viscous drag and Surface shear.
(3) Point load (P)
It is a force acting at a particular point which causes displacement.

PEC / MECH / FEA / NOL Page 9

ME6603-FINITE ELEMENT ANALYSIS

 Finite Element Modeling

It has two processes.
(1) Discretization of structure
(2) Numbering of nodes.

 CO – ORDINATES
(A) Global co – ordinates, (B) Local co – ordinates and (C) Natural co –
ordinates.

 Natural Co – Ordinate (ε)

pc
ε =  x2  x1 
 
 2 

Integration of polynomial terms in natural co – ordinates for two dimensional
elements can be performed by using the formula,


 !  ! !
 L1  L2  L3  dA 
 
X 2A
     !

PEC / MECH / FEA / NOL Page 10

ME6603-FINITE ELEMENT ANALYSIS

 Shape function
N1N2N3 are usually denoted as shape function. In one dimensional problem,
the displacement
u =  Ni ui =N1 u1
For two noded bar element, the displacement at any point within the element
is given by,
u =  Ni ui =N1 u1 + N2 u2
For three noded triangular element, the displacement at any point within
the element is given by,=  Ni ui =N1 u1 + N2 u2 + N3 u3

v =  Ni vi =N1 v1 + N2 v2 + N3 v3

Shape function need to satisfy the following

(a) First derivatives should be finite within an element; (b) Displacement should
be continuous across the element boundary.

 Polynomial Shape function

Polynomials are used as shape function due to the following reasons,
(1) Differentiation and integration of polynomials are quite easy.
(2) It is easy to formulate and computerize the finite element equations.
(3) The accuracy of the results can be improved by increasing the order of the
polynomial.

 Stiffness Matrix [K]

 B DBdv
T
Stiffness Matrix [K] =
V
 Properties of Stiffness Matrix
1. It is a symmetric matrix, 2. The sum of elements in any column must be equal
to zero, 3. It is an unstable element. So the determinant is equal to zero.

 Equation of Stiffness Matrix for One dimensional bar element

AE  1  1
[K] =  
PEC / MECH / FEA / NOL Page 11
ME6603-FINITE ELEMENT ANALYSIS

l  1 1 

 Finite Element Equation for One dimensional bar element

 F1  AE  1 1u1 
  1 1  u 
 2 
F l   2 

PEC / MECH / FEA / NOL Page 12

ME6603-FINITE ELEMENT ANALYSIS



 The Load (or) Force Vector {F}


Al 1
F e  
2 1

 Problem (I set)
1. A two noded truss element is shown in figure. The nodal displacements are
u1 = 5 mm and u2 = 8 mm. Calculate the displacement at x = ¼, 1/3 and ½.

 Trusses
It is made up of several bars, riveted or welded together. The following
assumptions are made while finding the forces in a truss,
(a) All members are pin joints, (b) The truss is loaded only at the joints, (c) The
self – weight of the members is neglected unless stated.

 Stiffness Matrix [K] for a truss element

 l2 lm l2  lm 
 2 


2
 lm  m 
K   e e  2
A E lm m
le l
  lm l2 lm 
 
 lm  m 2 lm m 2

PEC / MECH / FEA / NOL Page 13

ME6603-FINITE ELEMENT ANALYSIS

 Finite Element Equation for Two noded Truss element

 F1   l2 lm  l2  lm  u1 
 F    
 2  Ae E e  lm m2  lm  m 2  u 2 
    l 2  
 F3  le  lm l2 lm  u 3 
 
F4  lm  m2 lm m 2 u4

 Problem (II set)

1. Consider a three bar truss as shown in figure. It is given that E = 2 x 10 5
N/mm2. Calculate (a) Nodal displacement, (b) Stress in each member and
(c) Reactions at the support. Take Area of element 1 = 2000 mm2, Area of
element 2 = 2500 mm2, Area of element 3 = 2500 mm2.

 
 The Galerkin Approach  
 1  1
K    AE 
  
Stiffness Matrix
l   1 1 

 Types of beam
1. Cantilever beam, 2. Simply Supported beam, 3. Over hanging beam, 4. Fixed
beam and 5. Continuous beam.

PEC / MECH / FEA / NOL Page 14

ME6603-FINITE ELEMENT ANALYSIS

 Types of Transverse Load

1. Point or Concentrated Load, 2. Uniformly Distributed Load and 3. Uniformly
Varying Load.

 Problem (III set)

1. A fixed beam of length 2L m carries a uniformly distributed load of w (N/m)
which runs over a length of L m from the fixed end. Calculate the rotation at Point
B.

PEC / MECH / FEA / NOL Page 15

ME6603-FINITE ELEMENT ANALYSIS

Unit – III Two Dimension Problems – Scalar variable Problems

Syllabus
Finite element modeling – CST & LST elements – Elements equations – Load vectors and
boundary conditions – Assembly – Applications to scalar variable problems such as torsion,
heat transfer.

 Two dimensional elements

Two dimensional elements are defined by three or more nodes in a two
dimensional plane (i.e., x, y plane). The basic element useful for two dimensional
analysis is the triangular element.

 Plane Stress and Plane Strain

The 2d element is extremely important for the Plane Stress analysis and Plane
Strain analysis.
Plane Stress Analysis:
It is defined to be a state of stress in which the normal stress () and shear
stress () directed perpendicular to the plane are assumed to be zero.
Plane Strain Analysis:
It is defined to be a state of strain in which the normal to the xy plane and
the shear strain are assumed to be zero.

 Finite Element Modeling

It consists of 1. Discretization of structure and 2. Numbering of nodes.
1. Discretization:
The art of subdividing a structure into a convenient number of smaller
components is known as discretization.
PEC / MECH / FEA / NOL Page 16
ME6603-FINITE ELEMENT ANALYSIS

2. Numbering of nodes:
In one dimensional problem, each node is allowed to move only in  x
direction. But in two dimensional problem, each node is permitted to move in the two
directions i.e., x and y.

The element connectivity table for the above domain is explained as table.

Element (e) Nodes

(1) 123
(2) 234
(3) 435
(4) 536
(5) 637
(6) 738
(7) 839
(8) 931

 Constant Strain Triangular (CST) Element

A three noded triangular element is known as constant strain triangular (CST)
element. It has six unknown displacement degrees of freedom (u 1v1, u2v2, u3v3).

PEC / MECH / FEA / NOL Page 17

ME6603-FINITE ELEMENT ANALYSIS

 Shape function for the CST element

Shape function N1 = (p1 + q1x + r1y) / 2A
Shape function N2 = (p2 + q2x + r2y) / 2A
Shape function N3 = (p3 + q3x + r3y) / 2A

 Displacement function for the CST element

 u1 
   v1 
   
 0  u2
u(x, y)  N1 0 N2 0 N3
Displacement function u =  X
    
v(x, y)   0 N1 0 N2 0 N 3 v2 
u3
 
 v3 



 Strain – Displacement matrix [B] for CST element

 q1 0 q2 0 q3 0 
1 
Strain – Displacement matrix [B] = 0 r1 0 r2 0 r3 
2A  
 r1 q1 r2 q2 r3 q3
Where, q1 = y2 – y3 r1 = x3 – x2
q2 = y3 – y1 r2 = x1 – x3
q3 = y1 – y2 r3 = x2 – x1
 Stress – Strain relationship matrix (or) Constitutive matrix [D] for two
dimensional element

 
 1 v v v 0 0 0 
 
 v 1 v v 0 0 0 
 
 v v 1 v 0 0 0 
E  
[D] = 
1  v  1  2v   0 1 2v 
0 0 0 0 
 2 
 0 1 2v 
0 0 0 0 
 2 
 1  2v 
 0 0 0 0 0 
PEC / MECH / FEA / NOL Page 18
ME6603-FINITE ELEMENT ANALYSIS

 Finite Element Equation for Two noded Truss element
 Stress – Strain relationship matrix for two dimensional plane stress problems
The normal stress z and shear stresses xz, yz are zero.

 
1 v 0
E  

 v 1 0 
[D] = 1  v 2 
1  v 
 0 0 
 2 



 Stress – Strain relationship matrix for two dimensional plane strain
problems
Normal strain e z and shear strains exz, eyz are zero.

 
(1 v) v 0 
E  v 
(1 v) 0 
[D] =
1 v  1 2v   1 2v 
 0 0
 2 



 Stiffness matrix equation for two dimensional element (CST element)
Stiffness matrix [k] = [B]T [D] [B] A t

 q1 0 q2 0 q3 0 
1 
0 r1 0 r2 0 r3 
[B] =
2A 
 r1 q1 r2 q2 r3 q3
For plane stress problems,

 
1 v 0 
E 
v 1 0 
[D] = 1  v 2 
1  v 
0 0 
  2 

For plane strain problems,

PEC / MECH / FEA / NOL Page 19

ME6603-FINITE ELEMENT ANALYSIS

  Finite Element Equation for Two noded Truss element
 
(1 v) v 0 
E  v 
(1 v) 0 
[D] =
1 v  1 2v   1 2v 
 0 0
 2 



 Temperature Effects
Distribution of the change in temperature (ΔT) is known as strain. Due to the
change in temperature can be considered as an initial strain e0.
σ = D (Bu - e0)

 Galerkin Approach
Stiffness matrix [K] e = [B]T [D][B] A t.
Force Vector {F}e = [K]e {u}

 Linear Strain Triangular (LST) element

A six noded triangular element is known as Linear Strain Triangular (LST)
element. It has twelve unknown displacement degrees of freedom. The displacement
functions of the element are quadratic instead of linear as in the CST.

PEC / MECH / FEA / NOL Page 20

ME6603-FINITE ELEMENT ANALYSIS

 Finite Element Equation for Two noded Truss element

 Problem (I set)
1. Determine the shape functions N1, N2 and N3 at the interior point P for the
triangular element for the given figure.

The two dimensional propped beam shown in figure. It is divided into two CST
elements. Determine the nodal displacement and element stresses using plane stress
conditions. Body force is neglected in comparison with the external forces.
Take, Thickness (t) = 10mm,
Young’s modulus (E) = 2x105 N/mm2,
Poisson’s ratio (v) = 0.25.

3. A thin plate is subjected to surface traction as in figure. Calculate the global

stiffness matrix.

PEC / MECH / FEA / NOL Page 21

ME6603-FINITE ELEMENT ANALYSIS

 Finite Element Equation for Two noded Truss element

 Scalar variable problems

In structural problems, displacement at each nodal point is obtained. By using
these displacement solutions, stresses and strains are calculated for each element. In
structural problems, the unknowns (displacements) are represented by the
components of vector field. For example, in a two dimensional plate, the unknown
quantity is the vector field u(x, y), where u is a (2x1) displacement vector.

 Equation of Temperature function (T) for one dimensional heat conduction

Temperature (T) = N1T1 + N2T2
 Equation of Shape functions (N1 & N2) for one dimensional heat conduction
lx
N1 =
l
x
N2 =
l

 Equation of Stiffness Matrix (K) for one dimensional heat conduction

 1  1
K c   Ak  1 1 
l  

PEC / MECH / FEA / NOL Page 22

ME6603-FINITE ELEMENT ANALYSIS

  Finite Element Equation for Two noded Truss element


 Finite Element Equations for one dimensional heat conduction
Case (i): One dimensional heat conduction with free end convection

 Ak  1  1 0 0  T1   0

 
   hA     hT A 
 l  1 1   0 1   T2  1
Case (ii): One dimensional element with conduction, convection and internal heat
generation.

 Ak  11 hPl 2 1  T1   QAl PhT l 1

 
   1 2 T   
 l 1 1  6     2  2 1



 Finite element Equation for Torsional Bar element

 M 1x  GJ  1  11x 
    1 1   
 M 2 x  l   2 x 

  1
GJ  1
Where, Stiffness matrix [K] =  

PEC / MECH / FEA / NOL Page 23

ME6603-FINITE ELEMENT ANALYSIS
l  1 1 



 Problem (II set)
1. An Aluminium alloy fin of 7 mm thick and 50 mm long protrudes from a wall,
which is maintained at 120°C. The ambient air temperature is 22°C. The heat
transfer coefficient and thermal conductivity of fin material are 140 W/m2K and 55
W/mK respectively. Determine the temperature distribution of fin.

2. Calculate the temperature distribution in a one dimension fin with physical

properties given in figure. The fin is rectangular in shape and is 120 mm long, 40mm
wide and 10mm thick. Assume that convection heat loss occurs from the end of the
fin. Use two elements. Take k = 0.3W/mm°C, h = 1 x 10-3 W/ mm2 °C, T=20°C.

PEC / MECH / FEA / NOL Page 24

ME6603-FINITE ELEMENT ANALYSIS

Unit – IV AXISYMMETRIC CONTINUUM

Syllabus
Axisymmetric formulation – Element stiffness matrix and force vector – Galarkin approach
– Body forces and temperature effects – Stress calculations – Boundary conditions –
Applications to cylinders under internal or external pressures – Rotating discs

 Elasticity Equations
Elasticity equations are used for solving structural mechanics problems. These
equations must be satisfied if an exact solution to a structural mechanics problem
is to be obtained. The types of elasticity equations are
1. Strian – Displacement relationship equations

u v u v  u w
ex  ; ey  ;  xy   ;  xz   ;
x y y x z x
v w
 yz   .
z y
ex – Strain in X direction, ey – Strain in Y direction.
 xy - Shear Strain in XY plane,  xz - Shear Strain in XZ plane,

 yz - Shear Strain in YZ plane

2. Sterss – Strain relationship equation

1  v  v v 0 0 0 
  v 
 x 
 1  v  v 0 0 0   e x 
   
 y   v v 1  v  0 0 0   e x 
 z    
1 2v   

E
  0 0 0 0 0   e x 
 
 xy  1  v 1  2v  
 2   xy 


 0 1 2v 

 yz  0 0 0 0  yz 
    2   
 zx  1 2v   zx

 0 0 0 0 0 
 2 

PEC / MECH / FEA / NOL Page 25

ME6603-FINITE ELEMENT ANALYSIS


σ – Stress, τ – Shear Stress, E – Young’s Modulus, v – Poisson’s Ratio,
e – Strain,  - Shear Strain.
3. Equilibrium equations

x  xy xz y yz xy

   Bx  0 ;    By  0
x y z y z x
z xz  yz
   Bz  0
z x y
σ – Stress, τ – Shear Stress, Bx - Body force at X direction,

B y - Body force at Y direction, Bz - Body force at Z direction.

4. Compatibility equations
There are six independent compatibility equations, one of which is

e
2
e2y  2xy
x
 
y 2 x 2 xy .

The other five equations are similarly second order relations.

 Axisymmetric Elements
Most of the three dimensional problems are symmetry about an axis of rotation.
Those types of problems are solved by a special two dimensional element called
as axisymmetric element.

PEC / MECH / FEA / NOL Page 26

ME6603-FINITE ELEMENT ANALYSIS

 Axisymmetric Formulation
The displacement vector u is given by

 u 
u r, z    
 w 



The stress σ is given by

 r 
 
  
Stress,     
 z 





 rz 

The strain e is given by


 er 
e 
 
Strain, e    
 ez 
 rz
 Equation of shape function for Axisymmetric element
Shape function,

1 1r  1 z  2  2 r 2 z  3  3 r 3 z

N1  ; N2  ; N3 
2A 2A 2A
α1 = r2z3 – r3z2; α2 =r3z1 – r1z3; α3 = r1z2 – r2z1
ᵝ1 = z2-z3; ᵝ2 = z3-z1; ᵝ3 = z1-z2
ᵞ1 = r3-r2; ᵞ2 = r1-r3; ᵞ3 = r2-r1
2A = (r2z3 – r3z2)-r1(r3z1 – r1z3)+z1(r1z2 – r2z1)

PEC / MECH / FEA / NOL Page 27

ME6603-FINITE ELEMENT ANALYSIS

 Equation of Strain – Displacement Matrix [B] for Axisymmetric element

 u1 
 1 0 2 0 3 0   w 
 2  3z   
1
 z 2z
1  1  1  1  2  3
 3  0  u 2 
B    r r
0
r r
0
r r  
2A
 0 1 0 2 0  3  w2 
 1 1 2 2 3  3  u 3 
 
w3

r1r 2 r 3
r
3

 Equation of Stress – Strain Matrix [D] for Axisymmetric element

1 v v v 0 
 v 1 v v 0 
 
D   E
0 
1 v 1 2v   v v 1 v
12v 
 0 0 0
 2
 Equation of Stiffness Matrix [K] for Axisymmetric element


K   2rABT DB
r1r 2 r 3
r ; A = (½) bxh
3
 Temperature Effects
The thermal force vector is given by

 f t  2rABDet

PEC / MECH / FEA / NOL Page 28

ME6603-FINITE ELEMENT ANALYSIS

 F1u 
  F w 
  1 
  F u 
 f t   2 
 F2 w
 F3u 
 


 F
 3  
w

 Problem (I set)

1. For the given element, determine the stiffness matrix. Take E=200GPa and
υ= 0.25.

2. For the figure, determine the element stresses. Take E=2.1x105N/mm2 and
υ= 0.25. The co – ordinates are in mm. The nodal displacements are u 1=0.05mm,
w1=0.03mm, u2=0.02mm, w2=0.02mm, u3=0.0mm, w3=0.0mm.

3. A long hollow cylinder of inside diameter 100mm and outside diameter

140mm is subjected to an internal pressure of 4N/mm2. By using two
elements on the 15mm length, calculate the displacements at the inner radius.

PEC / MECH / FEA / NOL Page 29

ME6603-FINITE ELEMENT ANALYSIS

UNIT – V ISOPARAMETRIC ELEMENTS FOR TWO DIMENSIONAL

CONTINUUM

Syllabus
The four node quadrilateral – Shape functions – Element stiffness matrix and force
vector – Numerical integration - Stiffness integration – Stress calculations – Four node
quadrilateral for axisymmetric problems.

 Isoparametric element
Generally it is very difficult to represent the curved boundaries by straight edge
elements. A large number of elements may be used to obtain reasonable
resemblance between original body and the assemblage. In order to overcome this
drawback, isoparametric elements are used.

If the number of nodes used for defining the geometry is same as number of nodes
used defining the displacements, then it is known as isoparametric element.
 Superparametric element
If the number of nodes used for defining the geometry is more than number of
nodes used for defining the displacements, then it is known as superparametric
element.

PEC / MECH / FEA / NOL Page 30

ME6603-FINITE ELEMENT ANALYSIS

 Finite Element Equation for Two noded Truss element

 Subparametric element
If the number of nodes used for defining the geometry is less than number of
nodes used for defining the displacements, then it is known as subparametric
element.

PEC / MECH / FEA / NOL Page 31

ME6603-FINITE ELEMENT ANALYSIS

 Finite Element Equation for Two noded Truss element

 Equation of Shape function for 4 noded rectangular parent element

 x1 
 y 
 2 
 x1 
 
x  N1 0 N2 0 N3 0 N4 0   y 2 
u 
    
 y  0 N1 0 N2 0 N3 0 N 4   x3 
 y3 
 
 x4 
y 
 4 


N1=1/4(1-Ɛ) (1-ɳ); N2=1/4(1+Ɛ) (1-ɳ); N3=1/4(1+Ɛ) (1+ɳ); N4=1/4(1-Ɛ) (1+ɳ).

 Equation of Stiffness Matrix for 4 noded isoparametric quadrilateral element

1 1 T

K   t   B DB J 

1 1

J J 12 
J    11 ;
 J 21 J 22 

J11   1  x1  1  x2  (1  )x3  (1  )x4  ;

1
4

J 12 1  1   y1  1   y 2  (1  ) y3  (1  ) y4  ;

4
J 21  1  1  x  1  x  (1  )x  (1  )x 
1 2 3 4 ;
4
J 22  1
 1   y1  1   y2  (1  ) y3  (1  ) y4 ;
4

PEC / MECH / FEA / NOL Page 32

ME6603-FINITE ELEMENT ANALYSIS

 Finite Element Equation for Two noded Truss element

   
   N1 N 2 N 3 N 4 
    0 0 0 0 
    
 N N N 3 N 4 
 J 22  J 12 0 0   1 0 2
0 0 0 
   
 J 21 J 11  
1  
B    0 0
 N1 N 2 N 3 N 4 
J
 J 21 J 11 J 22  J 12   0 0 0 0 
     
 0 N1 N 2 N 3 N 4 
0 0 0
   


 
1 v v 0 
 E 
[D]   v 1 0 
(1 v ) 
2
1 v , for plane stress conditions;
 0 0 
 2 
 
 
1 v v 0 
 E
[D]   v 1 v 0 
(1 v)(1 2v)  1 2v , for plane strain conditions.
 0 0 
 2 

 Equation of element force vector

 
F e  [N ]  ;
F x
T


 Fy 
N – Shape function, Fx – load or force along x direction,
Fy – load or force along y direction.

 Numerical Integration (Gaussian Quadrature)

The Gauss quadrature is one of the numerical integration methods to calculate the
definite integrals. In FEA, this Gauss quadrature method is mostly preferred. In
this method the numerical integration is achieved by the following expression,
1 n

 f (x)dx   w f (x )
1 i 1
i i

Table gives gauss points for integration from -1 to 1.

PEC / MECH / FEA / NOL Page 33

ME6603-FINITE ELEMENT ANALYSIS

Number of Location Corresponding Weights

Points xi wi
n

1 x1 = 0.000 2.000

2 1
x1, x2 =   0.577350269189
3 1.000

5
x1, x3  
3
 0.774596669241  0.555555
3 5 9

x2=0.000 8
 0.888888
9

4 x1, x4 =  0.8611363116 0.3478548451

x2, x3 =  0.3399810436 0.6521451549

 Problem (I set)

x
1

1. Evaluate I  cos
1
2
dx , by applying 3 point Gaussian quadrature and

compare with exact solution.

 x 1 
1

2. Evaluate I   3e  x 
2

 dx , using one point and two point


1 
x 2 
Gaussian quadrature. Compare with exact solution.
3. For the isoparametric quadrilateral element shown in figure, determine the
local co –ordinates of the point P which has Cartesian co-ordinates (7, 4).

PEC / MECH / FEA / NOL Page 34

ME6603-FINITE ELEMENT ANALYSIS

4. A four noded rectangular element is in figure. Determine (i) Jacobian

matrix, (ii) Strain – Displacement matrix and (iii) Element Stresses. Take
E=2x105N/mm2,υ= 0.25, u=[0,0,0.003,0.004,0.006, 0.004,0,0]T, Ɛ= 0, ɳ=0.
Assume plane stress condition.

PEC / MECH / FEA / NOL Page 35