MEEN 673

*********************************
NONLINEAR FINITE ELEMENT
ANALYSIS
*********************************

J. N. Reddy
e-mail: [email protected]
Texas A&M University
College Station, TX 77843-3123 USA

JN Reddy
 INTRODUCTORY REMARKS

 What we do as engineers?
 Role of numerical simulations
 Classical variational methods
and their major drawback
 The main features of FEM
 Terminology used

JN Reddy Introduction: 2
 GENERAL INTRODUCTION
(MOTIVATION FOR THE COURSE)

Engineering is the discipline, art, and

profession of acquiring and applying
technical, scientific, and mathematical
knowledge to design and implement materials,
structures, machines, devices, systems, and
processes that safely realize a desired objective.

JN Reddy Introduction: 3
 GENERAL INTRODUCTION
(WHAT IS ENGINEERING?)

Engineering is a problem-solving discipline,

and solution requires an understanding
of the phenomena that occurs in the system.

The study of natural phenomena involves

 developing mathematical models,
 conducting physical experiments,
 carrying out numerical simulations, and
 designing and building systems to achieve
a functionality in most economical way.

JN Reddy Introduction: 4
 GENERAL INTRODUCTION
ROLE OF ANALYSIS IN DESIGN AND
MANUFACTURING

Analysis is an aid to design and manufacturing,

and not an end in itself.
Analysis steps:
 identifying the problem and nature of the
response to be determined,
 selecting the mathematical model,
 selecting a solution method, and
 evaluating the results in light of the
design parameters.
JN Reddy Introduction: 5
 GENERAL INTRODUCTION
BENEFITS OF NUMERICAL SIMULATIONS
 Mathematical Model Development and
Computer Simulations continues to be a
major component of engineering analysis,
design and manufacturing (CAE and CAM).

Computer (numerical) simulations

reduce/replace prototype testing and hence
reduce time and product costs

Numerical simulations facilitate

investigations into the use of alternative
materials and configurations.
JN Reddy Introduction: 6
 GENERAL INTRODUCTION
WHY SHOULD YOU TAKE THE COURSE?
 A good understanding of the phenomena
modeled and the computer method used to
simulate the process is essential for the
analyst to aid the development and
manufacturing of complex systems.
Engineer-scientists with good background in
engineering and science subjects as well as in
computer modeling techniques will continue
to have excellent opportunities to contribute
to the science and technology for the benefit
of the mankind.
JN Reddy Introduction: 7
 NUMERICAL SIMULATION OF A
TYPICAL PHYSICAL PROCESS

Physical
System
Assumptions Laws of physics
concerning (conservation
the system principles)
Mathematical
Model
(BVP, IVP)
BVP – Boundary value problems FEM – Finite Element
(equilibrium problems) Method
IVP – Initial value problems FDM  Finite Difference
Numerical
(time-dependent problems) Method
Simulations
BEM – Boundary Element
Method
Numerical method Computational
(FEM,FDM,BEM,etc.)
device

JN Reddy Introduction: 8
 BASIC FEATURES OF THE
