Finite Element Method
Finite Element Method
"Finite element" redirects here. For the elements of a poset, see compact element.
Visualization of how a car deforms in an asymmetrical crash using finite element analysis
Differential equations
Scope
show
Fields
Classification
show
Types
show
Relation to processes
Solution
show
General topics
show
Solution methods
People
show
List
v
t
e
Basic concepts[edit]
The subdivision of a whole domain into simpler parts has several advantages: [2]
1. dividing the domain of the problem into a collection of subdomains, with each
subdomain represented by a set of element equations to the original problem
2. systematically recombining all sets of element equations into a global system of
equations for the final calculation.
The global system of equations has known solution techniques, and can be calculated from
the initial values of the original problem to obtain a numerical answer.
In the first step above, the element equations are simple equations that locally approximate the
original complex equations to be studied, where the original equations are often partial differential
equations (PDE). To explain the approximation in this process, the finite element method is
commonly introduced as a special case of Galerkin method. The process, in mathematical
language, is to construct an integral of the inner product of the residual and the weight
functions and set the integral to zero. In simple terms, it is a procedure that minimizes the error of
approximation by fitting trial functions into the PDE. The residual is the error caused by the trial
functions, and the weight functions are polynomial approximation functions that project the
residual. The process eliminates all the spatial derivatives from the PDE, thus approximating the
PDE locally with
History[edit]
While it is difficult to quote a date of the invention of the finite element method, the method
originated from the need to solve complex elasticity and structural analysis problems
in civil and aeronautical engineering. Its development can be traced back to the work by A.
Hrennikoff[3] and R. Courant[4] in the early 1940s. Another pioneer was Ioannis Argyris. In the
USSR, the introduction of the practical application of the method is usually connected with name
of Leonard Oganesyan.[5] It was also independently rediscovered in China by Feng Kang in the
later 1950s and early 1960s, based on the computations of dam constructions, where it was
called the finite difference method based on variation principle. Although the approaches used by
these pioneers are different, they share one essential characteristic: mesh discretization of a
continuous domain into a set of discrete sub-domains, usually called elements.
Hrennikoff's work discretizes the domain by using a lattice analogy, while Courant's approach
divides the domain into finite triangular subregions to solve second order elliptic partial differential
equations that arise from the problem of torsion of a cylinder. Courant's contribution was
evolutionary, drawing on a large body of earlier results for PDEs developed by Rayleigh, Ritz,
and Galerkin.
The finite element method obtained its real impetus in the 1960s and 1970s by the developments
of J. H. Argyris with co-workers at the University of Stuttgart, R. W. Clough with co-workers at UC
Berkeley, O. C. Zienkiewicz with co-workers Ernest Hinton, Bruce Irons[6] and others at Swansea
University, Philippe G. Ciarlet at the University of Paris 6 and Richard Gallagher with co-workers
at Cornell University. Further impetus was provided in these years by available open source finite
element programs. NASA sponsored the original version of NASTRAN, and UC Berkeley made
the finite element program SAP IV[7] widely available. In Norway the ship classification society Det
Norske Veritas (now DNV GL) developed Sesam in 1969 for use in analysis of ships.[8] A rigorous
mathematical basis to the finite element method was provided in 1973 with the publication
by Strang and Fix.[9] The method has since been generalized for the numerical modeling of
physical systems in a wide variety of engineering disciplines, e.g., electromagnetism, heat
transfer, and fluid dynamics.[10][11]
History[edit]
While it is difficult to quote a date of the invention of the finite element method, the method
originated from the need to solve complex elasticity and structural analysis problems
in civil and aeronautical engineering. Its development can be traced back to the work by A.
Hrennikoff[3] and R. Courant[4] in the early 1940s. Another pioneer was Ioannis Argyris. In the
USSR, the introduction of the practical application of the method is usually connected with name
of Leonard Oganesyan.[5] It was also independently rediscovered in China by Feng Kang in the
later 1950s and early 1960s, based on the computations of dam constructions, where it was
called the finite difference method based on variation principle. Although the approaches used by
these pioneers are different, they share one essential characteristic: mesh discretization of a
continuous domain into a set of discrete sub-domains, usually called elements.
Hrennikoff's work discretizes the domain by using a lattice analogy, while Courant's approach
divides the domain into finite triangular subregions to solve second order elliptic partial differential
equations that arise from the problem of torsion of a cylinder. Courant's contribution was
evolutionary, drawing on a large body of earlier results for PDEs developed by Rayleigh, Ritz,
and Galerkin.
The finite element method obtained its real impetus in the 1960s and 1970s by the developments
of J. H. Argyris with co-workers at the University of Stuttgart, R. W. Clough with co-workers at UC
Berkeley, O. C. Zienkiewicz with co-workers Ernest Hinton, Bruce Irons[6] and others at Swansea
University, Philippe G. Ciarlet at the University of Paris 6 and Richard Gallagher with co-workers
at Cornell University. Further impetus was provided in these years by available open source finite
element programs. NASA sponsored the original version of NASTRAN, and UC Berkeley made
the finite element program SAP IV[7] widely available. In Norway the ship classification society Det
Norske Veritas (now DNV GL) developed Sesam in 1969 for use in analysis of ships.[8] A rigorous
mathematical basis to the finite element method was provided in 1973 with the publication
by Strang and Fix.[9] The method has since been generalized for the numerical modeling of
physical systems in a wide variety of engineering disciplines, e.g., electromagnetism, heat
transfer, and fluid dynamics.[10][11]
Technical discussion[edit]
The structure of finite element methods[edit]
A finite element method is characterized by a variational formulation, a discretization strategy,
one or more solution algorithms, and post-processing procedures.
Examples of the variational formulation are the Galerkin method, the discontinuous Galerkin
method, mixed methods, etc.
A discretization strategy is understood to mean a clearly defined set of procedures that cover (a)
the creation of finite element meshes, (b) the definition of basis function on reference elements
(also called shape functions) and (c) the mapping of reference elements onto the elements of the
mesh. Examples of discretization strategies are the h-version, p-version, hp-version, x-
FEM, isogeometric analysis, etc. Each discretization strategy has certain advantages and
disadvantages. A reasonable criterion in selecting a discretization strategy is to realize nearly
optimal performance for the broadest set of mathematical models in a particular model class.
Various numerical solution algorithms can be classified into two broad categories; direct and
iterative solvers. These algorithms are designed to exploit the sparsity of matrices that depend on
the choices of variational formulation and discretization strategy.
Postprocessing procedures are designed for the extraction of the data of interest from a finite
element solution. In order to meet the requirements of solution verification, postprocessors need
to provide for a posteriori error estimation in terms of the quantities of interest. When the errors of
approximation are larger than what is considered acceptable then the discretization has to be
changed either by an automated adaptive process or by the action of the analyst. There are some
very efficient postprocessors that provide for the realization of superconvergence.