- Home
- Documents
*Introduction to Finite Element Analysis (FEA) or Finite ... mech410/old/2_Lecture_Notes/6_2_FEA_...*

If you can't read please download the document

View

224Download

2

Embed Size (px)

Introduction to Finite Element Analysis y(FEA) or Finite Element Method (FEM)

Finite Element Analysis (FEA) or FiniteFinite Element Analysis (FEA) or Finite Element Method (FEM)

The Finite Element Analysis (FEA) is a numerical method for solving problems of engineering and mathematical physicsengineering and mathematical physics.

Useful for problems with complicatedUseful for problems with complicated geometries, loadings, and material propertieswhere analytical solutions can not be obtained.

The Purpose of FEAIn Mechanics Courses Analytical Solution Stress analysis for trusses, beams, and other simple

structures are carried out based on dramatic simplificationpand idealization: mass concentrated at the center of gravity beam simplified as a line segment (same cross-section)

D i i b d th l l ti lt f th id li d Design is based on the calculation results of the idealizedstructure & a large safety factor (1.5-3) given by experience.

In Engineering Design - FEAIn Engineering Design - FEA Design geometry is a lot more complex; and the accuracy

requirement is a lot higher. We need To understand the physical behaviors of a complexy

object (strength, heat transfer capability, fluid flow, etc.) To predict the performance and behavior of the design;

to calculate the safety margin; and to identify theweakness of the design accurately; andweakness of the design accurately; and

To identify the optimal design with confidence

Brief History

Grew out of aerospace industryPost-WW II jets, missiles, space flightNeed for light weight structuresRequired accurate stress analysisParalleled growth of computers

Developments1940s - Hrennikoff [1941] - Lattice of 1D bars,

- McHenry [1943] - Model 3D solids, y [ ]- Courant [1943] - Variational form, - Levy [1947, 1953] - Flexibility & Stiffness

1950 60 A i d K l [1954] E P i i l f1950-60s - Argryis and Kelsey [1954] - Energy Principle for Matrix Methods, Turner, Clough, Martin and Topp [1956] -2D elements, Clough [1960] - Term Finite Elements

1980s Wide applications due to:

Integration of CAD/CAE automated mesh generation and graphical display of analysis resultsand graphical display of analysis results

Powerful and low cost computers

2000s FEA in CAD; Design Optimization in FEA; Nonlinear2000s FEA in CAD; Design Optimization in FEA; Nonlinear FEA; Better CAD/CAE Integration

FEA ApplicationsFEA Applications

Mechanical/Aerospace/Civil/Automotive EngineeringStructural/Stress Analysis

Static/DynamicStatic/DynamicLinear/Nonlinear

Fluid FlowHeat TransferElectromagnetic FieldsSoil MechanicsAcousticsBi h iBiomechanics

DiscretizationDiscretization

Complex Object Simple Analysis(Material discontinuity,

Complex and arbitrary geometry)

Real Word

Simplified(Idealized) Physical

Mathematical Model

Discretized(mesh) Modely

ModelModel

Discretizations

Model body by dividing it into an equivalent system of many smaller bodies

it (fi it l t ) i t t d tor units (finite elements) interconnected at points common to two or more elements (nodes or nodal points) and/or boundary(nodes or nodal points) and/or boundarylines and/or surfaces.

Elements & Nodes - Nodal Quantity

Feature

Obtain a set of algebraic equations to solve for unknown (first) nodal quantity(di l t)(displacement).