FINITE ELEMENT METHOD (FEM)

 Divide whole into parts (finite element mesh)

 Set up the `problem’ over a typical part
(derive a set of relationships between
primary and secondary variables)
 Assemble the parts to obtain the solution to
the whole

JN Reddy Basic Concepts: 9

 A BOUNDARY VALUE PROBLEM (BVP)
  u    u 
 a11   a22   f  0 in 
x  x  y  y 
 u   u 

u  uˆ or a11  
 nx  a22  n  qˆ on 
 y 
y n
 x 

y  nyˆj
n̂
nx ˆi

ˆ  nx ˆi  ny ˆj
n
x
JN Reddy Introduction: 10
 Exact and Approximate Solutions
An exact solution satisfies (a) the differential equation at
every point of the domain and (b) boundary conditions on
the boundary. An approximate solution satisfies the
differential equation as well as the boundary conditions in
some “acceptable sense” (to be made clearer shortly).
We seek the approximate solution as a linear combination of
unknown parameters ci and known functions i ( x , y ) that
Satisfy the boundary conditions:
N
u( x , y )  uN ( x , y )   ci i ( x , y )
i1

We determine uN ( x , y ) such that the above expression

satisfies the differential equation in a weighted-residual
sense.
JN Reddy Basic Concepts: 11
 Determining Approximate Solutions
Suppose that i is selected to satisfy the boundary conditions
exactly. Then substitution of uN(x,y) into the differential
equation

  uN    uN 
 a11   a22   f  R(x , y, ci )  0 in 
x  x  y  y 

Then ci are determined such that the residual (or error in the
differential equation), R(x,y,ci), is zero in the weighted-
residual sense:
0   w R dxdy, i  1,2, , N
i


where w are linearly independent set of weight functions

i

JN Reddy Basic Concepts: 12

 WEIGHTED-INTEGRAL METHODS
for the Numerical Solution of Differential Eqs.

0   w R( x , y, c , c ,..., c ) dxdy, i  1,2, , N

i 1 2 N


Collocation method:
0  R( xi , yi )  wi ( x , y )   ( x  xi , y  yi )  Dirac delta

Least-squares method:
R R
0 R dxdy  w ( x , y ) 
c c
i

i i

Galerkin Method:
0    R dxdy  w ( x , y )   ( x , y )
i i i


JN Reddy Basic Concepts: 13

 MAJOR DRAWBACK OF CONVENTIONAL
VARIATIONAL METHODS
y
  u    u 
 a11   a22   f  0 in 
x  x  y  y 

x u  0 on 
2b


2a

 
(the four line
(interior) segments)

JN Reddy Introduction: 14
 MAJOR DRAWBACK OF CONVENTIONAL
VARIATIONAL METHODS
y N
u( x , y )  uN ( x , y )   ci i ( x , y )
i1

  1 ( x , y )  ( a 2  x )(b2  y2 )
2

x
2b
or

x y
2a  1 ( x , y )  cos cos
y 2a 2b

 1 ( x , y)  ?
x
JN Reddy Introduction: 15
 Finite Element Discretization

Elements
Nodes
•
• •
•
• •
Domain, Ω •
• Ωe • •
• • • •
•
°°°° •
° •
Boundary, Γ •

(a) Given domain (b) Finite element mesh

Ωe Domain,
Boundary flux Ωh

(c) Typical element with (d) Discretized domain

boundary fluxes
JN Reddy
Basic Concepts: 16
 Approximation of a curved surface
with a plane
u(x,y) − Temperature profile

Domain (Triangular element)

JN Reddy Basic Concepts: 17

 FEM Terminology

 Element A geometric sub-domain of the region being

simulated, with the property that it allows a unique (1)
representation of its geometry and (2) derivation of the
approximation (interpolation) functions.

 Node A geometric location in the element which plays a

role in the derivation of the interpolation functions and it
is the point at which solution is sought.

 Mesh A collection of elements (or nodes) that replaces

the actual domain.

JN Reddy Basic Concepts: 18

 FEM Terminology (continued)

• Finite Element Model A set of algebraic equations

relating the nodal values of the primary variables (e.g.,
displacements) to the nodal values of the secondary
variables (e.g., forces) in an element.
 Finite element model is NOT the same as the finite
element method. There is only one finite element
method but there can be more than one finite element
model of a problem (or mathematical model).
 Numerical Simulation Evaluation of the
mathematical model (i.e., solution of the governing
equations) using a numerical method and computer.

JN Reddy Basic Concepts: 19

 Major Steps of Finite Element
Model Development
 Begin with the governing equations of the problem

 Develop its weak form over a typical finite element

 Approximate the solution over each finite element

 Obtain algebraic relations among the quantities of

interest over each finite element (i.e., finite element
model)

JN Reddy Basic Concepts: 20