Lecture - Finite Element Method Theory
Lecture - Finite Element Method Theory
• The Finite Element Method (FEM) is a numerical method for solving engineering problems by finding
approximate solutions of the partial differential equations that represent the problem. Therefore it is a
mathematical representation of a physical system and the solution of that mathematical
representation.
• In other words, FEM is a method for dividing up a very complicated problem into small elements that
can be solved in relation to each other.
• FEM is useful for problems with complex geometries, loadings and material properties where analytical
solutions cannot be obtained.
2
How does Finite Element Method work?
• FEM uses a system of points called nodes which form a grid called a mesh.
• This mesh is programmed to contain the material and structural properties which define how much the
structure will react to certain loading conditions.
• Nodes are assigned at a certain density throughout the material depending on the anticipated stress
levels of a particular area.
• Regions which will receive large amounts of stress usually have a higher node density than those which
experience little or no stress. Points of interest may consist of: fracture point of previously tested
material, fillets, corners, complex detail and high stress areas.
• The mesh acts like a spider web in that from each node, there extends a mesh element to each of the
adjacent nodes. This web of vectors is what carries the material properties to the object, creating many
elements.
3
Types of FEM Engineering Analysis
Structural Analysis:
4
Types of FEM Engineering Analysis
Vibrational Analysis:
5
Types of FEM Engineering Analysis
Fatigue Analysis:
• Fatigue analysis helps designers to predict the life of a material or structure by showing the effects of
cyclic loading on the specimen.
• Such analysis can show the areas where crack propagation is more likely to occur.
• Failure due to fatigue may also show the damage tolerance of the material.
6
Types of FEM Engineering Analysis
Heat transfer Analysis:
• Heat transfer analysis models the conductivity or thermal fluid dynamics of the material or structure.
• This may consist of steady-state or transient transfer.
• Steady-state transfer refers to constant thermal properties in the material that yield linear heat diffusion.
7
Finite Element Analysis/Method
FEA requires three main steps:
• Pre-Processing – A finite element mesh is developed to divide the given
geometry into subdomains for mathematical analysis and the material properties
are applied and also the boundary conditions.
• Solving Matrices (solver) – In this phase, governing matrix equations are derived
and the solution for the primary quantities is generated.
9
Discretization
The body is modelled by dividing it into an equivalent system of many smaller
bodies or units (finite elements) interconnected at points common to two or more
elements (nodes or nodal points) and/or boundary lines and/or surfaces.
Example
10
FEA Pre-Processing
12
Boundary Conditions
Loads
Loads are used to represent inputs to the system.
They can be in the forms of forces, moments,
pressures, temperature, or accelerations.
Constraints
Constraints are used as reactions to the applied
loads. Constraints can resist translational or
rotational deformation induced by applied loads.
13
Boundary Conditions
Degrees of Freedom
Spatial DOFs refer to the three translational and three rotational
modes of displacement that are possible for any part in 3D
space. A constraint scheme must remove all six DOFs for the
analysis to run.
Elemental DOFs refer to the ability of each element to transmit
or react to a load. The boundary condition cannot load or
constrain a DOF that is not supported by the element to which it
is applied.
14
Boundary Conditions
A solid face should always have at least three points in
contact with the rest of the structure. A solid element
should never be constrained by less than three points and
only translational DOFs must be fixed.
Accuracy
15
Boundary Conditions
Linear Static Analysis
Boundary conditions are assumed constant from
application to final deformation of system and all loads
are applied gradually to their full magnitude.
Dynamic Analysis
The boundary conditions vary with time.
Non-linear Analysis
The orientation and distribution of the boundary
conditions vary as displacement of the structure is
calculated.
16
Boundary Conditions -Example
Excessive Constraints
Model of the chair seat with patches representing the tops of
the legs.
Patch 1
Patch 2
Patch 3 Patch 4
17
Boundary Conditions -Example
It may appear to be acceptable to constrain each circular patch
in vertical translation while leaving the rotational DOFs
unconstraint. This causes the seat to behave as if the leg-to-
seat interfaces were completely fixed (over-constrained).
A more realistic constraint scheme would be to pin the
center point of each circular patch (translational), allowing
the patch to rotate. Each point should be fixed vertically,
and horizontal constraints should be selectively applied so
that in-plane spatial rotation and rigid body translation is
removed without causing excessive constraints.
Patch 1
Patch 2
Patch 3 Patch 4
18
Boundary Conditions -Example
• Constraining the center point of patch 1 in all 3
translational DOFs.
• Constraining x and y translations of the center point of
patch 2.
• Constraining z and y translation of the center point of
patch 3.
• Constraining just the y translation of the center point of
patch 4.
This scheme allows in-
Patch 1
plane translation induced
Patch 2
by bending of the seat Patch 3 Patch 4
without rigid body
translation or rotation.
19
Summary of Pre-Processing
20
Solving the Model - Solver
Once the mesh is complete, and the properties and
boundary conditions have been applied, it is time to solve
the model. In most cases, this will be the point where the
computer does the work.
21
Post-Processing, Displacement Magnitude
Unexpectedly high or low displacements (by order of magnitude)
could be caused by an improper definition of load and/or elemental
properties.
22
Post-Processing, Displacement Animation
Animation of the model displacements serves as the best means of
visualizing the response of the model to its boundary conditions.
23
Post-Processing, FEA of a connecting rod
24
Post-Processing, Stress Results
The magnitude of the stresses should not be entirely unexpected.
25
Post-Processing, thermal analysis
26
View Animated Post-Processing
Displacements
No Review Boundary
Does the shape of deformations make sense?
Yes Conditions
View Displacement
Fringe Plot
28
Advantages of Finite Element Analysis/Method
• Can readily handle complex geometry.
• Can handle complex analysis types like vibration, heat transfer, fluids etc.
• Can handle complex loading:
1. Node-based loading (point loads)
2. Element-based loading (pressure, thermal, inertial forces)
3. Time or frequency dependent loading
• Can handle complex restraints: Indeterminate structures can analyzed
• Can handle bodies of nonhomogeneous materials: Can handle bodies comprised
of non-isotropic materials: Orthotropic and Anisotropic
• Special material effects are handled such as temperature dependent properties,
plasticity, creep, swelling etc.
29
Disadvantages of Finite Element Analysis/Method
• A specific numerical result is obtained for a specific problem.
• The FEM is applied to an approximation of the mathematical model of system.
• Experience and judgement are needed in order to construct a good finite element
model.
• A powerful computer and reliable FEM software are essential.
• Input and output data may be large and tedious to prepare and interpret.
• Numerical errors such as the limitation of the number of significant digits,
rounding-off occur very often.
• Fluid elements with boundaries at infinity can be computed and treated by using
boundary element method.
30
Software Packages for Finite Element Method
• ANSYS
• NASTRAN
• PATRAN
• NISA/DISPLAY III
• LS DYNA
• HYPERMESH
• CATIA
• Pro-E(CREO)
• SOLID WORKS
• COSMOS
31
The End
32