S (Secondary quantities (stresses and strains) are expressed in terms of nodal values of primary quantityvalues of primary quantity

Object

Elements

Nodes

Displacement Stress

Strain

E l f FEA 2DExamples of FEA - 2D

E l f FEA 3DExamples of FEA 3D

E l f FEAExamples of FEA

Advantages

Irregular BoundariesGeneral LoadsDifferent MaterialsDifferent MaterialsBoundary ConditionsVariable Element SizeVariable Element SizeEasy ModificationDynamicsNonlinear Problems (Geometric or Material)

E l f FEAExamples for FEA

FEA in Fuel Cell Stack Designg

H2Air (O2) 2( 2)

A PEM Fuel Cell

Compression of A PEM Fuel Cell StackCooling

PlateCooling

Plate

Gas/OxidantDelivery Plate

MEA MEA

O 2 H 2 H 2O 2H 2O H 2O

e -

e -

e -

e -

e -

e -

StackPressure

StackPressure

e ee

Heat

Load

FEA on the Fuel Cell Stack and End Plate

Fuel Cell in the First Ballard Prototype BusFuel Cell in the First Ballard Prototype Bus

One of the key design challenge - getting rid of the low grid y g g g g gheat using an inefficient stainless radiator

A Novel Fuel Cell Stack DesignTri-stream, External-manifolding, Radiator Stack

(TERS)

+Obj ecti ves:

Adv ant age s:

A new PE MFC sta ck arch itec ture for tr ansp orta tion ap plic atio ns w ith low cos ts a nd h igh syst em pow er den sity

Lo w c om pon ent , ma nuf actu ring , se alin g,

a nd o pe ratio n c ostsC om pac t de sign wit h in tegr ated en gin e a nd r adia tor; hig h sy stem po wd er d ens ity M ulti- func tion al p ane ls m ad e ex tern al m ani fold ing feas ible for a PE M s tac k

Ambient Air Out

H2 InAir In

Ambient Air Out

Air Out

H2 Out

Ambient Air In Multi-functional Panel

The Multi-functional Panel

Heat transfer and rejection Deformation: compensation to thermal and

hydro expansionhydro expansion Electrical conductivity

Design Objective: Ideal Compression Force and Deformation - StiffnessDesign Objective: Ideal Compression Force and Deformation - Stiffness

- cannot be achieved without modern design tool: FEA & Optimization

Stiffness Analysis and Design Optimization of the Panel

Design Variables: Shape

H i ht Height Wavelength Thickness

S f fi i h Surface finish Cuts

Principles of FEAPrinciples of FEA

e

StiffnessStiffness Equation

(first nodal quantity)(first nodal quantity)

(second nodal quantities)

A Si l Stiff E tiA Simple Stiffness Equation

Kx F=

KF

x

Simplest

Stiff E ti f O S iStiffness Equation of One Springxi xjK

ifi fj

j

ii ij i ik k x ffk k

=

(( )

)i j iK x x f

K f =

j jji jj x fk k

By using the unit displacement method, we can express the

( )i j jxK x f =

By using the unit displacement method, we can express the stiffness coefficients kij etc. in terms of the spring coefficient K

KkkKkk ==== KkkKkk jiijjjii ==== ,

Th Th S i S tThe Three-Spring System

x x xx1

F1

x2

F2

x3

F3K3K2K1

No geometry influence

1 2 3

No geometry influence Simple material property (K) Simple load condition Simple constraints

(boundary condition) ixi

f

xj

fj

K

fi fj

Identify and Solve the Stiffness EquationsIdentify and Solve the Stiffness Equationsfor a System of Finite Elements

x1 x2 x3

KKK F1 F2 F3K3K2K1

11131211 Fxkkk

=

22232221

F

Fxx

kkkkkk

33333231 Fxkkk Medium

An Elastic Solid > A System of SpringsAn Elastic Solid --> A System of Springs

Task of FEA: To identify and solve the stiffness equations for a system of finite elements.equations for a system of finite elements.

Real-Complex

An Elastic Solid A System of SpringsAn Elastic Solid A System of Springs

The actual solid (plate, shell, etc.) is discretized into a number of smaller units called elements.

These small units have finite dimensions - hence the d fi it l tword finite element.

The discrete equivalent spring system provides an approximate model for the actual elastic bodyapproximate model for the actual elastic body

It is reasonable to say that the larger the number of elements used the better will be the approximationelements used, the better will be the approximation

Think spring as one type of elastic units; we can use other types such as truss, beam, shell, etc.other types such as truss, beam, shell, etc.

Real-Complex

A System of Springs under A Number y p gof Forces

The systems configuration will change.

We need to measure deflections at several points to characterize such changes.

A system of linear equations is introduced.

11 12 1 1 1nk k kk k k

FF

xx

L

L21 22 2 22nk k k

k k

F

k F

x

x

=

M M M

L

MM M

1 2n nn nn nk k k Fx Real-Complex

Procedure for Carrying out Finite y gElement Analysis

To construct the stiffness equations of a complex system made up of springs, one need to develop the stiffness equation of one spring and use the equation as a building block Stiffness equation of one spring/block

Way of stacking blocks

3 nodes2 elements

3 nodes2 elementse e e s

W f St ki Bl k /El tWay of Stacking Blocks/Elements

Compatibility requirement: ensures that the displacements at the shared node of adjacent elements are equal.elements are equal.

Equilibrium requirement: ensures that elemental forces and the external forces applied to the

t d i ilib isystem nodes are in equilibrium. Boundary conditions: ensures the system satisfy

the boundary constraints and so on.the boundary constraints and so on.

Applying Differen