MEEN 673
*********************************
NONLINEAR FINITE ELEMENT
ANALYSIS
*********************************
J. N. Reddy
e-mail: jnreddy@tamu.edu
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 )
i1
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 )
i1
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