FEM
Collected by
Jyoti Anischit
MSc Engineering Geology.
III semester.
Course Content
• Introduction, Definition of Finite Element Method,
Differential Equation and Weak Form, Variational
Principal,
• Ritz‐Galerkin Method (approximate function, Galerkin
Method and Ritz Method), Finite Element Method (1‐D
• Problem): Construction of approximate function,
Element matrix, Total element matrix and simple
example,
• Finite Element Method (2‐D Problem): Construction of
approximate function, Element matrix & total element
• matrix, simple example and Gauss’s method of
elimination.
Finite Element Method – Introduction
• The Finite Element Method (FEM) is a numerical
method of solving systems of partial differential
equations (PDEs)
• It reduces a PDE system to a system of algebraic
equations that can be solved using traditional linear
algebra techniques.
• In simple terms, FEM is a method for dividing up a
very complicated problem into small elements that
can be solved in relation to each other.
“…The laws of nature are written in the language of
mathematics. These often take the form of ordinary or
partial differential equations.
The electronic digital computer is an amazingly fast
calculating tool; but it can handle only arithmetic.
The differential equation governing the physical
phenomena… have therefore to be reduced to a system of
simultaneous equation, before the computer can solve
then by a series of arithmetical operations.
Finite Element techniques do precisely this.”
– G. S. Ramaswamy
in “Design and Construction of
Concrete Shell Roofs”
Where FEM started from…
5
Where FEM started from…
1) The first paper in FEM:
• M. J. Turner, R. W. Clough, H. C. Martin & L. J. Topp,
“Stiffness and Deflection Analysis of Complex
Structures”, J. Aeronautical Science 23 (9), pp. 805-823,
Sept. 1956.
2) The coining of the name ‘Finite Element Method’:
• R. W. Clough, “The Finite Element Method in Plane Stress
Analysis”, Proc. 2nd ASCE Conf. On Electronic Computation,
Pittsburg, Sept. 1960.
6
Where FEM started from…
3) Introduction of Isoparametric elements that made
FEM so very versatile:
• B. M. Irons & O. C. Zienkiewicz, “The Isoparametric Finite
Element System – A New Concept in Finite Element
Analysis”, Proc. Conf. Recent Advances in Stress Analysis,
Royal Aeronautical Society, London, 1968.
4) SAP-IV -- the first FEM Package :
• K. J. Bathe, E. L. Wilson & F. E. Peterson, “SAP IV – A
Structural Analysis Program for Static and Dynamic
Response of Linear Systems”, Report No. 73/11, Earthquake
Engineering Research Center, June 1973.
7
F.E.M.
• In finite element method, the structure to be
analyzed is subdivided into a mesh of finite-
sized elements of simple shape, and then the
whole structure is solved with quite easiness.
Finite Sized Element
Rectangular Body Circular Plate
Finite Sized Elements
• The rectangular panel in the rectangular body
and triangular panel in the circular plate are
referred to an ‘element’.
• There’re one-, two- and three-dimensional
elements.
• The accuracy of the solution depends upon
the number of the finite elements; the more
there’re, the greater the accuracy.
Finite Element of a Bar
• If a uniaxial bar is part of a structure then it’s
usually modeled by a spring element if and
only if the bar is allowed to move freely due to
the displacement of the whole structure. (One
dimensional element)
Bar
Spring element
Types of Elements
• Here goes the examples of two- and three-
dimensional finite sized elements.
Triangle
Rectangle
Hexahedron
Node
• The points of attachment of the element to
other parts of the structure are called nodes.
• The displacement at any node due to the
deformation of structure is known as the
nodal displacement.
Node
Why F.E.M.?
Simple trusses can be solved by just using the
equilibrium equations. But for the complex
shapes and frameworks like a circular plate,
equilibrium equations can no longer be applied as
the plate is an elastic continuum not the beams or
bars as the case of normal trusses.
Hence, metal plate is divided into finite
subdivisions (elements) and each element is
treated as the beam or bar. And now stress
distribution at any part can be determined
accurately.
Why F.E.M.?
Simple Bar Analysis
• Consider a simple bar made up of uniform
material with length L and the cross-sectional
area A. The young modulus of the material is
E. L
• Since any bar is modeled as spring in FEM thus
we’ve:
F1 x1 k x2 F2
Simple Bar Analysis
• Let us suppose that the value of spring constant is
k. Now, we’ll evaluate the value of k in terms of
the properties (length, area, etc.) of the bar:
We know that:
i.e.
Also: i.e.
And i.e.
Simple Bar Analysis
• Now substituting the values of x and F is the
base equation of k, we’ll have:
But
Hence, we may write:
Simple Bar Analysis
• According to the diagram, the force at node x1
can be written in the form:
• Where x1 – x2 is actually the nodal displacement
between two nodes. Further:
• Similarly:
Simple Bar Analysis
• Now further simplification gives:
• These two equations for F1 and F2 can also be
written as, in Matrix form:
• Or:
Simple Bar Analysis
• Here Ke is known as the Stiffness Matrix. So a
uniform material framework of bars, the value of
the stiffness matrix would remain the same for all
the elements of bars in the FEM structure.
Further Extension
k1 k2
F1 F2 F3
x1 x2 x3
• Similarly for two different materials bars joined
together, we may write:
;
Differential Equation
A Differential Equation is an equation containing the derivative of one or more
dependent variables with respect to one or more independent variables.
For example,
Classifiation by Type:
Ordinary Differential Equation
If a Differential Equations contains only ordinary derivatives of one or more
dependent variables with respect to a single independent variables, it is said to be an
Ordinary Differential Equation or (ODE) for short.
For Example,
Partial Differential Equation
If a Differential Equations contains partial derivatives of one or more dependent
variables of two or more independent variables, it is said to be a Partial Differential
Equation or (PDE) for short.
For Example,
Classifiation by Order:
The order of the differential equation (either ODE or PDE) is the order of the highest
derivative in the equation.
For Example,
Order = 3
Order = 2
Order = 1
General form of nth Order ODE is
= f(x,y,y1,y2,….,y(n))
where f is a real valued continuous function.
This is also referred to as Normal Form Of nth Order Derivative
So, when n=1, = f(x,y)
when n=2, = f(x,y,y1) and so on …
Importance
• FEM has become very familiar in subdivision of
continuum. It gives reliable and accurate results if
the number of elements are kept greater.
• Modern computer technology had helped this
analysis to be very easy and less time consuming.
• Large structures under loadings are now easily
solved and stresses on each and every part are
now being determined.
Galerkin method
• Galerkin suggested that the residue should be
multiplied by a weighting function that is a
part of the suggested solution then the
integration is performed over the whole
domain!!!
• Actually, it turned out to be a VERY GOOD idea
Galerkin Method
• Engineering problems: differential equations with
boundary conditions.
• Generally denoted as: D(U)=0; B(U)=0
• Our task: to find the function U which satisfies the
given differential equations and boundary conditions.
• Reality: difficult, even impossible to solve the problem
analytically
• In practical cases we often apply approximation.
• One of the approximation methods:
• Galerkin Method, invented by Russian mathematician
Boris Grigoryevich Galerkin.
1D Rod Elements
• To understand and solve 2D and 3D problems we
must understand basic of 1D problems.
• Analysis of 1D rod elements can be done using
Rayleigh-Ritz and Galerkin’s method.
• To solve FEA problems same are modified in the
Potential-Energy approach and Galerkin’s
approach
1D Rod Elements
• Loading consists of three types : body force f ,
traction force T, point load Pi
• Body force: distributed force , acting on every
elemental volume of body i.e. self weight of
body.
• Traction force: distributed force , acting on
surface of body i.e. frictional resistance,
viscous drag and surface shear
1D Rod Elements
Element -1 Element-2
1 T e
• Element strain energy U e q [k ]q
2
• Element stiffness matrix Ee Ae 1 1
[k ]
e
1 1
le
• Load vectors e Ae le f 1
f
– Element body load vector 2 1
– Element traction-force vector e Tle 1
T
2 1
Bar application
2u
EA 2 F x 0
x
n
u x ai i x
i 1
n
d 2 i x
EA ai 2
F x Rx Applying Galerkin method
i 1 dx
n
d 2 i x
EA ai j x 2
dx j x F x dx
i 1 Domain
dx Domain
In Matrix Form
d 2 i x
EA j x dx
ia j
x F x dx
2
Domain dx Domain
Solve the above system for the “generalized
coordinates” ai to get the solution for u(x)
Same conditions on the functions are
applied
• They should be at least twice differentiable!
• They should satisfy all boundary conditions!
• Let’s use the same function as in the
collocation method:
x
x Sin
2l
Substituting with the approximate solution:
n
d 2 i x
EA ai j x 2
dx j x F x dx
i 1 Domain
dx Domain
x x
2 l
EA a1 Sin Sin dx
2l 0 2l 2l
x
l
Sin fdx
0 2l
2
l 2l f 16 l 2 fl 2
EA a1 a1 0.52
2l 2 EA 3
EA
Substituting with the approximate solution:
(Int. by Parts)
n
d 2 i x
EA ai j x 2
dx j x F x dx
i 1 Domain
dx Domain
d 2 i x
j x 2
dx
Domain
dx
d i x
l
d j x d i x
j x dx
Zero! dx 0 Domain dx dx
2
l 2l f 16 l 2 fl 2
EA a1 a1 0.52
2l 2 EA 3
EA
What did we gain?
• The functions are required to be less
differentiable
• Not all boundary conditions need to be
satisfied
• The matrix became symmetric!
The Finite Element Method
2nd order DE’s in 1-D
Objectives
• Understand the basic steps of the finite
element analysis
• Apply the finite element method to second
order differential equations in 1-D
The Mathematical Model
• Solve:
d du
a cu f 0
dx dx
0 x L
• Subject to:
du
u 0 u0 , a Q0
dx x L
Step #1: Discretization
• At this step, we divide
the domain into
elements.
• The elements are
connected at nodes.
• All properties of the
domain are defined at
those nodes.
Step #2: Element Equations
• Let’s concentrate our
attention to a single
element. d du
• The same DE applies on a cu f 0
dx dx
the element level, hence,
we may follow the x1 x x2
procedure for weighted u x1 u1 , u x2 u2 ,
residual methods on the
element level! du du
a Q1 , a Q2
dx x x1 dx x x2
Conclusion
• Good at Hand Calculations, Powerful
when applied to computers
• Only limitations are the computer
limitations
