finite element method (fem): an overview
TRANSCRIPT
![Page 1: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/1.jpg)
FINITE ELEMENT METHOD (FEM): AN
OVERVIEW
Dr A Chawla
![Page 2: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/2.jpg)
ANALYTICAL / MATHEMATICAL SOLUTIONS
• RESULTS AT INFINITE LOCATIONS
• CONTINUOUS SOLUTIONS • FOR SIMPLIFIED SITUATIONS ONLY • EXACT SOLUTION
![Page 3: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/3.jpg)
NUMERICAL (FEM) SOLUTIONS
• APPROXIMATE SOLUTIONS
• VALUES AT DISCRETE LOCATIONS • FOR COMPLEX GEOMETRY MATERIAL PROPERTIES LOADING BOUNDARY CONDITIONS
![Page 4: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/4.jpg)
THE FINITE ELEMENT METHOD
• A METHOD OF PIECEWISE APPROXIMATION
• BY CONNECTING SIMPLE FUNCTIONS • EACH VALID OVER A SMALL REGION / ELEMENT • A PROCESS OF DISCRETIZATION
![Page 5: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/5.jpg)
ESSENTIAL STEPS IN FEM
• DISCRETIZATION
• SELECTION OF THE DISPLACEMENT MODELS
• DERIVING ELEMENT STIFFNESS MATRICES
• ASSEMBLY OF OVERALL EQUATIONS / MATRICES
• SOLUTIONS FOR UNKNOWN DISPLACEMENTS
• COMPUTATIONS FOR THE STRAINS / STRESSES
![Page 6: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/6.jpg)
DISCRETIZATION (Fig 1.1)
• SELECTING CERTAIN DISCRETE POINTS (NODES)
• FORMATION OF ELEMENT MESH 2D: 3/6 NODED TRIANGLES, QUADRILATERALS 3D: TETRAHEDRAL, PRISMATIC etc • ELEMENTS INTERCONNECTED AT THE NODES • DECIDE NUMBER, SIZE AND TYPE OF ELEMENT
![Page 7: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/7.jpg)
DISPLACEMENT MODELS (Fig 1.2)
• IF NODAL DISPLACEMENTS ARE KNOWN
• DISPLACEMENT WITHIN IS COMPUTED
• USING SIMPLE FUNCTIONS (eg. POLYNOMIAL) • INTRODUCES APPROXIMATION • MODEL SHOULD SATISFY CERTAIN BASIC REQUIREMENTS TO MINIMIZE ERRORS
![Page 8: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/8.jpg)
DERIVATION OF THE ELEMENT MATRICES
• EQUIVALENT FORCES AT THE NODES
• SPECIFY MATERIAL AND GEOMETRIC PROPERTIES • STIFFNESS RELATES NODAL DISPLACEMENT TO FORCES • DERIVE STIFFNESS MATRIX • (MATRIX OF INFLUENCE COEFFICIENTS)
![Page 9: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/9.jpg)
DERIVATION OF OVERALL EQUATIONS / MATRICES
• DISPLACEMENT AT A NODE TO BE SAME
FOR ALL ADJACENT ELEMENTS
• COMBINE ELEMENT MATRICES • DERIVE EXPRESSIONS FOR POTENTIAL ENERGY • ∏ = 1/2 QT K Q - QT F
![Page 10: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/10.jpg)
SOLUTIONS FOR UNKNOWN DISPLACEMENTS
• SPECIFY BOUNDARY CONDITIONS
• USE MINIMIZATION OF P.E. (say) • DERIVE SIMULTANEOUS EQUATIONS • KQ = F (Q’s ARE UNKNOWNS) • SOLVE USING NUMERICAL TECHNIQUES
1. FOR LINEAR PROBLEMS: MATRIX AGEBRA TECHNIQUES
2. FOR NON LINEAR PROBLEMS: MODIFY STIFFNESS / FORCE MATRIX AT EACH ITERATION
![Page 11: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/11.jpg)
COMPUTE STRESSES AND STRAINS
• DERIVE STRAINS FROM DISPLACEMENTS
• DERIVE STRESSES FROM STRAINS • USING SOLID MECHANICS PRINCIPLES
![Page 12: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/12.jpg)
FUNDAMENTALS OF MECHANICS (1D)
• Stress Strain Relations
ε = du / dx σx = E ε = E du / dx
• Force Equilibrium dσx / dx + f = 0 E d2u / dx2 + f = 0
SECOND ORDER DE TO BE SOLVED
![Page 13: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/13.jpg)
BOUNDARY CONDITIONS
u = 0 at x = 0 and u = 0 at x = L
![Page 14: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/14.jpg)
FOR BENDING PROBLEMS • EQUILIBRIUM EQUATION
d2M / dx2 + q = 0 ε = z d2w / dx2 M = σ I / y
• FOURTH ORDER DE
![Page 15: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/15.jpg)
BOUNDARY CONDITIONS (in
bending) • w, dw/dx, d2w / dx2 or d3w / dx3 AT THE
BOUNDARY • for instance w(0) = 0, dw / dx (0) = 0
![Page 16: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/16.jpg)
A GENERAL 3D CASE • DEFORMATIONS
u = [u v w]T • STRESSES
σ = [σx σy σz τyz τxz τxy]T • STRAINS
ε = [εx εy εz γyz γxz γxy]T = [δu/δx δv/δy δw/δz (δv/δz+δw/δy) ...]T
• FORCES
BODY FORCES [fx fy fz]T TRACTIVE FORCES [Tx Ty Tz]T POINT FORCES [Px Py Pz]T
![Page 17: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/17.jpg)
3D EQUILIBRIUM EQUATIONS
• BODY FORCES (equilibrium of a volume
element) δσx/δx + δτxy/δy + δτxz/δz + fx = 0 δτxy/δx + δσy/δy + δτyz/δz + fy = 0 δτxz/δx + δτyz/δy + δσz/δz + fz = 0
![Page 18: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/18.jpg)
• TRACTIVE FORCES σx nx + τxy ny + τxz nz = Tx τxy nx + σy ny + τyz nz = Ty τxz nx + τyz ny + σz nz = Tz
where [nx ny nz]T : surface normal
![Page 19: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/19.jpg)
MATERIAL BEHAVIOR
• LINEAR ISOTROPIC MATERIAL (σ - ε relation defined using two constants)
εx = (σx - ν σy - ν σz) / E • ORHOTROPIC (composites)
different properties in different directions upto nine constants to relate σ - ε For instance, composite materials
• OTHER MATERIALS non-linear isotropic (rubber) hypoelastic (incremental σ - ε relation) (geological materials) elasto-plastic (-do- with plasticity)
• ONLY σ - ε relation changes • FEM REMAINS SAME
![Page 20: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/20.jpg)
MINIMUM PE PRINCIPLE • BASIS OF FEM
• ∏ = 1/2 ∫σT ε dV - ∫uTfdV - ∫uTTdS - ∑ui
TPi
• AT EQUILIBRIUM ∏ IS A MINIMA
• FOR AN ASSUMED DISPLACEMENT FIELD
• δ∏ / δai = 0
![Page 21: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/21.jpg)
ERRORS IN FEM • WRONG ASSUMTIONS
• USER ERRORS • INAPPROPRIATE ELEMENT TYPE • DISCRETIZATION ERRORS • WRONG MESH SIZE • YIELDING / BUCKLING OVERLOOKED • WRONG SUPPORT CONDITIONS • LARGE VARIATIONS IN STIFFNESSES • PROGRAM BUGS + ROUNDING OFF • IMPROPER TRAINING WITH SOFTWARE
![Page 22: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/22.jpg)
SOME POSSIBLE ANALYSIS TYPES
• STATIC ANALYSIS
• DYNAMIC (MODAL / TRANSIENT)
• THERMAL / COMBINED STRESSES
• IMPACT STRESSES
• NON-LINEAR / PLASTIC MATERIALS
• COMPOSITE MATERIALS
• COMPLICATED LOADINGS AND
BOUNDARY CONDITIONS
![Page 23: FINITE ELEMENT METHOD (FEM): AN OVERVIEW](https://reader033.vdocument.in/reader033/viewer/2022051008/58a024521a28ab954a8c4645/html5/thumbnails/23.jpg)
TYPES OF APPLICATION AREAS
• STRUCTURAL ENGINEERING APPLICATIONS
• HEAVY ENGINEERING COMPONENTS • AUTOMOBILE PARTS • AEROSPACE ENGINEERING • NUCLEAR ENGINEERING • TURBINE BLADES / OTHER POWER PLANT COMPONENTS • AND MANY MORE