FEM Presentation

FEM
Dr.K.Vinoth Babu
FEM
• Existing Mathematical tools not significant to find even an approximate
solution
• In FEM we get approximate solution and also improve approximate solution
by more computation.
• Find solution of a complicated problem by replacing it by a simpler one.
• Process of dividing the problem domain in to an equivalent system smaller domains or units and
selecting discrete number of points is called discretization.
• One problem domain discritized solution can be obtained for each of smaller domain. Finally
subdomain with solution can be combined together to obtain solution for entire domain. Hence this
approach is called going from point to whole.
• solution find region is build up by many small interconnected subregion called “Finite elements”.
• These Elements are interconnected at specified points called nodes or nodal points.
•Since the actual variation of the field variable (Stress,
displacement, temperature, pressure, velocities) inside the
continuum is not known.
• We can assume that variation of field variables inside a FE can
be approximated by a simple function called “ Interpolation
model” (approximate functions)
•By solving the FE equations which are generally in the form of
Matrix equation the nodal values of field variables will be
known.
Example
• Convert irregular polygon in to number of regular polygons – discretization
• Identify the formula required to calculate the area – Shape function
• Find the size parameter for each polygon based on the size given in the
diagram - Boundary condition
• Arrive the area each regular polygon by solving it – Solution
• Assemble the area of individual polygons to get the total area of polygon
Same Procedure can be followed to arrive solution by FEM. They are discussed below
General Procedure
• Divide structure in to discrete elements (Discretization)
• Selection of element types
• Select approximate function/interpolation or displacement node
• Define the gradients of the unknown Quality and constitute relationships
• Derive element equation (Stiffness matrix x force vector)
• Assemble the element equations to obtain global or Total equations
• Solve from the unknown DOF
• Solve for secondary quantities

Finite Element Method (FEM)
• A continuous function of a continuum (given

domain ) having infinite degrees of freedom is
replaced by a discrete model, approximated by a
set of piecewise continuous functions having a
finite degree of freedom.
General Example
• A bar subjected to some f

excitations like applied force at x
one end. Let the field quantity f Subdomain e
flow through the body, which has

been obtained by solving Domain divided with subdomains
with degrees of freedom
governing DE/PDE, In FEM the
domain  is subdivided into sub f1
domain and in each sub domain a
f2
piecewise continuous function is f5 f6
f3
f4
assumed. x
General Steps of the FEM
• 1. Discretize & Select the Element Types
• 2. Select a Displacement Function
• 3. Define the Strain/Displacement & Stress/Strain
• Relationships
• 4. Derive the Element Stiffness Matrix & Equations
• 5. Assemble the Element Equations to Obtain the Global
• & Introduce Boundary Conditions
• 6. Solve for the Unknown Degrees of Freedom
• 7. Solve for the Element Strains & Stresses
• 8. Interpret the Results
Discretize & Select the Element Types
• Divide the body into equivalent systems of finite elements

with nodes and the appropriate element type
• Element Types:
– One-dimensional (Line) Element
– Two-dimensional Element
– Three-dimensional Element
– Axisymmetric Element
One Dimensional Element
Select a Displacement Function
• There will be a displacement function for each element

Pascal’s Triangle
Define Strain Displacement & Stress/Strain
Relationships
• For one-dimensional;
Deformation in the x-direction,
strain e is related to the
displacement u
du d[N]
e  {d}  [B]{d}
dx dx
• [B] – Matrix relating strain to 1
nodal displacement [ B]  [ 1 1]
L
• Hooke’s Law is used for the

stress/strain relationship
sxEex To Stiffness Matrix
Derive the Element Stiffness Matrix & Equations
• Virtual work principle of a deformable body in equilibrium is

subjected to arbitrary virtual displacement satisfying compatibility
condition (admissible displacement), then the virtual work done
by external loads will be equal to virtual strain energy of internal
stresses.
• dUe is the element internal energy
• dWe is the element external energy
• Please view the integration sheet
Stiffness Matrix
• fe – Element Force
• ke – Element Stiffness Matrix
• de – Element Displacement
• E – Young Modulus
• A – Cross Section Area
• L - Length
Assemble Equations for Global Matrix &
Introduce Boundary Conditions
• Combine each element stiffness matrix into one,
which is known as the global matrix
• This is done by combining each [ke] into their
proper location on the global matrix
{F} = [K]{D}
• Capital letters represent the same as the element
stiffness matrix, but for global matrix
Solve for Unknown DOF’s
• Using the global matrix with the boundary conditions, we can

