Module 1 Introduction To FEM
Module 1 Introduction To FEM
FiniteElementAnalysis(ENGR455)
Dr.AndreasSchiffer
AssistantProfessor,MechanicalEngineering
Tel:+971(0)24018204
[email protected]
TheFiniteElementMethod
The Finite Element Method (FEM) is a numerical method for
solving problems of engineering and mathematical physics.
In this method, the partial differential equations of a mathematical
model are discretized to obtain a set of simultaneous algebraic
equations.
The discretization is achieved by dividing the solution domain into
an equivalent system of smaller bodies or units (= finite elements).
Elements
Nodes
The elements are
interconnected at
points common to two
or more elements
(nodal points or nodes)
2
TheFiniteElementMethod
The FE solutions yields approximate values for the unknowns at
discrete points in space.
The FE method is important because for problems involving
complicated geometries, loadings and material properties, it is
generally not possible to obtain analytical mathematical solutions.
GeneralStepsofanFEanalysis
Step 1: Discretize the Problem and Select the Element Types
Divide the structure into small pieces, usually called meshing
Type of element depends on the nature of the problem
(structural) and the dimension (1D, 2D or 3D)
GeneralStepsofanFEanalysis
Step 2: Select a Displacement Function
Choosing a displacement function within each element
connecting the nodes.
Linear, quadratic and cubic polynomials are most common.
Step 3: Define the Stress-Strain Relationships
Describe the constitutive law relating stresses to strains in each
element
The simplest constitutive relation is Hookes law,
x = E x (in 1D).
Step 4: Derive the Element Stiffness Matrix and Equations
Write the system of equations describing the structural behavior
of an element.
Relating nodal forces to nodal displacements of the element
through an element stiffness matrix: f k d
5
GeneralStepsofanFEanalysis
Step 5: Assemble the Global System Equation
All individual elements are assembled using the method of
superposition (or direct stiffness method) to produce the global or
total system of equations of the problem.
F K d
Here, {F} is the vector of global nodal forces, [K] is the structure
global or total stiffness matrix, {d} is now the vector of known and
unknown structure nodal degrees of freedom (displacements).
It can be shown that at this stage, the global stiffness matrix [K] is
a singular square matrix because its determinant is equal to zero.
To remove this singularity, we must invoke certain boundary
conditions (or constraints or supports) so that the structure
remains in place instead of moving as a rigid body.
Step 6: Apply Boundary Conditions and Loading
Prescribe forces and displacements at nodes.
F K d
6
GeneralStepsofanFEanalysis
Step 7: Solve for the Unknown Degrees of Freedom
Involves finding the inverse of the global stiffness matrix [K]-1 .
Then the structures unknown nodal degrees of freedom {d} can
be calculated via
1
F K d
d K F
ApplicationoftheFEmethod
Impact analysis of an ice deflector ramp for the
railway industry
Experimental investigation
ApplicationoftheFEmethod
Impact analysis of an ice deflector ramp for the
railway industry
Plastic strains induced in the ramp
ApplicationoftheFEmethod
Side-crash simulation of a road car
v0
Modellingofdamageand
failurebyelementdeletion
DeformationoftheBframe
aftertheimpact
10
ApplicationoftheFEmethod
Deformation of a sandwich plate under blast loading
(Source: ABAQUS)
TNT
after1kgTNT:
after2kgTNT:
Deformationofthehoneycombcoreafter1kgTNT
11
ApplicationoftheFEmethod
Structural dynamics analysis of a steel bridge
12
ApplicationoftheFEmethod
Fluid-structure interaction in underwater blast loading
Freestandingrigid
plate
13
ApplicationoftheFEmethod
Crack propagation in a random heterogeneous material
FEsimulations:
Experiments:
Crackpropagation
ismodeledbyusing
theextendedFE
method(XFEM)
14