now eliminate some variables and solve for the unknowns,
i.e. displacements, end forces
Solve for Element Strains & Stresses
Interpret Results
• Solve for the stress using the equation below
• To interpret the results use the FBD with your found

values or use the computer program Algor
FEM Steps (Displacement Method)
• Discretize into finite elements, Identify nodes &

elements
• Develop element stiffness matrices [ke] for all
elements
• Assemble element stiffness matrices to get the global
stiffness matrix
• Apply kinematic boundary conditions
• Solve for displacements
• Finally solve for element forces and stresses by
picking proper rows
Example
• We combine the two stiffness matrices into the
global matrix.
• Now that the displacement at u2 has been obtained, the end forces and
stress values can be obtained by reverting back to the individual element
stiffness matrices
• For the stress, you only need to look at the individual node of the stifness
equation
Reactions
Element Forces
Element Stresses
Final Notes
• For this case, the calculated stresses in elements 1 &
2 are exact within the linear theory for 1-D bar
structures. Smaller finite elements will not help
• For tapered bars, averaged values of the cross-
sectional areas should be used for the elements.
• The displacements must be found first to find the
stresses, since we are using the displacement based
FEM
Assignment
• Write the displacement
functions for the
following elements:
• Analyze the bar shown

below for:
– (a) Displacement at B
– (b) End Forces
– (c) Average Stresses in bar
AB & BC
Applications of Finite Element Methods
 Structural & Stress Analysis

 Thermal Analysis
 Dynamic Analysis
 Acoustic Analysis
 Electro-Magnetic Analysis
 Manufacturing Processes
 Fluid Dynamics
 Financial Analysis
Applications: Aerospace Engineering
(AE)
35
Applications: Civil Engineering (CE)
36
Applications: Electrical Engineering
(EE)
37
Applications: Biomedical Engineering
(BE)
38
The Future – Virtual Engineering
39
Review of Basic Statics and Mechanics of Materials
• Static equilibrium conditions/free-body diagram
• Stress, strain and deformation
• Constitutive law – Hooke’s law
• Analysis of axially loaded bar, truss, beam and frame
• 2-D elasticity
Review of Matrix Algebra

• Matrix operation: addition, subtraction, multiplication
• Basic definitions and properties of matrix
• Inverse of matrix and solution of linear equations
• etc
1. Problem Description
• This exercise consists of an analysis of an electronics component cooling design
using fins:
• All electronic components generate heat during the course of their operation.
• To ensure optimal working of the component, the generated heat needs to be
removed. This is done by attaching fins to the device which aid in rapid heat
removal to the surroundings.
• For the sake of simplicity, we’ll assume that the electronic circuit is made of copper
with thermal conductivity of 386 W/mºK and that it generates heat at the rate of 1
W. The enclosing container is made of a steel with thermal conductivity of 17
W/mºK.
• The fins are made of aluminum with thermal conductivity of 180 W/mºK. There is
convection along all the boundaries except the bottom, which is insulated. The film
(convection) coefficient is h=50 W/m2ºK and the ambient temperature is 20ºC.

FEM Presentation

Uploaded by

Copyright:

Available Formats

FEM Presentation

Uploaded by

Document Information

Original Description:

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

FEM Presentation

Uploaded by

Copyright:

Available Formats

FEM

• Selection of element types

• Select approximate function/interpolation or displacement node

• Define the gradients of the unknown Quality and constitute relationships

• Derive element equation (Stiffness matrix x force vector)

• Assemble the element equations to obtain global or Total equations

• Solve from the unknown DOF

• Solve for secondary quantities

• A continuous function of a continuum (given

• A bar subjected to some f

one end. Let the field quantity f Subdomain e

flow through the body, which has

• Divide the body into equivalent systems of finite elements

• There will be a displacement function for each element

• Hooke’s Law is used for the

• Virtual work principle of a deformable body in equilibrium is

• Using the global matrix with the boundary conditions, we can

• Solve for the stress using the equation below

• To interpret the results use the FBD with your found

• Discretize into finite elements, Identify nodes &

• Analyze the bar shown

 Structural & Stress Analysis

Review of Matrix Algebra

You might also like